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Abstract 

In this thesis, we have tried to put forth some of the aspects of light-front (LF) field 
theory through their successful application in the Deep Inelastic Scattering (DIS). We have 

q\ ■ developed a LFQCD Hamiltonian description of the DIS structure functions starting from 

Bj or ken- Johnson-Low limit of virtual forward Compton scattering amplitude and using LF 
current commutators. We worked in the LF gauge A + = and used the old-fashioned 
LFQCD perturbation theory in our calculations. The importance of our work can be sum- 

O ■ marized from two different viewpoints. 

f-«« \ Firstly, from DIS point of view, our approach is convenient for it closely follows the 

physical intuitions of the parton model. Our approach addresses directly the structure func- 
tions, which are experimental objects, instead of its moments which come naturally in the 

^vJD . usual way (i.e., using QCD improved parton model or OPE method). More importantly, our 

approach has the potential of incorporating the non-perturbative contents of the structure 

Q ' functions. We have shown introducing a new factorization scheme that the non-perturbative 

contents of the DIS structure functions can be obtained by solving the LF bound state equa- 
tions (which seems viable due to the ongoing research activity towards this direction using 
the similarity renormalization scheme), while the perturbative contents can be extracted 
Q_i| by calculating the dressed parton structure functions as we worked them out explicitly to 

q; the leading order. In contrast, in the usual approaches nonperturbative information is only 

parametrized, putting the emphasis only on the (Revolution of the moments of the structure 
functions using perturbation theory. 

Simplicity of our approach becomes evident when we try to describe structure functions 

^ ■ in the context of the nucleonic helicity by defining new structure functions. The ambiguity 

of gauge invariance and interaction dependence in defining various parts of the helicity 
operator for quarks and gluons in the usual way are absent in our gauge fixed theory and 
the well known LF helicity operator seems to provide the consistent physical information as 
we have shown explicitly. Then, we proceed to calculate anomalous dimensions relevant for 
the (Revolution of such structure functions which agree with recent calculation using other 
methods. 

Secondly, our study is important in view of establishing LF field theory as a valid frame- 
work. Even though LF field theory is a self-consistent theory, there is still some doubt 
whether it is equivalent to the more familiar equal-time one, since the formulation and the 
methods used in two cases are so different, which we discussed in detail. Perturbative cal- 
culations of the dressed parton structure functions presented in this work show that, to the 
leading order in coupling, our approach yields equivalent results as in the usual way and 
hence the equivalence. Also, our investigation in the context of coupling constant renormal- 
ization in LFQCD hamiltonian methods shows the importance of Gallilean boost symmetry 
in understanding the correctness of any higher order calculation using (x + )-ordered LFQCD 
perturbation theory. 
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Chapter 1 



I. PROLOGUE 



The interest in the high energy community has become very much diverse and, at the 
same time, very intense. It is so much so, that it becomes very difficult sometimes for one 
working in a particular field to find a connection between his or her work with that of another 
one working in a different field, let alone understanding that. This is particularly true for 
one relatively new in the field, to whom it may appear that a lot of research activities are 
going on outside his or her future light cone (in a figurative sense, of course) and hence, 
is unable to find the causal connection. In spite of all these ever-growing diversity and 
intricacy in modern day research, some of the very basic questions that intrigued our mind 
have not changed over the years. Namely, what are the basic building blocks of matter and 
how they conspire to build up the nature as we see it? As a matter of fact, the answer to 
these questions did change from time to time. 

As far as the basic building blocks are concerned, all forms of matter were thought 
to be built out of what is known as chemical elements in the early days of Mendeleev, 
who successfully organized them by their various unique properties. But the profusion 
of new elements coming into being led to the suspicion that there might be something 
more fundamental out of which all these chemical elements were composed of. Around 
1932, Chadwick's discovery of neutron put the Particle Physics on the solid premise for the 
first time, with electrons, protons and neutrons being established as the fundamental and 
elementary building blocks of matter. In fact, that was also the last time when Particle 
Physics was so simple and complete (apparently). As history repeated itself, all kinds 
of new particles like mesons, hadrons and other leptons were discovered, which appeared 
no less elementary than proton or neutron. That again hinted towards further elementary 
constituents of matter. Now we know that the quarks and leptons and the vector bosons that 
mediate interactions are the most elementary fundamental particles. As it stands now, in 
accordance to the Standard Model, already the number of elementary particles (considering 
different colours as different species) amounts to sixty (Higgs!) |[]][] and the number can 
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easily be increased if the nature prefers to have the symmetry like Supersymmetry. The 
story, of course, does not end here and extends to the String Theory (supposed to be the 
theory of everything) where the most elementary building block is a string and all the 
present day particles are various excitation modes of the string. Only time can tell whether 
naming quarks and gluons and so on as the elementary constituents is premature. On the 
experimental side, one of the most important high energy experiment devised to partly 
answer this question is the Deep Inelastic Scattering (DIS) experiment where a high energy 
electron probes the interior of a hadron. Thus, proper understanding of DIS phenomena is 
very important in order to find evidence regarding the basic building blocks in nature. We 
will devote a major portion of our work towards understanding the DIS phenomena. 

Now we go back and concentrate a little bit on the other aspect of the question that 
we posed, namely, what actually controls the underlying behavior when these elementary 
constituents interact among themselves. From our experience, we know that all the basic 
phenomena are governed by four kinds of underlying interactions, viz., strong, electromag- 
netic, weak and gravitational interactions. Classical principles of physics, which are so 
successful in explaining, for example, planetary motion or the motion of a billiard ball, fails 
to describe phenomena when it comes to electrons and protons. The fact that simultaneous 
measurement of two variables, if they are conjugate to each other, can not yield results with 
arbitrary level of accuracy (as first pointed out by Heisenberg), revolutionized our under- 
standing at the microscopic length scale determined by the fundamental constant h, known 
as the Planck's constant. Thus, starting from atomic to further smaller length scale, the 
description of physical phenomena is given in terms of quantum mechanical system where the 
Quantum Principles take over from Classical Principles. Also, for particles travelling with 
speed comparable to that of light, Special theory of relativity (as discovered by Einstein) 
takes over from Gallilean relativity in classical physics. Now, as it turns out, if one tries to 
keep track of particles at very small length scale and at the same time allows them to have 
very high velocities (as is the case in high energy physics), then for the sake of consistency 
what emerges is known as Quantum field theory, which actually describe the physical world. 

In our work, we shall use quantum field theory, to be more precise, light-front Hamilto- 
nian formulation of Quantum Chromodynamics and apply it to DIS, in order to understand 
the phenomena as well as the theory itself. 



II. MOTIVATION AND THE ORGANIZATION OF THE THESIS 

It is easier in explaining our motivation to introduce first the very basic concepts involved 
in the light-front QCD. Let us briefly recall some of the features of quantum field theory, in 
order to introduce the light-front QCD itself. After the success of Quantum Electrodynamics, 
it is now strongly believed that a quantum mechanical system involving particles moving 
with relativistic speed and their interactions can be best described by a local quantum field 
theory 0. The description of such a system involves infinitely many degrees of freedom. 
Various particles are described by separate field operators whose behavior under the Lorentz 
transformations determines what kind of particle these are describing, i.e, scalar, vector or 
spinor. Interactions among the particles are assumed to be local and given in terms of 
a local Lagrangian or Hamiltonian. In fact, we now know that strong, weak and electro- 
magnetic interactions can be described by similar set of gauge theories. In particular, the 
Quantum Chromodynamics (QCD) which is based on a non-abelian S'?7(3)-colour gauge 
theory describes the strong interactions among the quarks and gluons, which are carrying 
special colour charges and are the basic fields that enter into the QCD Lagrangian. In our 
work, we shall be mostly concerned with QCD. 

Now, the basic algorithm in studying such a system is to start from a suitable local 
Lagrangian and get the equation of motions satisfied by the fields using variational principle 
and then quantize the system by assuming the values of the commutators of the fields and 
their conjugate momenta on a space-like hypersurface (the so called canonical quantization). 
This algorithm goes by the name of covariant Lagrangian formulation, since it starts from 
a Lagrangian and does not need to refer, in principle, to any particular Lorentz frame. Free 
field equation of motion can be solved exactly but these solutions are not very interesting for 
they describe particles which are doing nothing but simply existing by themselves. On the 
other hand, equation of motions involving non-trivial interacting fields are in general highly 
nonlinear and rarely solved exactly. Only approximate solutions can be obtained through 
iterative methods or by using perturbative techniques. While the scattering problem can be 
treated using the whole machinery of perturbation theory and well-studied, the bound state 
problems are rarely so, let alone the analytic nonperturbative solution (save a few exactly 
solvable models). In the case of QCD low energy bound state problems, we are unable to 
carry out perturbative calculations due to the large value of the coupling in contrast to the 
situation in QED, where the coupling is small and still allows using perturbative techniques. 



The only hope for non-perturbative treatment of QCD bound state problem seems to be 
the Lattice calculations which again has its own problem. Thus, there is a need for an 
alternative to the covariant Lagrangian description where the QCD bound state problem 
can be reasonably addressed. 

On the other hand, there exists a parallel Hamiltonian formulation of field theory where 
the starting point is an Hamiltonian obtained either from the Lagrangian or constructed 
otherwise. Equations of motions are the usual Hamiltonian equations of motion as in the 
classical theory, with the Poisson brackets replaced by the commutators. Quantization 
conditions are set as in the earlier case on a space-like hypersurface which serve as the 
initial conditions and the Hamiltonian is responsible for the time evolution, i.e., gives the 
field configurations at some later time, similar to the situation in nonrelativistic quantum 
mechanics. Although the Hamiltonian formulation is somewhat natural and easy to visualize 
especially due to the background experience in the study of quantum mechanics, one has to 
pay some price for that. A hint of noncovariance creeps in while defining the Hamiltonian, 
which must be specified on a particular time slice. In principle, both the formulations should 
yield identical results if one can obtain the nonperturbative solutions and which, at present, 
is not possible. In absence of the nonperturbative solutions, one has to rely on perturbation 
techniques which not only are different for different formulations but also depend on the 
particular frame in the second case. Thus, one has to pay special attention in pursuing the 
Hamiltonian formulation so that one does not confuse the real physics with what can just 
be an artifact of a particular Lorentz frame. 

In our work, we shall be concerned with the light front Hamiltonian formulation where 
quantization surface is chosen to be the light-like hypersurface and the light front Hamil- 
tonian does the job of time evolution. To be more precise, we define first the light-front 
co-ordinates {x + , x~ ,x x } in terms of the Minkowskian co-ordinates x^: 

x ± = x°±x 3 , x L = {x l ) x 2 }, 

where x + plays the role of light-front time. Inner product is defined in a different way (see 
later), so that Lorentz invariant scalar remains the same: 

x 2 = s% = x + x~ - (x 1 ) 2 = (:r ) 2 - (X 1 ) 2 - (x 2 ) 2 - (x 3 ) 2 . 

Quantization surface is now chosen to be x + = 0, called the light-front, instead of a; = in 
usual equal-time case and light-front hamiltonian P~ (which now is conjugate to x + ) does 
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the job of new ttme-evolution. It is hoped that in this description QCD bound states can 
be addressed in a more tractable way than in other existing description. We shall explain 
why it is so expected while introducing the light-front field theory in detail and reviewing 
its special features in Chapter 2, in the context of light-front QCD. There we also give an 
overview of the selected problems some of which we attempt to understand in our work. 

Now, that we know roughly what light-front field theory is (and details will come in 
due course), it is time to spell out how we want to make use of it and what has been our 
motivation in this work. In field theory calculations, as for example in the higher order 
perturbative calculations, the resulting expressions often become infinite. These infinities 
are removed from the theory by redefining the physical parameters in the theory. This is 
known as the renormalization of the initial bare theory and is unavoidable in field theory. As 
we shall see, the light-front power counting, which is quite different from that of the usual 
covariant description, makes the renormalization procedure very different in light front field 
theory. Renormalization procedure itself becomes somewhat complicated in light-front QCD 
hamiltonian framework due to lack of covariance. Also, due to the fact that light-front QCD 
essentially is the dynamics of a constrained system, as we shall see, it is particularly im- 
portant to work in the light-front gauge A + = (see later). All these differences, namely, 
quantization procedure, renormalization and gauge-fixed noncovariant formulation, lead to 
the suspicion whether the light-front field theory description is equivalent to the usual co- 
variant one. As is obvious, it is not possible to check this by comparing intermediate steps 
of some calculation. The best way seems to be to try and apply both in a physical phe- 
nomena where experiment can settle the issue. Keeping this objective in our mind, in this 
work, we try to formulate a light-front hamiltonian description of DIS where we have ample 
experimental evidence and theoretical understanding using the usual covariant formulation 
of QCD for comparison. 

There are, of course, other motivation for building up a parallel description of DIS. First 
of all, DIS is a light cone dominated process and it is strongly believed that light-front 
description may be ideally suited here. For example, as we know, the understanding of gross 
features in DIS are most suitably realized in terms of parton distributions as provided by 
Feynman's parton model 0], which is simple and closely follows our physical intuition. In 
order to define field theoretical generalization of these parton distributions, people have used 
light-front description where it can be done most conveniently (as we shall also see). But 
barring from a few attempt M, there is no real extensive study exploring these possibilities. 



Also the fact that constituent quark model (CQM) ||, where hadrons are thought of as 
composites of minimal number of quarks and anti-quarks and is so successful in describing 
non-relativistic properties of hadrons, is not realizable in QCD due to the complicated 
structure of the vacuum state in its covariant formulation. As we shall see in Chapter 2, the 
simplicity of light-front QCD vacuum, a unique property in light-front theory, gives us the 
hope of reconciling CQM with the QCD. 

Secondly, all the state-of-the-art calculations using OPE methods or QCD improved par- 
ton model to understand DIS up till now, concentrate solely on the evolution of the DIS 
structure functions using perturbative QCD. As a consequence of using usual equal-time 
covariant field theoretic formulation in describing the structure functions, what naturally 
comes up in this analysis is the moments of structure functions (instead of the structure 
functions themselves). This point is clearly explained in Chapter 3. On the other hand, in 
DIS experiment, it is the structure functions that are directly measured. Thus, for compar- 
ing theoretical analysis with the experimental data, one has to extrapolate the measured 
structure functions (to a domain of small x-values, which is still not accessible to DIS exper- 
iments) in order to compute the moments bringing in little uncertainty. Very little emphasis 
is put actually on studying the structure functions themselves. In contrast, our approach 
addresses the structure functions themselves. 

Thirdly, as we shall explain in detail in Chapter 3, we have used the equal light-front 
current commutators and light-front version of Bj or ken- Johnson-Low (BJL) theorem to get 
explicit expression for various structure functions closely following the current algebra meth- 
ods of pre-QCD era. These expressions are obtained without recourse to the perturbation 
theory, and hence contain perturbative as well as the nonperturbative information contained 
in the structure functions. In Chapter 4, we will show introducing a new factorization 
scheme, how nonperturbative and perturbative contents of the structure functions can be 
addressed in the same light-front Hamiltonian framework. As it turned out, nonperturbative 
dynamics involved in the structure functions can be obtained from the light-front bound- 
state equation and the perturbative pieces can be easily obtained by working out the dressed 
parton structure functions. This is particularly important from the viewpoint of ongoing 
development towards solving QCD nonperturbatively using newly proposed similarity renor- 
malization scheme in the light-front QCD Hamiltonian formalism, since any breakthrough 
here can now easily be incorporated to the DIS structure functions giving rise to a com- 
plete picture. This is in contrast to the usually employed methods, which are not properly 
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equipped to address the nonperturbative information contained in the structure functions. 
Instead, the nonperturbative informations are parametrized and the main theme of these 
studies has been the (Revolution of the structure functions in terms of moments, assuming 
its form at a lower energy scale Qq. 

Lastly, since this approach is going to be quite different from the commonly used ones, we 
can expect some new results and clear understanding of certain phenomena which otherwise 
are hard to explain. As far as the light-front QCD is concerned, such an example, is shown 
in Chapter 5, which shows the implication of certain special symmetry in light-front theory 
in simplifying the higher order old-fashioned time-ordered perturbative calculations usually 
employed in light-front QCD. As far as DIS is concerned, the examples are scattered in 
Chapter 3, 4, and 6 as well as some of the references mentioned therein. Specifically, in 
Chapter 6, we shall discuss in detail, how our approach can be most naturally extended 
towards understanding the distribution of nucleonic total helicity among its constituents 
and gives us new insights to the helicity of a composite dynamical system such as the 
hadrons. 

We reemphasize the fact that our original works are presented in Chapter 4, 5 and 6; 
and also partially in Chapter 3. 
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Chapter 2 



In this Chapter, we are going to introduce the basic features of light-front field the- 
ory, in the context of light-front QCD (LFQCD), its advantages and the problems that one 
must understand for a successful practical application of the theory. Here we shall only 
be concerned with the Hamiltonian formulation of LFQCD. In the usual Hamiltonian for- 
mulation of field theory, quantization conditions in the form of commutator of dynamical 
fields and their conjugate momenta are specified on the space-like hypersurface x° = and 
the equal-time Hamiltonian is responsible for the time-evolution of the system, i.e., gives 
the field configurations at a later time. Dirac |IJ first showed that consistent field theory 
can be formulated by specifying the quantization conditions on a light-like hypersurface 
x + = x° + x 3 = (called a light front) and using a different (light front-) Hamiltonian for a 
new time (x + ) evolution. This formulation (in Dirac's terminology, "front form") is known 
as the light-front Hamiltonian field theory. On the other hand, in an attempt to see what 
happens if one traveling at the speed of light tries to formulate field theory, by boosting 
the covariant field theory results to a so called infinite momentum frame (IMF), following 
observation was made 01. Although the Lorentz transformation required to arrive at IMF is 



evidently singular (7 = 1/Jl — '^ — > ooasf^c), the singularity cancels in the calculation 
of physical objects (like Poincare generators) and results in an effective coordinate change 
given by 

x {±) =x°±x\ x L = (x 1 ,x 2 ), (2.1) 

same as the light-front co-ordinate we defined in Chapter 1. Thus, one can see the fact that 
what one obtains after going through singular limiting procedure in IMF is built in quite 
naturally in the light front field theory. Naively, using quantum mechanics analogy, one can 
expect for particles moving with relativistic speed where the particle world lines cluster near 
the light cone, the quantization on a light front and evolving the system along x + -axis by 
light-front Hamiltonian may be helpful. For this reason light-front field theory is sometimes 
also referred as the field theory of Infinite Momentum Frame. But, we should reemphasize 
that the formulation here is as prescribed by Dirac and has no connection with any singular 
limiting procedure 0. 
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1. LF dispersion relation 

Light-front co-ordinates given in eq. ( f2.1[ ) are defined in terms of x M and should not be 
thought of as a Lorentz transformation. The inner product between two four-vectors is 
defined on the light front as 

a + b~ crb + , , , . 

a.b = — + — -a ± -b ± . (2.2) 

We define the light-front four momenta as 

k {±) = k° ± k 3 , k ± = (k 1 ,k 2 ), (2.3) 

with k~ being conjugate to x + is the light front energy and k + is the longitudinal momentum. 
With the above definitions, the dispersion relation, i.e., the relation between light-front 
energy k~ and the spatial components of momenta (k + , k -1 ), for an on mass-shell particle of 
mass m, is given by, 

k 2 , +m 2 . t . 

First thing that one notices is the fact that there is no square root involved in contrast 



to the usual case E = y k 2 + m 2 . This may provide great simplification if one tries to 
solve eigenvalue equation like H | ip) = E | ip). Secondly, the numerator in eq.(^.4[) being 
always positive implies that the particles with positive light-front energy (k~) always carries 
positive longitudinal momentum (k + ). As usual, the particles with negative k~ which must 
have negative k + are mapped to antiparticles with positive k~ and k + . As a consequence, we 
always have k + > for real particles. This has great significance for the light front vacuum 
which we shall discuss shortly. Finally, k~ becomes large for the large value of k x as well 
as very small values (near zero) of k + . This makes light front renormalization aspects very 
different from the usual one. We shall see later in this Chapter how these aspects crucially 
dictate the build up of light-front theory. These drastic differences in the basic dispersion 
relation on the light-front makes light-front theory so different from the usual covariant 
theory. 

2. Simplicity of LF vacuum 

We first emphasize the uniqueness of light front vacuum as dictated by the above disper- 
sion relation. Vacuum state is always an eigenstate of the longitudinal momenta P + | 0) = 0. 
Now, the condition k + > implies that the vacuum | 0) is either devoid of particles or, at 
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most can have particles with longitudinal momenta exactly equal to zero, k + = (so called 
zero modes). Now, zero modes will have infinite energy unless both k -1 = and m — 0. This 
makes it sensible to replace the zero modes by an effective interaction, since this exactly 
is the strategy used when renormalizing the divergences away. Thus, we always consider a 
cut-off theory where longitudinal momentum is restricted to be k + > e for all the states. 
With this prescription, | 0) becomes completely devoid of any particle and therefore, an 
eigenstate of the full interacting hamiltonian with zero eigenvalue, i.e., in cut-off theory 

P+ ] 0) = => P~ | 0) = 0. (2.5) 

Thus, the light-front vacuum has become trivially simple. At the same time, it puts a 
restriction on a state with finite P + , which can now contain at most P + /e constituents. 

This is very different from the equal-time case where the vacuum has highly complicated 
structure, since it can contain, in principle, infinite number of particles moving with positive 
and negative momenta adding up to zero. This actually makes the the idea of minimal 
constituents of CQM impracticable in covariant QCD. On the other hand, this complicated 
vacuum structure is supposed to be responsible for spontaneous chiral symmetry breaking 
or confinement in QCD. It seems that with the trivial vacuum structure in light-front theory 
after removing the zero modes, we may lose these important aspects in our theory. It should 
be emphasized that we have not simply removed the zero modes from our theory. The 
longitudinal momentum cut-off (e) should be removed from the theory at the end of any 
calculation by adding necessary counter terms in the effective Hamiltonian to render the 
observables independent of e. Thus, we expect to get back all the effects of zero mode as an 
effective interaction in the Hamiltonian through renormalization. 

3. Poincare generators in LF 

Before we proceed further, let us first highlight the dynamical structure of any light- 
front theory, which again is unique and very different from the equal-time version of the 
theory. Any dynamical system can be described by ten dynamical variables (Hamiltonian 
(P°), three linear momenta (P), three angular momenta (J) and three boosts (K)). As is 
well-known, in field theory they become operators that generate the corresponding changes 
of the state vectors and known as Poincare generators, which satisfy Poincare algebra. In 
equal time theory, we know that six of them are kinematical operators {P, J} (i.e., do not 
depend on the dynamics) and the rest are dynamical {P°,K}. In light-front theory, boost 
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operators become kinematical. Longitudinal boost is like a scale transformation and the 
transverse boosts behave like Gallilean boosts in the nonrelativistic theory. On the other 
hand, two rotations (about transverse axes) which are kinematical in equal-time case become 
dynamical in light-front theory, other than the Hamiltonian itself. 

Let us consider an example which elucidates further the difference between light-front 
case and the usual equal-time one. Specifically, we consider the effect of boost on the space- 
time co-ordinates along 3-axis {K 3 ) with a speed v: 

x° = 1 (x°-(3x 3 ) } x 3 = 1 (x 3 -3x°), x 1 ' 2 = x 1 ' 2 , (2.6) 

where 3 — - and 7 = —^ — . It is clear from the above equation that under K 3 , the 
quantization surface x° = (in equal-time case) evolves to something else, which requires 
dynamical information and hence, K 3 is a dynamical generator. Introducing the parameter 
(j> such that 7 = cosh0 and d'y = sinh0, we see that, 

x+ = x° + x 3 = e~* x + , x~ = x°-x 3 = e 4, x~, x 1 ' 2 = x 1 ' 2 . (2.7) 

It clearly shows that K 3 , which is known as generator of longitudinal boost in light-front, 
behaves like a scale transformation. In particular, it keeps the quantization surface x + = 
invariant. Therefore, it is a kinematical generator in light-front theory. 

Light-front Poincare generators are obtained in the same way as in equal-time case. 
Namely, starting from Lagrangian density we obtain the energy momentum stress tensor 
T^ v and then integrate over a space-time hypersurface. Only difference being the following: 
the role of time x° in equal-time is replaced by the new time x + in the light-front and 
consequently, the integration surface is changed to a light-like surface which is normal to 
the new time x + -direction. Thus 

P" = - fdx-d 2 x x T + ", (2.8) 

M"" = - J dx~d 2 x^ [x v T + i* - x»T +v ] . (2.9) 

Note that M^ u is antisymmetric and hence has six independent components. In light-front 
dynamics P~ is the Hamiltonian and P + and P % with (i = 1,2) are the longitudinal and 
transverse momenta. M + ~ = 2K 3 and M +l = E % are the boosts. M 12 = J 3 and M~ % = F % 
are rotations. For details of the Poincare algebra in light-front see Ref. M. For our purpose, 
we notice that the boost generators form a closed algebra among themselves: 
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[E\ E 2 ] = 0, [K 3 , E l ] = iE\ (2.10) 

and 

[J 3 , EP] = ie ij E j , (2.11) 

which is similar to the generators of non-relativistic dynamics in a plane where K 3 has, 
of course, no role. This suggests the fact that there are built-in non-relativistic structures 
in the relativistic light-front dynamics, as was also evident from the dispersion relation 
in eq. (|2.4j ). In Chapter 5, we will see the consequences of such underlying structure in 
simplifying the relativistic field theory calculations using old-fashioned perturbation theory, 
which is appropriate for the light-front QCD Hamiltonian formulation. 

To sum up, we see that kinematical subgroup of the Poincare group enlarges and contains 
seven generators in light-front theory. This may prove easier in understanding a dynamical 
system, since now we can fix more variables of the system, irrespective of any knowledge 
regarding the real dynamics. Moreover, since different set of generators are kinematical in 
light-front compared to the equal-time theory, it is worth pursuing this theory, for certain 
things difficult to study in equal-time may just become simpler here. One such example is 
the feasibility of representing the QCD-bound states in terms of just a few boost invariant 
multi-particle wave-function in the Fock-space expansion, which we discuss next. 

4. Fock expansion for the bound state 

Since the Fock-states form a complete basis, any state vector, in principle, can be ex- 
panded in terms of that basis introducing corresponding amplitude for each Fock-basis. For 
example, the bound state of a hadron on light-front can be simply expanded in terms of the 
Fock states as 

\PS) = J2 dxid 2 K ±i \n,XiP + ,XiP± + K ± i,Xi)^(xi,K, ± i,Xi), (2.12) 

n,Xi 

where n represents n constituents contained in the Fock state \n, XiP + , X{P± + «_u, Aj), Aj is 
the helicity of the i-th constituent. J denotes the integral over the space: 

^Xi = l, and 5Z K -Li = ' ( 2 - 13 ) 

i i 

while X{ is the fraction of the total longitudinal momentum carried by the i-th constituent, 
and K±i is its relative transverse momentum with respect to the center of mass frame: 
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Pi± - XiP± 



(2.14) 



with p^~, Pj_i_ being the longitudinal and transverse momenta of the i-th constituent. 
$~ (xi, K±i, Aj) is the amplitude of the Fock state |n, XiP + ,XiP± + «_lj,Aj), i.e., the multi- 
parton wave function, which is boost invariant and satisfies the normalization condition: 



^2 dxid 2 K ±i \^(x h K ±i ,\i 

n.Xi 



(2.15) 



and is, in principle, determined from the light-front bound state equation, 
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(2-16) 



Here ii/^ is the interaction part of the light-front QCD Hamiltonian given later. 

In any practical application of this Fock-expansion we may face two problems. Firstly, 
as we know, each Fock-state is obtained by operating various creation operator(s) on the 
vacuum of the theory. Now, it can so happen that the vacuum already has a complicated 
structure (as is the case in equal-time theory), which may contain arbitrary number of 
particles and thereby, needs a Fock-expansion in itself. This, in effect, render the Fock- 
expansion in equal-time theory meaningless for any practical application. This is not the 
case in light-front theory due to the simplicity of the vacuum. Specifically, in the cut-off 
theory that we are going to use, it is trivial as mentioned earlier and the Fock-expansion 
can be used meaningfully. Also, the restriction k + > makes the Fock-space smaller. 

Secondly, the expansion is still infinite and it is impossible to solve the bound state 
equation, eq.( |2.16 ), which is an infinite dimensional coupled equation. To make any practical 
calculation viable using Fock-expansion, one needs to truncate the expansion at a suitable 
maximum particle number (Tamm-Dancoff truncation, TDF), with the hope that a first 
few terms in the expansion may give useful information. This truncation in light-front 
theory violates rotational invariance about the two transverse axes while in equal-time theory 
boost invariance is lost, since by restricting the particle number the possibility of creating 
or annihilating particles under these operations are restricted. It is argued 0] that the 
restoration of rotational invariance in the light-front case could be easier than the boost 
invariance in equal-time, since the rotation forms a compact group compared to the boosts 
which are noncompact. Thus, the simplicity of vacuum and the kinematical nature of boost 
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transformation enable us to make use of Fock-expansion in a practical calculation in the 
light-front QCD and we shall see such applications later in our work. 

5. Light-front QCD: Two component formalism 

Here we briefly introduce the two component formalism of light-front QCD in order to 
introduce the basic features of the hamiltonian field theory. For details see Refs. . We 
start from the QCD Lagrangian 

C = -~Tr(F^F^) + $(i^D» - m)i> , (2.17) 

where F^ = d^A u -d u A ll -ig{A^, A% A>* = J2 a ^ a T a is a (3 x 3) gluon field colour matrix, 
and T a are the generators of the SU(3) colour group: [T a , T b ] = if abc T c and Tr(T a T b ) = \5 ab . 
The field variable ip describes quarks with three colours and Nf flavours, D^ = (d^ — igA 11 ) is 
the covariant derivative, and m is an (Nf x Nf) diagonal quark mass matrix. The Lagrange 
equations of motion are: 

d^ va + gf^c A b^vc + g $YT a i; = , (2.18) 

(i7„0" - m + g^A")^ = 0. (2. 19) 

We always work in the light-front gauge A + = 0. We also define ip = ip + + ip~, where 
^ = A if) with A = ^7°7 ± . (For details of notation and convention, see Appendix A.) 
Now, in terms of these variables and in this gauge, we get the following from the above 
equations of motion, 

-(8 + ) 2 A- a = d + d { A ia + gf abc A ib d + A ic + 2^ +t T<V+ , (2.20) 

id + ip- = [a 1 - ■ (id x + gA L ) + 7 °m]^ + , (2.21) 

which are constrained equations, since these do not involve time derivative d~ . Thus, in 
light-front variables and in light-front gauge, A~ a and ip~ are constrained fields. Their 
dynamics is constrained by the dynamics of the rest of the fields A ia and tp + , which are 
dynamical fields in light-front. Therefore, the dynamics involved in LFQCD is that of a 
constrained system, which is a general feature in light-front field theory. 

To quantize such a system one may proceed according to the Dirac's method. Alter- 
natively, if one can solve the constrained fields in terms of the dynamical fields and write 
down the Hamiltonian of the system completely in terms of these dynamical fields, then the 
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canonical quantization procedure goes through considering only the dynamical fields. We 
shall follow the second path here. Conjugate momenta of the dynamical fields are given by 

**+=w^) = ¥*' < 2 ' 23) 

**" = W^r, = -¥ ' (2 ' 24) 

Now, we separate the time derivative terms in the Lagrangian, which helps us identifying 
the Hamiltonian of the system. We rewrite the QCD Lagrangian as 

C = {^F +ia (d~A ia ) + l -^ + \d~i) + ) - (d-^)i) + ) 



H- 



{A- a C a + ity-tC + Cty-}, (2.25) 



where 



H = \\{E~ a f + (B- a ) 2 ] + \{^{a x ■ {id 1 + gA 1 ) + 7 °m} + H.c.) 

+ [-d + {E~ a A- a ) - d\E ia A- a )} (2.26) 



and C = C a = are exactly identical to the constrained equations, eq.( |2.20|) and eq.( p.21| ) 



In eq.( |2.26|) we have defined E a = —\F + a and B a = F l2a as the longitudinal component 



of electric and magnetic colour fields respectively. The reason for writing the Lagrangian 
in the above form is to make the Hamiltonian density and the constraints manifest, where 
constrained fields A~ a and ip~ serves as the Lagrange's multiplier. Note that the Hamiltonian 
density 7i depends on the constrained fields A~ a and ip~ . To obtain TC in terms of the 
dynamical fields alone, one needs to solve the constrained equations to eliminate A~ a and 
ip~ in favour of A ta and ip + . In order to solve the constrained equations, we require a suitable 
definition of inverse longitudinal derivatives (^p) and (^p) . We use the following: 



^>/oo = 


1 r+oo 

- - / dy~e(x~ 
4 J-oo y V 


-y')f(y-), 


^r)V(x-) = 


= o / dy~ \x~ 

O J — oo 


-y-\ 2 f(y 



(2.27) 

which amounts to assuming antisymmetric boundary conditions for the dynamical fields 
at the longitudinal infinities, x~ — ► ±oo. See Ref. || for detail discussion on boundary 
conditions. We define two component quark fields £ as follows: 



v>~ 



r 



£r)[* i (id i + gA i )+im]Z 



(2.28) 



Thus, the LFQCD Hamiltonian can now be expressed in terms of two component dynamical 
quark and gluon fields: 



H = f dx-d 2 x ± (H + H 



int) i 



(2.29) 



where 



l-.,~.-., a) |Ct / -OT 2 + m 2 > 



Ho = -(d l A> a )(d l A> 
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H 



int 



i^-aaa \ i*-aaa ~t~ r^-aaaa T r\-aaaa T /T 



L 99S ' ' L 999 



9933 



l m 



L 9999 



(2.30) 

(2.31) 



Here interaction Hamiltonian density is split into various pieces, which give rise to various 
interaction vertices in light-front perturbation theory. Notice that 'H qqqq gives rise to four 
quark interaction which is not present in the covariant version. For complete expressions 
of various interactions see Ref . || . Here we give only those expressions which we shall use 
later on. 



n 



= g#{ - 2 (i-) (^ • A 1 ) +a ± -A ± (ip) (a 1 ■ d x + m) 
+ {^) (^ ■ d ± - m)a ± ■ A^, 
1i ggg = g f ahc [d i A 3 * A il A jc + (d 1 • A ±a ) (J-^\ A jb A jc ] . 



(2.32) 
(2.33) 



The quantization conditions are specified on the light-front as follows: 

[A ia (x), A jb (y)] x+=y+ = -iSaSij^x- - y-)5 2 (x ± - y L ) , 
{i(x),i\y)}x+= y + =5 3 (x-y), 



(2.34) 
(2.35) 



where S 3 (x — y) = 8(x~ — y~)S 2 (x — y L ). Eq.( [2.34j ) shows that the commutator between 
gluon fields themselves are nonvanishing, which is very different from the usual equal-time 
case. Also the presence of e(x~ — y~) makes it nonlocal without violating causality as will 
be discussed later on. Note the nonlocality is only in the longitudinal direction x~ and is 
one of the very important consequences of quantizing the theory on the light-front. 

In the interaction picture, the equations of motion of the dynamical fields are those of 
free fields and given by 
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d-A ia (x) = -[A ta (x),H ] 

1 roc 

= - / dy-e(x- -y-)(d ± ) 2 A ia (x + ,y-,x ± ), (2.36) 

4 J — oo 

1 r°° 
= - / dy-e(x~ -y-)[(d x ) 2 -m^+jT,^) , (2.37) 

4 J — oo 

and their solutions are 

^ (X) = ^/2^^ [ ^ a(g ' A)e "^ + i7 - C - ] ' (2 ' 38) 

£(*) = E Xa / 2 ^ + } f^ + P»(P, A)e-*» + o%, _A)e*»] , (2.39) 

with q~ = ^-f- and jT = p - + m . In eq. fl2.38|) and eq. Q2.39Q , A is defined to be 



A = < for gluons, A = < \. ,_ for quarks. (2.40) 

The gluon polarization vectors are e\ = 75(1, i) and e^ = -75(1, — i)- The quark spinors 
are simply the eigenstates of a spin-1/2 non-relativistic particle, Xi = n an d X-± = 



. The creation and annihilation operators in eq.(2.38) and eq.(2.3ii) satisfy the basic 



commutation relations 

[a(g, A), at (g\ A')] = 2(2tt) 3 q + 5 3 (q - </)<W 
{6(p,A),6t(p',A')} = Mp,A),rft(p',A')} 

= 2(27r)V<5 3 (p-p')oAA' (2.41) 

The above few paragraphs introduced the basic features of light-front field theory in the 
context of LFQCD and the two component formalism we use for LFQCD. The perturbative 
calculations in LFQCD are that of old-fashioned time ordered Hamiltonian perturbation 
theory, for a overview of which we refer Ref. ||. We will mention them as and when used. 

7. Renormalization Aspects 

In light-front field theory in the Hamiltonian framework, the renormalization is a more 
complicated issue mainly due to the noncovariant structure of the theory and quite different 
compared to the usual covariant one. This is due to the fact that the power counting in light 
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front is very different. For a detail discussion on light-front power counting, see the Ref . || . 
Here we notice the fact that only transverse directions x x carry the mass dimension, while 
the longitudinal direction x~ has no mass dimension. Thus, one has to treat transverse 
and longitudinal directions separately in determining the superficial degree of divergence 
of a divergent integral by power counting, in contrast to the covariant case where all the 
space-time directions are treated democratically. This is also evident in the single particle 
dispersion relation k~~ = - — k + m , which shows that there are two sources of divergences: 
k + — > 0+ and k —>■ oo. The divergence coming from k + —>■ 0+ is referred as infrared (IR) 
divergence, whereas k L — > oo is known as the ultraviolate divergence (UV) in light-front 
theory. 

For the above reason, dimensional regularization, which is so elegant and commonly used 
in covariant theory, is of very little importance in light-front theory. Only in the transverse 
direction, one may use dimensional regularization. But this is now no better than putting a 
simple cut-off to the transverse momenta, as shown, at least, in the context of wave-function 
renormalization in the Ref. ||. In fact, it turns out that to regularize the UV divergences 
in light-front theory cut-off regularization is the most convenient method. IR divergences 
are also regularized by putting a small longitudinal momentum cut-off, which is equivalent 
to using principal value prescription for the integration over longitudinal momenta. Also 
the fact that the light-front theory being gauge fixed and noncovariant, leads to new type 
of divergences like quadratic divergences (if we are using cut-off instead of transverse di- 
meninsional regularization) in mass renormalization or mixed divergences involving both IR 
and UV ones. To remove these divergences one has to add counter terms to the canonical 
Hamiltonian, which are often nonlocal and help restoring the invariance of the theory that 
might be broken in the process of manipulation. For detail discussion on this subject, which 
is still an unsettled issue, see the Refs. ||, 0, ||, ||. Another method specially designed to 
address the bound state problem in light-front, is that of similarity renormalization intro- 
duced by Glazek and Wilson, where first an effective Hamiltonian is obtained perturbatively, 
by giving a similarity transformation to the original Hamiltonian. A discussion on that is 
beyond the scope of this work (see for an overview, Refs. ||, ||10||). 

For our purpose, we shall use small longitudinal momentum cut-off and an UV cut-off 
for transverse momenta in light-front perturbative calculations in the context of DIS, as 
discussed later on. 
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COMMENTS 

In this Chapter, we provided an overview of some of the special features of light-front 
field theory and what can possibly be the advantages in describing a dynamical system in 



this language [11]. It all emerged from the fact that we defined our coordinates and momenta 
as well as the inner product in a different way and subsequently, quantized the theory on a 
light front while using a different "time" evolution for the system. This makes light-front 
power counting very different and hence, the renormalization. In effect, the formulation 
here becomes so different from the usual covariant field theory that it is not obvious apriori 
whether the light-front theory is equivalent to the usual one. Lot of investigations are 
currently on to establish the equivalence |T3. The best way seems to be the application of 
the theory in a physical problem where experiment can settle the issue. 

At the same time such an application, for example in DIS, gives us the opportunity to 
study QCD in a new but consistent way, which is worth pursuing in the quest for a better 
understanding, especially for the case of bound states. Feasibility of describing the hadrons 
in terms of Fock states gives the hope of reconciling constituent quark model (which has 
been so successful in explaining the hadronic spectra) with QCD. Perturbative analysis in 
this context is also interesting not only to see the predictive power of such a theory but also 
to understand the Hamiltonian renormalization itself, which is so important for building up 
a consistent light-front Hamiltonian field theory description. 

Towards this goal, we have taken up the project of investigating DIS, in an attempt 
to obtain the nonperturbative and perturbative picture involved therein, in one consistent 
language, as discussed in the rest of this thesis. 
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Chapter 3 



The basic goal in the deep inelastic scattering experiments is to probe the interior of 
a hadron target, supposedly a bound state in QCD that is readily available in the nature. 
Although the hadrons are what we have at our disposal, only quarks and gluons are the 
true dynamical degrees of freedom in QCD. Details of how these quarks and gluons conspire 
to form a low energy bound state hadron are not known and generically attributed to 
some non-perturbative confinement mechanism. On the other hand, due to asymptotic 
freedom, quarks and gluons are accessible in the very high energy in the sense that their 
dynamics can be understood in the perturbation theory. Experiments like DIS give us the 
opertunity to parameterize the non-perturbative effects in terms of some structure functions 
and study their behavior due to the quark gluon dynamics in the high energy regime as can 
be calculated using perturbative QCD. 

Historically, the observation of Bjorken scaling in the early SLAC experiments on DIS 
prompted the prediction of point like constituents in the hadron and gave birth to the quark- 
parton model as a valid framework to interpret the data in terms of parton distributions 
inside the hadron. Thus, it turned out that deep inelastic structure functions can be used 
to measure hadron's parton distributions. Later observation of logarithmic violation of 
scaling indicated that the non-abelian gauge theory of QCD might be the correct theory of 
the strong interactions. More and more accurate measurements now left very little doubt 
regarding the unique description of the structure functions in terms of perturbative QCD. 
Parton distributions may still be measured, but one must account for their evolution with 
Q 2 . The very fact that the interpretation in terms of parton distributions was successful 
in explaining the early data can now be attributed to the asymptotic freedom, one of the 
unique feature of the non-abelian gauge theories like QCD. 

Thus, at present, the goal of studying DIS is twofold - studying strong interactions in 
terms of perturbative QCD and measuring various parton distributions inside the hadrons. 
Keeping these two things in mind, we attempt to formulate a description of DIS in the 
light front hamiltonian framework which at the one end follows closely the intuitive partonic 
interpretation and at the same time, takes into account QCD with its full glory. In this 
Chapter, we first review the basic ingredients of DIS and thereby introducing the notations. 
And then we mention very briefly how one usually goes about to deal with them using 
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perturbative QCD, emphasizing the need for building up an alternative description of DIS 
what we are up to. Then we will discuss the light front current algebra and the BJL 
theorem extended to the light front framework, which happen to be our starting point 
towards building up the alternative description. Lastly, we will discuss how we obtain 
various structure functions as the Fourier transform of hadronic matrix elements of the 
bilocal currents. 

3.1. A BRIEF OVERVIEW ON DEEP INELASTIC STRUCTURE FUNCTIONS 

We begin with a brief review of the basic ingredients of lepton-nucleon deep inelastic 
scattering (DIS): 

e(k) + h(P)^e(k')+X(P + q), (3.1) 

where we have specified the four momenta of the particles explicitly and q = k — k' is the 
momentum transfer in the process through the virtual photon. The inclusive cross section 
for the above scattering process is given by 

* 1 a2E '-L^, (3.2) 



dttdE' 2M q 4 E ^ 
where E (E') is the energy of the incoming (outgoing) lepton, L^ v is the leptonic tensor, 

L ^ = 9 IZN^' sh^ik', s')u(k', s')-f u u{k, s)} 

s' 

= Ayk^ky + k v k^) Zg^ u k ■ k Zie^po-q s , (3.o) 

and W^ v is the hadronic tensor which contains all the hadronic dynamics involved in DIS 
process, 

W " V = h J ^ e* q 'HPS\mO,J u (0)]\PS) , (3.4) 

where P and S are the target four-momentum and polarization vector respectively (P 2 = 
M 2 ,S 2 = —M 2 ,S ■ P = 0), q is the virtual-photon four momentum, and J^(x) = 
J2a e aip a { x )l lJ " l l J ct{x) the electromagnetic current with quark field ip a {x) carrying the flavor 
index a and the charge e a . 

The above hadronic tensor can be decomposed into independent Lorentz invariant scalar 
functions: 
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(3.5) 



V V 

Here, in the first step, we have parametrized W^ v in terms of four scalar functions WVs with 
{i = 1, 2, 3, 4} (as is usual in the parity conserving cases), which are known as the structure 
functions. These structure functions are again functions of two independent scalar variables 
present in the problem, Q 2 = —q 2 (the negative of momentum transfered square) and the 
Bjorken scaling variable x = ^j with v = P ■ q. In the next step, it is reparametrized in 
terms of experimentally more accessible structure functions which can be written in terms 
WiS as follows. 



Fl(x,Q 2 ) 
F 2 (x,Q 2 ) 

9i(x,Q 2 ) 



- W 1 + [M 2 - 



vW 2 {x,Q 2 



{P.qf 



]W 2 



v 



W 3 (x,Q 2 ) + uW 4 (x,Q 2 



g T {x, Q 2 ) = gi (x, Q 2 ) + g 2 (x, Q 2 ) = uW 3 {x, Q 2 



(3.6) 

(3.7) 
(3.8) 
(3.9) 



Fl{x,Q 2 ) and F 2 (x,Q 2 ) contained in the symmetric part of the W^ u are known as the 
unpolarized structure functions, since only the symmetric part of hadronic tensor contributes 
in the scattering from unpolarized target. Whereas g\(x,Q 2 ) and gT(x,Q 2 ) are known as 
the longitudinal and transverse polarized structure functions respectively. The longitudinal 
and transverse polarization vector components are given by 

S ■ q 



s, 



fiL 



On, J i 



iiT , 



s, 



fJ .T 



On in 



(3.10) 



As mentioned earlier, these structure functions provide a probe to explore various aspects 
of the intrinsic structure of the hadrons. It may be worth noting that in the literature, g\ 
and g 2 are usually used to characterize the longitudinal and transverse polarized structure 
functions. However, g 2 is not really a transverse polarized structure function. It also has no 
clear physical interpretation. Only gx which can be directly measured when the target is 
polarized along the transverse direction characterizes the full information on the transverse 
polarization structure. 
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A. Scaling and scaling violation 

As noted, DIS structure functions are scalar functions of x and Q 2 . But in the DIS 
regime (i.e., Q 2 — > oo and v — ► oo with x fixed), the Q 2 dependence of the structure 
functions fades away and roughly speaking depends only on x. This phenomena as was 
first noted by Bjorken from current algebra approach and also was in reasonable agreement 
with the early SLAC data, is known as Bjorken scaling. In the parton model QXJ , where the 
hadron is supposed to be a cluster of collinearly moving, non-interacting, massless, point- 
like particles known as parton, electron-hadron scattering is viewed as the incoherent sum of 
electron-parton scattering and the scaling of the structure functions comes out automatically. 
Partons are actually quarks and gluons in reality as described by QCD and are not free. 
Thus in an asymptotically free theory like QCD scaling is expected to be violated making 
structure functions Q 2 -dependent as was also confirmed in the later accurate measurement. 
In this subsection we briefly mention how Q 2 -dependence is usually addressed using QCD. 
Our aim here is to provide the basic picture and highlight only those aspects which make 
these approaches very different from that we are going to adopt and the motivation behind 
building up such an alternate approach in the first place. 

QCD improved parton model. This approach of incorporating QCD into the DIS pic- 
ture is based on the factorization of the cross-section into 'soft' and 'hard' part which is 
deeply rooted into the parton model. The basic idea behind the intuitive parton model in 
understanding DIS is the following. If we consider the electron-hadron scattering mediated 
by a virtual photon with high energy and momentum transfer (which is the case in DIS), 
two major things happen to the hadron when looked from the centre-of-mass frame. The 
hadron gets Lorentz contracted in the direction of collision and its internal dynamics gets 
almost frozen due to time dilation. Then the virtual photon sees the hadron like a 'pan- 
cake' composed of a bunch of nearly non-interacting partons. Each of these partons may 
be thought of as carrying a definite fraction x of the parent hadron's momentum satisfying 
< x < 1, since otherwise one or more partons would have to move in the opposite direction 
to that of the hadron, an unlikely configuration. Thus the hadron could be modeled in terms 
of these parton distributions and the cross-section of the electron-hadron scattering can be 
computed approximately by summing all possible electron-parton cross-section folded with 
the probability of finding such a parton inside the hadron. So, 

eH(x,Q 2 ) = Y^ / dyfa/H(y)0- e a(x/y,Q 2 ) (3.11) 

a Jx 
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where fa/niy) gives the probability to find the parton of type a or the parton distribution 
inside the hadron and the summing over all possible y starts from x (the Bjorken variable) 
by simple kinematics. Eq.( |3.11| ) expresses the basic theme of factorization where all the 
long-distance nonperturbative effects of the dynamics known as the 'soft' part, are dumped 
into the unknown parton distributions and the 'hard' part a ea (x/y, Q 2 ) can be computed in 
the perturbation theory. 

In the QCD improved parton model calculations this 'hard' part gets QCD correction 
from the diagrams, some of which are shown in the figure below. 
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Parton Picture 



QCD corrections to order a s : 

(a)real and (b) virtual gluon emission. 



and the statement of factorization changes to the following. 

a eH {x, Q 2 ) = Y, / dyfa/H(y, n)(r ea (x/y, Q 2 , fi 2 , a s (/i 2 )) + 



(3.12) 



Notice that only the short distance QCD effects can be calculated utilizing asymptotic free- 
dom. Details of how one proceeds and calculates are unimportant for our discussion and can 
be found, for example, in Ref. f2|. It is important, however, to notice that the cross-sections 
in the figure above become singular due to gluon mass going to zero (infrared singularity) or 
due to the possibility of collinear emission of gluon (mass singularity) if quarks and gluons 
are assumed to be massless (which one generally does in such calculations). So, one needs to 
regularize them. Infrared singularities get canceled among the real and virtual contributions, 
while the mass singularity is dumped into the unknown parton distribution functions bring- 
ing in the renormalization scale \x in the picture as shown in eq.( |3.i2j ). Obviously, one needs 
to choose this scale /x to be large enough to ensure the validity of perturbative calculation 
by keeping a s (n 2 ) small. The Q 2 -dependence of the hard scattering cross-section turns out 



to be like a s lnQ 2 and since a s 



InQ 2 



one has to sum all the terms of the form (a s lnQ 



2\n 



This is generally done by what is known as leading log approximation (LLA) j3|. One thus 
eventually arrives at the Altarelli-Parisi equation, which can also be obtained otherwise us- 
ing more intuitive but less-rigorous method as was first obtained by Altarelli and Parisi [H . 
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The Q 2 -evolution of the structure functions are then studied usually by taking moments or 
using convolution method (see ref. ||). One important point in the calculation of these hard 
processes is that the cross-section (and hence the structure functions) become renormaliza- 
tion scheme dependent due to the presence of finite terms (Q 2 -independent terms) which are 
different in different renormalization scheme. This ambiguity does not bother people simply 
because all the studies so far are directed towards how the structure functions evolve with 
Q 2 (where this finite Q ^independent terms play no role) and not the structure functions 
themselves. The ellipsis in eq.( p.l2|) stands for the higher twist terms which are suppressed 
by 4a and needs separate consideration (beyond the purview of this simple factorization 
technique) for their studies. 

OPE method. The OPE method || starts from considering the virtual Compton scattering 
amplitude T^ u defined by 

T^ = i fd 4 ze iq - z (PS\T(J^z).r(0))\PS), (3.13) 

and uses the optical theorem 

W<* v = — ImT^ (3.14) 

to get predictions for the structure functions. As is well known, the products of currents 
only in the region near the light-cone z 2 ~ contribute to the DIS. So, one uses the OPE 
for the products of currents in T^ v at § 2 ^ -oo or z 2 -> 0. Pictorially it can be naively 
depicted as the following. 

O(0) 






T» v Parton Picture 

The handbag diagram here corresponds to the parton model where applying OPE (naturally 
in the free theory) gives perfect scaling with the matrix elements of the local operators O(0) 
between the target states remain unknown and are related to the parton distributions that 
has to be inferred from the experiment. Now, adding appropriate radiative gluon corrections 
to the above diagram and using the renormalization group equation, one can calculate all 
the leading InQ 2 corrections to all orders in a s . 
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To elucidate further the basic theme, we consider the light-cone behavior of the product 
of two scalar operators A and B. OPE near light-cone is given by the following expansion 

A(z)B(0)=J2cn^K 1 -.z, n Or-" n (0), (3.15) 

i,n 

where the sum is over i, the various types of local operators O(0) that may contribute and 
n denotes the spin of the operator 0(0), determined by its Lorentz transformation property. 
The string of local operators Of 1 " "^"(0) are considered to be non-singular, local, symmetric 
and traceless operators to ensure definite spin (n). The expansion parameter C™(z 2 ), the 
so called Wilson coefficients, are c-number singular functions at z 2 ~ and controls all the 
singular behavior of the product of the currents. From naive dimensional argument C"s 
may be taken to behave as 

Cnz 2 ) z ~° (Lf A+dB -^- n ™ 2 , (3.16) 

where do t denotes the naive mass dimension of the appropriate operator in eq. fl3.15j ). Thus 
the strongest singularity in the expansion is obtained for the operator with minimum twist 
t defined as 

T = d 0l -n, (3.17) 

whereas less-singular terms do not contribute, as we shall see shortly, to the leading power 
behavior in Q 2 . It should be emphasized that the Wilson coefficients are process independent 
and can be calculated using perturbative QCD, while all the specific information regarding 
the particular process is buried in the matrix elements of the local operator, implying 

(P\A(z)B(Q)\P) =J2C?(z 2 )z» 1 ....z, n (P\Or--'* n (0)\P)- (3-18) 

i,n 

Now we go back to the Compton amplitude and use the above information there. We 
suppress all the obvious Lorentz indices as well as avoid the unnecessary complication coming 
out of the target spin in the the following, in order to highlight the basic logic. 



T(x,Q 2 ) = i I d*ze iq - z {P\T{J{z)J{0))\P) 

3 J 'd*ze i «-*J2C?(z 2 )z» 1 ....z» n (P\Or"" n mP) 



i.n 



an . 

£2g m ....2g M „— - - _/ dSe^Or(^)(P|Of-^(0)|P) (3.19) 



d{iq 2 ) n 
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The last step is obtained by replacing z u s by -^- and using the relation 

d d d n 

2q fll ....2q Un + trace terms, (3.20) 



and ignoring the trace terms for their contributions will be suppressed for large Q 2 . The 
most general Lorentz structure for the matrix elements of Oi is given as 

(P|Of •~"»(0)|P) = 4 n (P w ....P Mn - M 2 g^ P„....P, n + permutations) (3.21) 

where M is the target mass and the terms proportional to g^ 1 ^ 2 are the so called trace terms. 
Afs are the numbers containing the nonperturbative information of the process. Thus, using 



eq. (|3.21|) we get, 

where we have multiplied Q 2n in the numerator and denominator. Defining the Fourier 
transform of Wilson coefficient as 

an r 

C?(Q 2 ) = (Q 2 ) n g^ J d'ze^C^z 2 ) (3.23) 

we get, 

T(x, Q 2 ) = J2 C?{Q 2 )x~ n A n t + 0{x- n+2 M 2 /Q 2 ) . (3.24) 

i.n 

The terms suppressed by -^ are either coming from target mass effect (the trace terms) or 
due to higher twist contributions not included in the simple handbag kind of diagram shown 
here. Now, these Wilson coefficients C™(Q 2 ) (which also depends on the renormalization 
scale \x and the coupling #(//)) can be shown to obey the renormalization group equation, 

(/^ + Pj~ g - lo?)cnQ 2 /fi 2 , g{»)) = , (3.25) 

which allows us to calculate the leading Q 2 -dependence in eq.( [3.24| ) to all orders in a s 
provided we know the first non-trivial (1-loop) order anomalous dimension (70™) for the 
operator o^ 1 '"^" and QCD /3-function ||. For instance, the solution of eq. (|3.25|) , which is 
of the form 

CnQ 2 /v 2 ,g(v)) = C?(l,g(Q 2 ))exp[ - / ' 7o? df , (3.26) 

1 Jo J 

can be used in eq.( |3.24|) . 
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After everything said and done for T^ v , we can use the optical theorem to pass over to 
the structure functions. Now, this passing over to the structure functions (W^) from T^ v 
is not that straightforward. Notice that the physical region for DIS is < x < 1 and clearly 
the expression that we have obtained for T(x,Q 2 ) in eq. Q3.24p diverges for this physical 
region. What is therefore needed is an analytic continuation of T(x, Q 2 ) in complex x and 
get the result corresponding to the physical region as a limit of the analytic function in the 
unphysical region. This naturally leads to the moments of structure functions, instead of 
the structure functions themselves, as is shown below. Now, T(x,Q 2 ) (eq. (|3.24|) ) is good 
enough as a function of complex x for it is analytic as |x| — > oo, and has a cut from —1 to 
+1, since this region is connected to particle production in the elastic or inelastic scattering. 
Thus, we can only isolate the coefficient of x~ n by taking the Mellin transform: 



±-. J c dxx n ~ x T[x, Q 2 ) = £ C? {Q 2 )A n i 



(3.27) 



where the contour C is chosen as in the figure below. 

Im x 
C 




Re x 



Using the analyticity of T(x, Q 2 ), the contour C can be shrunk to the cut, to obtain 

— [ / + dxx n - x T(x + ie, Q 2 ) - l + dxx n - x T(x - ie, Q 2 )] 
2m 7-i 7-i 

= ^— l + dxx 11 - 1 [T(x + ie, Q 2 )-T*(x + ie, Q 2 )] 

= / + dxx n - l 2ilmT(x + ie, Q 2 ) 

2ni 7-i 

= / + dxx n - l 2W(x,Q 2 ) 

= 4 / + dxx n ~ l W(x,Q 2 ) or, 0. (3.28) 

7o 

The result in the last line depends on the crossing symmetry (x — * —x) of Ws and the value 
of n (whether odd or even). So we get for the nonvanishing moments, 

r+l 



4 I dxx n ~ 1 W{x, Q 2 ) = J2 C^(Q 2 )A^ . 



(3.29) 
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To get back the structure functions one has to invert this equation which is done very seldom. 
Instead, one uses the solution for Cf (as given in eq. (|3.26|) ) in eq.( |3.24| ) and takes ratios of 



moments at different Q 2 , in order to get rid off Afs and compare moments at various scale 
Q 2 (see, for example, Ref. |J for details). 

Thus, we see that the methods that are usually employed in studying DIS structure 
functions mainly concentrate on the studies of moments and the Q 2 -evolution thereof, us- 
ing perturbative QCD. Structure functions themselves are either ambiguous or not properly 
addressed at all, whereas these are the objects that really go into the cross-section which 
is measured in the experiment. These methods are framed in such a way that the nonper- 
turbative informations contained in the parton distributions are separated from the very 
beginning and paid very little attention to. In QCD improved parton model, which enjoys 
the assumptions of collinearity and massless partons of the original model, this separation is 
done through factorization, while OPE does the job in the other method. Also, the intuitive 
meaning of the parton distributions gets buried in the OPE method which involves cum- 
bersome mathematics. On the other hand, a complete understanding of hadrons crucially 
depends on the knowledge of nonperturbative QCD informations contained in these parton 
distributions. So, an approach which treats nonperturbative as well as perturbative descrip- 
tions of DIS structure functions in the same framework (which is really missing) will be of 
great importance in understanding the hadrons better. This is exactly what we attempt to 
build up in our work. 

Towards building up this alternate approach, in the next section, we first review the basic 
ideas behind the current algebra (in particular, the characteristics of the light cone version 
of it) which plays a crucial role in our approach. We then show how one can get relatively 
simple expressions for the structure functions using these ideas and light-cone version of the 
Bj or ken- Johnson-Low theorem. We should emphasize that our approach closely follows the 
one proposed earlier before the advent of QCD 0, but now built within the framework of 
QCD and hence, free of most of the assumptions employed earlier as we shall see below. 

3.2. LIGHT FRONT CURRENT COMMUTATORS AND BJL THEOREM 

In the pre-QCD era, current algebra approach was proposed to study the strong interac- 
tions based on various observed symmetries of these interactions without touching upon the 
real dynamics that are involved. In particular, it was introduced to circumvent two major 
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difficulties that hindered the progress in particle physics. Firstly, the lack of proper knowl- 
edge in those precise laws which govern the processes except electromagnetism and secondly, 
an inability to solve any of the realistic models which had been proposed to explain the dy- 
namics. It was first proposed by Gellman || and then extended and used heavily by him 
and others to produce various sum rules and low energy theorems which had experimental 
consequences. The basic idea comes from the fact that electromagnetic and weak interac- 
tions of hadrons could be successfully described in terms of a current-current interaction 
Lagrangian with experimental quantities like decay width or scattering cross-sections being 
intimately related to the matrix elements of operators involving these currents. Although the 
form of these currents are unknown (or, at best, model dependent), a knowledge of current 
commutators can be exploited to obtain sum rules which heavily constrain the experimental 
quantities. As is obvious, one can have exact knowledge of these commutators in a model 
(for example, the quark model) but can only postulate in reality. 

To illustrate further, let us consider the currents J£ (x) corresponding to some approxi- 
mate internal symmetry observed in the strong interaction processes. The charges Q a (£)'s 
that generate the symmetry are defined as / d 3 xJ®{x). If the symmetry is exact, then the 
currents are conserved (i.e., d^J£(x) = 0) and the charges Q a 's are time-independent and if 
we assume that the currents themselves transform in a known fashion under the symmetry 
transformations, we have the following relations, 

[Qa,Qb] = ifabcQc (3.1) 

[Qa,J£(0)]= ifabJm (3.2) 

where f a b c are the structure constants defining the characteristic algebra of the symmetry 
group. It was postulated that even if the symmetry is not exact in reality (thereby making 
the charges time- dependent), the equal time versions of eqs. ( |ff. i| , |# ■ °A ) would still remain valid. 



One just needs to go one step further and assume the local version of eqs.( |3.2j ) in the form 
of an equal time closed algebra among the currents themselves. 

[J° a (x), Jr(0)],o =0 = ifabcJm^^) (3-3) 

Such assumptions can be generalized for spatial components of the vector currents and also 
be extended for the axial vector currents if assumed to be present. Notice the presence of 
5 3 (x) in the RHS of eg. ( PT3"D which reflects the fact that if x ^ (i.e., two points are separated 
by a space-like distance), the commutator vanishes due to causality. Details of how one gets 
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various sum rules and low energy theorems by assuming such algebra of currents are, of 
course, beyond the scope of present discussion and can be found in the Ref. 0. 

One of the most important point found in such studies in current algebra was the follow- 
ing. The assumption of the local version of the algebra in terms of currents as in eq. Q3.3| ) 
entails further complication due to the Schwinger terms which may be present there. The 
exact form of the Schwinger terms are unknown except that they are like total spatial- 
derivatives and thus compatible with eq. (|3.2| ). Thus, the algebra no longer remains closed. 
This ambiguity makes the current algebra approach less predictive than it was thought to 
be and lot of work has been done towards how to live with them. That discussion is again 
beyond the scope of the present discussion. For our purpose, it is sufficient to notice that, 
(i) by postulating equal time commutators of hadronic currents in the form of a closed alge- 
bra one could predict certain sum rules having experimental consequences, (ii) the algebra 
is compatible with causality, (iii) the presence of Schwinger terms further complicates the 
issue rendering the algebra no more closed and (iv) the exact form of the currents as well as 
the schwinger terms are not known. 

Now, one introduced another concept of going to infinite momentum frame by taking 
infinite momentum limit (p — * oo) of the matrix elements of equal time current commutators 
in deriving, for example, the fixed mass sum rules like Fubini-Dashen and Gellman sum rule. 
Although the limiting procedure was not always straight-forward and free of ambiguity, it 
was realized that the light cone behavior of the current commutators played the important 
role there. It was also shown that the DIS structure functions directly measure the matrix 
elements of the current commutators on the light cone. These observations gave a new 
direction to the current algebra studies. Thus, it was suggested and shown that assuming 
the algebra of new charges defined on the light cone (Q a (x + ) = J dx~d 2 x ± J + (x)) and hence 
the light cone current algebra, one could directly obtain the fixed mass sum rules without 
going through the cumbersome p — •» oo limit. Details of light cone current algebra and their 
consequences can be found in the ref. |7],[J]. Here we want to emphasize the most striking 
difference between equal-time and light cone algebra of currents as dictated by causality. In 
the light cone version of the current algebra that corresponds to eq.( |3.3| ), we have 

[#(aO,W)W=o= if abc J?(x)6(x-)6 2 (x ± )+ d+S^ + d^S^ , (3.4) 

where, the Schwinger terms 5^,'s are explicitly shown. As pointed out earlier, current 
commutator vanishes if the points concerned are separated by a space-like distance, i.e., 
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x 2 < 0. Thus, restricting to x° = in the equal-time case, causality forces the commutator 
to be local, i.e., nonvanishing only when x = as depicted by the 5-function in eq.( |3.3| ). On 
the other hand, restricting to x + = in the case of light cone algebra, x 2 = x + x~ — x\ = 
is maintained ii x x = 0, irrespective of a; - . Implying that the causality only enforces locality 
in x x in the light cone commutators and non-locality in x~ creeps in. We will see shortly 
how this non-locality, which is one of the unique feature of LF field theory, is reflected in the 
equal-x + commutator of currents giving rise to generalization from local to bilocal currents 
in the context of DIS, and the important role played by it. 

After the advent of QCD as the underlying theory of strong interactions, current alge- 
bra approach was mostly abandoned in favour of operator product expansion (OPE) and 
perturbative calculation of Wilson co-efficients. It was largely due to the fact that most of 
the current algebra sum rules turned out to be invalid due to the perturbative QCD correc- 
tions, except those which are protected by some conservation laws. In our way of addressing 
the problem, we closely follow the current algebra methods employed earlier taking QCD 
as the guideline. For example, now since we know that QCD is the underlying theory, we 
know the exact form of the currents in terms of the field variables and are able to calculate 
the equal-:c + commutators of currents (instead of postulating), using that among the field 
variables themselves. Thus, we attempt to supplement the current algebra approach to DIS 
by incorporating QCD, as will be discussed here and in the later Chapters. In this section, 
we first show how equal-x 4 " current commutators come into the DIS picture and then derive 
the relevant commutators. 

A. An expansion in inverse power of light-front energy of the virtual photon 

The hadronic tensor W^ v is given in terms of hadronic matrix elements of the current 
commutator as in eq. fl3.4p . Notice that it is not equal-time or equal-a; + commutator. To see 
how DIS structure functions can be related to equal-x + commutators of current, we start 
with T^, the well known forward virtual photon-hadron Compton scattering amplitude: 

T^ = i I d 4 £e tq< (PS\T(J^t)r(Q))\PS). (3.5) 

As it is already noted, the hadronic tensor is related to the forward virtual-photon hadron 
Compton scattering amplitude as 

W"" = — Mr"". (3.6) 

2n 
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Similar to the case of hadronic tensor, T^ v can be parametrized as 



■yfl „V ^ 



S a T z {x,Q 2 ) + P G S-qT^x,Qi 2 ) . (3.7) 



T- = ( - gT + q -f-)Ux, Q 2 ) + {f - - 2 q») {? ~ - 2 f)T 2 (x, Q 



- ie» vXa q x 



Using the optical theorem, we have 



,+Wi^,cr 



T i (x,Q 2 ) = 2J dq'^ ]l^ + > , i = 1,2,3,4. (3.8) 



— oo 



The above relations give us the opportunity to connect any information regarding T^ u to 
the structure functions themselves and as we shall see below, these structure functions can 
really be connected to the light-front bilocal currents through the l/q~ expansion of T^ v 
and using equal-x + current commutators. 

An expansion of T^ v in l/q~ was originally proposed by Jackiw et al. |7[] based on BJL 
theorem. The general expansion in l/q~ is given by (see Appendix B for a derivation) 

T " v = - E (-) / dc#£Le i **(PS\[{ide) n j'>{€) } r(o)y =0 \PS) , (3.9) 



n=0 



<1 



where q~ = q° — q 3 , the light-front energy of the virtual photon, and d~ = 2-jS+ is a light- 
front time derivative and (£ + , £~ , £ l ) are the light-front space-time coordinates. The above 
expansion shows that the time-ordered matrix element in T^ u can be expanded in terms of 
an infinite series of equal light-front time (i.e., equal-x + ) commutators. 

For large Q 2 and large v limits in DIS, theoretically without loss of generality we can 
always select a Lorentz frame such that the light-front energy q~ of the virtual photon 
becomes very large. Explicitly, in terms of light front variables, we can choose q + to be 
negative and finite for the virtual photon. Also, keeping q % to be finite, one can get large 
space-like q 2 (Q 2 — > oo) by taking q~ — > oo. Thus, DIS regime can simply be obtained for 
q~ — > 00, such that, 

Q 2 ~ -q + q~ -»• 00, v ► 00, x — . (3.10) 

Notice that x is positive and finite, since q + is negative while both q + and P + are finite. 
Then, in the DIS regime, only the leading term in the above expansion of T^ u given in 
eq.(|3.9|) is dominant, i.e., 

q J 
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Here we have assumed the fact that the convergence of the expansion for large q~ is not 
spoiled by the integrals of the matrix elements that occur in the subsequent terms of the 
expansion. As a result, the leading contribution to the deep inelastic structure functions 
is determined by the light-front current algebra. Of course, as mentioned earlier, we can 
compute the light-front current commutator directly and exactly from QCD (where QCD 
should be quantized on the light-front time surface £ + = £° + £ 3 = with the light-front 
gauge A+ = 0) [0,0]. Hence, all the subsequent derivations are exact within the light- 
front QCD and without further assumptions or approximations of the collinear and massless 
partons that were used in the derivations as discussed earlier. 

At this point, we should emphasize the following facts which necessitate the use of light 
front description in this context. Firstly, the above exercise with the BJL expansion can be 
performed in terms of 4j as well, thereby obtaining usual equal-time current commutator 
instead. But, now taking q° — > oo limit gives q 2 > 0, i.e., time-like q 2 which is unphysical 
for DIS. One needs to circumvent this problem by going to complex g°-plane and taking 
iq° — > oo or otherwise, bringing in complexities. Secondly, as is well known, DIS is a light- 
cone dominated process. A knowledge of the current commutator on the x + = 0-surface, 
which shares a whole line with the light-cone, is rich in information and suitable for DIS 
compared to that on x° = 0, which touches only the tip of the light-cone. Thirdly, it comes 
out, as a consequence of using light-front commutator and will be shown shortly, that the 
structure functions in the inelastic scattering is directly related to the matrix element of 
bilocal currents in a similar way as the form factor in elastic scattering is to that of local 
currents. On the other hand, since DIS needs dynamical information starting from current 
commutators defined on x° = 0-surface, one can only get an infinite set of relations between 
each moment of structure function and the corresponding term in the BJL expansion J?|. 
Lastly, before proceeding further, it should be reemphasized that we do not need to bother 
about the intricacies associated with the current algebra assumptions and their validity, 
since the current commutators that we shall be using comes directly as a consequence of 
QCD, believed to be the underlying theory. 

B. Light front current commutator 

Our next task is to evaluate the light front current commutators. Here we shall outline 
the derivation of only those current commutators which may be necessary later on. The 
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hadronic current that takes part in the electromagnetic interaction with the virtual photon 
is given by, 

r(x) = 5> a ^ a (aOWa(aO , (3-12) 

a 

where ij) a (x) is the quark field carrying the flavor index a and the charge e a . We also define 
the axial vector current for later purpose: 

TO =E e «W^iW • ( 3 - 13 ) 

a 

Various components of the vector currents explicitly expressed in terms of dynamical and 
constrained quark fields are as follows. 

J+(x) = £e a 2V>Ux)V> + a(x) (3.14) 

a 

J\%) =E e « [ip+ a (x)a i ip- a (x) + ^l a (x)aV +a (x)] (3.15) 

a 

J ~i x ) = X! 2e a ^ a (:r)^_ a (:E) • (3-16) 

In the literature, sometimes J + is denoted as the good component of the current for it does 
not contain any interaction, J 1 as the bad component for it contains interaction through 
constrained field ?/>-« and so on. Degree of badness is referred by the number of constrained 
fields present and is supposed to be connected with the twist of these operators REM. 

For our purpose, we need to calculate the commutator between J + (x) and J~(0) on the 
light front x + = 0. As usual in light front field theory, the quark field ip(x) is decomposed 
into dynamical and constrained components, 

ip{x) = ip+(x) + xp-(x) , ij)±(x) = A+ipix) , A ± = - 7 V- (3.17) 

We use the basic commutation relation on light-front for the dynamical quark fields, 

{^ + (x) , i>t(y)} x+=y+ = A + 5(x~ - y-)6 2 (x ± - y ± ) , (3.18) 

where flavour indices are implicit. The minus component ip_(x) is determined from ip+(x) 
using the constraint equation: 

iJ-(x) = -^(ia ± -D ± +(3m q )ij + (x). (3.19) 

Here we have already used the light-front gauge A+ = 0, and D± = d± — igA± is the 
transverse component of the covariant derivative, a\ = 7°7\ (3 = 7 . 
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Because of the above special property of quark (or more generally fermion) fields on 
the light-front, the light-front current explicitly depends on interaction of the theory, which 
is very different from the usual equal-time formulation. In other words, the fundamental 
interaction is manifested explicitly in the light-front current commutators. 

From Eqs. fl3.l8 ) and ( |3.19| ), we have 



A+ _ _ r -, 

{ip+(x) , ip*_(y)} x+ + = — -e{x -y ) ia±- D* ± - (3m S 2 (x± - y±) . 

Ai L ~ J 



(3.20) 



Thus, using eq.( |3.18| ) and eq.( |3.20D , after a tedious but straightforward calculation, one can 
find that 



J + (x) , J~(y) 



x~ r =y 



EW 



-e(x - y )5 2 (x ± - y ± )V a (x\y) 



+di[-e(x- -y-)5 2 (x ± -y ± )[v:(x\y) 



+ ie ij A{{x\y) 



— h.c. 



(3.21) 



where V^ and A^ are defined as the bilocal vector and axial vector currents, which are 
straightforward generalization of the corresponding local currents: 



V£(x\y) = i; a (x)^ a (y), 

A^(x\y)=^ a (x)Y^My)- 



(3.22) 
(3.23) 



As we can see, the light-front current commutators are very different from the equal- 
time current commutators. Here the commutator is indeed given by terms involving spatial 
derivatives (Schwinger terms in the current algebra language). These space-derivatives come 
from the non-locality of ip-(x) on the light-front. In the equal-time formulation, there is 
no such nonlocality involved in connection with the fermion field. Therefore one cannot 
derive such a commutator from the naive canonical equal-time commutators. Note that the 
nonlocality is only in the longitudinal direction x~ as argued earlier in the current algebra 
context, but now it comes out of direct calculation. As we will soon see in the next section it 
is these nonlocalities that lead to the simple expressions of the structure functions in terms 
of bilocal current matrix elements. This is an essential feature in the present approach that 
make the light-front current algebra specially useful in the exploration of the deep inelastic 
structure functions. Also note the fact that these bilocal currents explicitly depend on the 
interaction. But the way it is written, the form of the commutator is the same as in the 
case of free theory. Thus, the scaling of the structure functions which we know is exact in 
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the free theory, will come out unaffected unless we really grind these bilocal operators and 
try to calculate their matrix elements. 

The commutators for other current components, for example, J + and J 1 can also be 
found straightforwardly, 



J+(x) , J\y) 



x~ r =y 



E^k" 



-~e(x -y )5 2 (x ± -y ± )[V^x\y)-te^Ai(x\y) 



-- e(x~ - y~)5 2 (x ± - y ± ) [g i3 V+(x\y) + id j A+{x\ 



h.c. \ . 



(3.24) 



Thus, one can use eq.( |3.21 ) to extract the structure functions and then use eg. 03.241 ) to make 
a consistency check. 



3.3. THE GENERALIZED EXPRESSIONS FOR DEEP INELASTIC 

STRUCTURE FUNCTIONS 

Now, the Compton scattering amplitude in the large q~ limit can be immediately ex- 
pressed in terms of the hadronic matrix elements of the bilocal vector and axial vector 
currents. We consider the (H — ) component the Compton amplitude and using eq. fl3.21| ) we 
get, 

1 



T+- larg =«" --L [dCe^-/h(C)(PS\Y,el{ l -q + V-(C\0) 

1 J a *• 



- ^IK(riO) + ie ij Ai(C\0)]} - h.c. \PS) . 



(3.1) 



Notice that here we have first used partial integration which brings in q +,t in the place of 
d +)l in the commutator (eq.( |3.21| )) and neglected the surface term and then used the <5 2 (^ _L ) 
to perform the ^-integration. We introduce the form factors for the bilocal current matrix 
elements using Lorentz covariance, 

(PS\V£{£\0) - V£(Q\t)\PS) = P"V la {P\i- P) + eV 2a (P\Z ■ P) , (3.2) 

(PS\A^\0) + A£(0\O\PS) = S»A la (P 2 , £ • P) + P^ ■ SA 2a (P 2 , £ • P) 

+ eS-ZA 3 a(P\ti-P). (3.3) 



Using the definition 



we find that 



<c 



P *9+r/2 



7r J-oo q n 



(3.4) 
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/+ 



nq J-oo q'^ — q + J-oo „ 12 



+ ^q + rV 2a + ~e i] qi S 3 A la + P. 



2 . , ■ - ■ 2 - «L^" ia ' 1J 2 



5+t 



■A 



2d 



(3.5) 



where the bilocal current form factors are determined from eq. (|3.2|) and eq.(|3] 



v la = —(ps\^ a (ch + MQ)-^ a (oh + MO\PS) 



p f 



(PS\^ a (cWMo) - i> a (o)fMC)\PS) 



V 2a = ^(PS\MO{l- - ^T7 + )^(0) - h.c.\PS) 
= ±(PS\^ a (C)(l- - ^il%(0) - h.c.\PS) , 
Ma = ^r(PS\^ a (C)(l l - ^:7 + W(0) + h.c.\PS) 



Ss, 



P+ 



A 



■2a 



~ 2 -(PS\iJ a (C)(Y - ^l + )l,M0) + h.c.\PS) . 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



P+^-S't 

Notice that, since £ + '~ L = in the above expressions, it follows that the matrix elements of 
the plus and transverse components of the bilocal current yield the same form factor V\ a as 
is evident from eq.( |3.6j ) and eq.( |3.7|) (and similarly for V2 a )- 

Now, let us pick up the same (H — ) component of the hadronic tensor eq. (|3.5| ) in the 
large q~ limit: 



W + ~ = -F L + (P ± f— - 2P ± ■ q ± , 



v 
+ 2ie l] q 



2 r J2_ 
V 



q 2 



c 9\ , c 9t 

jjl r OjT — 

V V 



(3.11) 



where Sjt = Sj — S + ^f and Sjl = Sj — Sjt = 5 ,+ p+, and v = jP + q~ in the large q~ limit. 
To find the deep inelastic structure functions, we compare eq.( p.ll|) with eq.( |3.5|) through 
eq. (|3.8|) . Thus, we obtain the structure functions as given below (here in the following we 
have used the notation r\ = |P + £~). 

^^ = ^Jd V e~^^elV la (3.12) 

= ^Jdve-^^el(PS\^ a (rh + M0)-^ a (0h + MC)\PS) (3.13) 
= ^Jdve-^Y,el(PS\^ a (nilM0)-^ a (0h l ±Mr)\PS} ■ (3.14) 
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Note, the last equality is found for the first time in our work and it is discussed thoroughly 



in the work by A. Harindranath and W. Zhang [13|. We shall come back to it later. 

irP+q- 



Fl(x, Q 2 ) = --§— f d V e-^ Yl e l [( p ~ - ft) Via + rV 2a 

/( A (J J ,*. A 



x (7" " 7^7 + )^(0) - h.c.\PS) . (3.15) 

Here the first equality may be reduced to the same expression obtained by the collinear 
expansion in the Feynman diagrammatic method up to the order twist-four [I3|. But it 



is obtained directly here in the leading order in the 1/q expansion without involving the 
concept of twist expansion. The polarized structure functions come out to be, 

9i(*,Q 2 ) = ^jdve' irix j:el(A la + l -P + CM a ) (3.16) 

= ^ f dne^^£(PS$ a {rft + 7d>a(P) +?(0) 7 + 75^(DI^), (3.17) 

9t(x,Q 2 ) = ^Jdrje-^^elA^ (3.18) 

J dve-^^eKPS^iC)^ ~ ^7 + )75^(0) + h.c.\PS) . (3.19) 



8ttS 1 t 



Thus we have obtained the structure functions as the Fourier transform of the various 
matrix elements of bilocal currents, eqs. ( 3.13 - 3.19|) being the main results of this section. 



The above results are derived without recourse to perturbation theory, and also without the 
use of concept of collinear and massless partons. They are also the most general expressions 
for the leading contribution (in the l/q~ expansion, not the leading contribution in terms of 
twists) to the deep inelastic structure functions in which the target is in an arbitrary Lorentz 
frame. Some of these expressions have not ever been obtained in previous works. Note that 
the above expressions of the structure functions make sense only in the A + = gauge, 
which we have used. Otherwise, bilocal expressions should always involve a path-ordered 
exponential to ensure the gauge invariance. Since the bilocality here is in the longitudinal 
direction, only in A + = gauge the exponential factor is one and we have the simple form. 
Notice from the above expressions of the structure functions which are derived within the 
framework of light front QCD, it appears that the scaling is exact as was predicted earlier 
from two seemingly disconnected works. In the pre-QCD era Bjorken predicted scaling 
using current algebra approach and taking infinite momentum limit, which in our approach 
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is built in. On the other hand, scaling was obtained through the parton model where partons 
were treated as non-interacting particles. Now we can clearly see the connection, since the 
interaction is buried in the bilocal operators and unless we pull them apart, their form looks 
the same whether we consider the free theory or the interacting one. Also the hadronic state, 
between which these bilocal operators are sandwiched, has the substructure, resolution of 
which depends on the energy of the probe. More and more substructures are resolved as we 
increase Q 2 for the probe and the structure functions become Q 2 -dependent. In this work, 
our effort is to unmask this Q 2 -dependence which is hidden in the hadronic bound states 
IPS') by describing them in terms of multi-parton wave functions, as will be discussed in the 
next Chapter. 
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Chapter 4 



In the previous Chapter, we have seen how various DIS structure functions are related 
to the Fourier transform of matrix elements of the bilocal vector and axial vector currents 
sandwiched between the target states. As they stand, it seems that these results are similar 
to what is obtainable using light front current algebra approach where form of the current 
commutators are always assumed and the scaling appears to be exact. This is due to the 
fact that, as mentioned earlier, the interactions are buried in the bilocal operators. Also, the 
target state substructures, which depends on the energy scale of the probe, are not manifest 
in these relations. Obviously, what is important for DIS structure functions are the matrix 
elements of the bilocal operators renormalized at the physical energy scale Q of the probe. 
Since there is no reference of the underlying dynamics and the necessary renormalization 
in the current algebra approach, many of the current algebra predictions are invalidated 
in QCD perturbation theory. The only exceptions are the sum rules protected by certain 
conservation laws. In this Chapter, we shall see how this renormalization procedure can be 
carried out expanding the target state in terms of multi-parton wave-functions and using 
the old-fashioned perturbation theory appropriate to light-front QCD. First we shall see 
in Sec. 4.1 how partonic interpretation becomes apparent by introducing multi-parton wave- 
functions, where the partons are not necessarily collinearly moving or massless as postulated 
originally. Next we shall propose there a new factorization scheme which enables us (for 
certain structure functions) to separate the soft and hard parts of the dynamics. It shows 
how nonperturbative and perturbative aspects of the dynamics can be dealt with within 
the same light-front Hamiltonian framework as well as the importance of the dressed parton 
structure functions, which we calculate explicitly next in Sec. 4. 2. Then in Sec. 4. 3 we shall 
calculate the ^-structure function for a meson-like target starting from the scratch, which 
shows among other things how factorization (we introduced) is realized. Lastly, we shall 
discuss in Sec. 4. 4 the physical interpretation of the various structure functions in the light 
of sum rules they satisfy. 

4.1. UNRAVELLING THE COMPLEXITIES OF STRUCTURE FUNCTIONS 
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A. Multi-parton wave functions 

In our formulation, the structure functions are proportional to the simple hadronic matrix 
elements of the bilocal currents that are separated only in the longitudinal direction. In this 
formulation, no time evolution or propagation is explicitly involved in the matrix elements. 
Hence, unlike the OPE or the perturbative field theory descriptions of parton model, all the 
perturbative and nonperturbative dynamics here are completely carried by the structure of 
the target's bound state. To unravel this dynamics we first expand the target state in terms of 
Fock-states introducing multi-parton wave-functions. As mentioned earlier, this expansion, 
which is closer to the real physical picture probed in the experiments, is meaningful due 
to the triviality of light-front vacuum. Thus, the bound state of a hadron on light-front is 
given as 



\ PS ) = ^2 dxid 2 K, ±i \n,XiP + ,XiP ± + K± i: Aj)<&f(xi, K± i: A») 

n,Xi 



(4.1) 



where n represents n constituents contained in the Fock state \n, x t P + , XiP± + n± iy Aj), A, is 
the helicity of the i-th constituent, /' denotes the integral over the space: 



J^Xi = 1, and J2^±i = 



(4.2) 



while Xi is the fraction of the total longitudinal momentum carried by the i-th constituent, 
and K±i is its relative transverse momentum with respect to the center of mass frame: 



Pi _ p 



(4.3) 



with pf,Pi± being the longitudinal and transverse momenta of the i-th constituent. 
^(xi, n±i, \i) is the amplitude of the Fock state \n, XiP + ,XiP± + «_u,Ai), i.e., the multi- 
parton wave- function, which is boost invariant and satisfy the normalization condition: 



53 / dxid 2 K ±i \$n(xi, K ±i , \i)\ 2 = 1, 

n,\i 



(4.4) 



and is, in principle, determined from the light-front bound state equation (long-hand version 
of H\PS) = E\PS) in light-front), 



/If 2 _ V^ K i± + m i 



i=l 



. 7 ', 



qqqg 


= 







(qqq\H int \qqq) (qqq\H int \qqqg) 
{qqqg\H int \qqq) 



qqq 
qqqg 



■ (4-5) 
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Here Hi nt is the interaction part of the light-front QCD Hamiltonian given in Chapter 2. 
Of course, solving this infinite dimensional coupled equation is a huge task and we shall 
shortly see how it can, at least, be made feasible by introducing a new factorization scheme. 
Thus, the complexities of the structure functions carried by hadronic bound states are now 
translated into the language of multi-parton wave functions on the light-front, rather than 
composite operators in OPE. 

Explicitly, let us look at the structure function F 2 (x,Q 2 ) as obtained in the previous 
Chapter: 

F 2 (x,Q 2 ) 1 



x An 

1 



f dr ] e- lvx V 1 

d V e-^(PSmCh + m - ${0h + iKC)\PS) P , (4.6) 



AtiP^ 

= -^ J ' dne-^iPS^icWm -^(oh^(r)\PS) P , (4-7) 

and for the illustration purpose, consider only the relation involving the plus component 
(usually called the "good" component), 

^(f)7^(0) = 2^(r)^(0). (4.8) 

This has no explicit dynamical dependence, and has the lowest mass dimension in light-front 
(a twist-two operator in OPE language). The corresponding matrix element has straight- 
forward parton interpretation. It can be easily shown, using the Fock-expansion of the state 
IPS') as well as the dynamical fields ip + present, that on the light-front eq.( [4.6|) reduces to the 
sum of expectation values of various quark (parton) number operators, which immediately 
leads to the fact that F 2 is proportional to the parton density distributions q a (x,Q 2 ). 

F2(x,Q 2 " 



X 



£e 2 g Q (*,Q 2 ), (4.9) 



q a (x,Q 2 ) = fd 2 k ± (PS\J2bl(x,k ± ,X)b a (x,k ± ,X)\PS) 

J A 

= Jd 2 n ± J2j dx i d i K ±i \$% ta (x,Xi,K 1A ,Xi)\ 2 , (4.10) 

n,\i 

where the Q 2 -dependence is carried by the multi-parton wave functions with the active 
parton renormalized at the scale Q 2 , J' denotes the integral in the right-hand-side over the 
space of eq.(PO) except for the active parton (x, n±) = (k + /P + ,p± — xP±). (Here we have 
omitted the antiquark contributions for simplicity, see Sec. 4. 3 for the complete description). 
With this consideration it is straightforward to derive the logarithmic corrections which is 
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same as that obtained in the QCD improved parton model or in the OPE, as will be shown 
later in this Chapter. In this case, all the three descriptions are almost the same. The 
only difference here in our framework is that the perturbative QCD dynamics is transferred 
from the composite operators in OPE language to the scale-dependent multi-parton wave- 
functions on the light-front, which enables us to describe the nonperturbative dynamics in 
the same framework, as we discuss next in the context of new factorization scheme that we 
are going to introduce. 

B. A scheme for the evaluation of soft and hard contributions to deep inelastic 

structure functions 

As shown in the previous Chapter, all the derivations and discussions of the deep inelastic 
structure functions in the l/q~ expansion are rigorously carried out within light-front QCD 
and without recourse to perturbation theory. The remaining problem is how to evaluate 
various matrix elements of the bilocal currents. These matrix elements contain both hard 
and soft quark-gluon dynamics. As we have just seen, all the hard and soft dynamics 
probed through the structure functions are completely carried by the target's bound state 
in the present formulation. This is the main advantage of this formalism that allows us to 
explore the perturbative and nonperturbative contributions to the structure functions in the 
same framework. Let us first introduce in detail the new factorization scheme for such an 
exploration. 

We have already seen how the hadronic bound state is formally expressed in terms of 
Fock space expansion on the light-front (see eq.( |4.1|) ), and it is determined in principle by 
the light-front bound state equation given by eq.( [4.5|) . However, the difficulty in determining 
the wave-functions by solving eq. (|4.5|) is that the QCD Hamiltonian contains more than one 
energy scale. At different energy scales, QCD Hamiltonian can exhibit different aspects of the 
dynamics. So, let us roughly divide the quark and gluon dynamics into two energy domains, 
namely, high energy and low energy. In the high energy domain, the dynamics is controlled 
by the renormalized QCD Hamiltonian with all the constituents carrying momenta greater 
than a scale /if act (~ lGeV), which we call the factorization energy scale. This high energy 
QCD Hamiltonian describes all the hard dynamics of quarks and gluons and determines the 
hard contributions to the structure functions which can be calculated in the perturbation 
theory. In the low energy domain, the effective QCD Hamiltonian is still unknown but such 
a low energy QCD Hamiltonian should fairly determine the low energy structure of the 
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hadrons and is responsible for the soft contributions to the structure functions. 
Schematically, we may write the QCD Hamiltonian on the light-front for DIS as 

' Hq CD = J k 2 >fl 2 dkfd 2 ki±TiQ CD (ki) for hard contributions, 



ttLF 
n QCD 



Hqcd — I d'k^d 2 ki±H < Q CD (k i ) for mixed hard and soft modes , 

Hqcd — Ik 2 </x 2 dkfd 2 ki±HQ CD (ki) for soft contributions , 

(4.11) 

where Hq CD represents the canonical light-front QCD Hamiltonian (with density 'Hqcd as 
given earlier) in which the transverse momenta of all the quarks and gluons are restricted to 
be /if act < k\ < Q 2 (i.e., hard partons), and Hq CD denotes a low energy effective light-front 
Hamiltonian in which all the constituents have the transverse momentum k\ < /Xf act (soft 
partons). In addition, we also introduce a Hamiltonian Hq CD which depends only on the 
interaction part and which mixes the hard and soft partons. The source of such a mixing 
term can easily be traced out as follows. If one tries to separate all the momenta integrals 
in the light-front Hamiltonian into two parts, one ends up with a term involving canonical 
Hamiltonian density where all the momenta integrals are restricted in the high momenta 
domain and another term where all are restricted in the low momenta domain; plus there will 
be all kinds of mixing term where some of the integrals are in one domain while the others 
are not. But, we should keep in mind that the above argument is naive. It is so because 
even though this way of separating the integral gives right Hq CD , and Hq CD is one of the 
mixed term which becomes effective according to the process, Hq CD is still a different object. 
This low energy Hamiltonian is, in principle, obtained by integrating out all the modes with 
k\ > /if act from the canonical light-front QCD Hamiltonian. This also explains the presence 
of the superscript C (for canonical) in Hamiltonian density while writing down Hq CD and 
Hqcd onr y) as i n eq. (|4.11| ). Writing the light-front QCD Hamiltonian in such three parts 
will make the discussion of the perturbative and nonperturbative QCD contributions to DIS 
structure functions much more transparent, as we will see next. 
Now, the target bound state can be expressed by 

\PS) = U h \PS,fi 2 act ), (4.12) 

with 

U h = T+ exp | - I f_^ dx + (Hg CD + H% CD ) j , (4. 13) 
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tL I DC. .2 \ _ r L + m I DC. .2 



Hfc D \PS, t4J = p, \PS, 14*) ■ (4-14) 



P 



In eg. ( j4.13|) , H and # contain the interaction parts only and the mixed Hamiltonian 
Hq CD is active only in the extreme right in the time-ordered expansion of the evolution 
operator Uh (as explained later). In other words, the hard and the soft dynamics in the 
bound states are determined separately by Hq CD and Hk CD but these two contributions 
are connected by Hq CD through the action of Uh on the state \PS, /if act ). On the other 
hand, the soft dynamics, contained in \PS, A*f act ), must be solved nonperturbatively from 
eq. (|4.14j) , and the key point to solve eq.( [4.14| ) is to find the low energy effective Hamiltonian 
Hq CD . A practical procedure to find Hq CD on the light-front may be the use of similarity 
renormalization group approach plus a weak-coupling treatment developed recently |1]-§|. 
Indeed, a major effort on the study of light-front QCD is underway at present H. 

To see how the perturbative and nonperturbative QCD contributions can be separately 
evaluated in the present formalism and how these two contributions are connected by Hq CD , 
we substitute eqs.( }4.12| - f4.14D into the expressions of structure functions. Denote the struc- 



ture functions simply by Pi =: {Fl, P2, 9i,9t}, 

P(x,Q 2 ) ~ [drie^^£(PS$ a (C)r<<l>a(0)±h.c.\PS), (4.15) 

where Tj involves the Dirac 7-matrices (see the expression for structure functions derived in 
the previous Chapter). After putting complete set of states in appropriate places, it follows 
that 

F t (x,Q 2 ) = I '%-^e^ £ (PS^) act \ni)(n 2 \PS^) act ) 

a ni,n,2 

x{n x \U^ a {C)ViM^) ± h.c]u h \n 2 ) , (4.16) 

where |ni), \n 2 ) are a complete set of quark and gluon Fock states with momentum k 2 < /i£ act 
only as dictated by orthonormality of the states. This is indeed the generalized factorization 
theorem in the light-front Hamiltonian formulation. The hard contribution is described by 
the matrix element, 

(^il^pUnrV/aO) ± h.c]u h \n 2 ) , (4.17) 

which can be evaluated in the light-front time-ordered perturbation theory |J. The physical 
picture corresponds to the multi-parton forward scattering amplitude with all the internal 
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partons carrying a momentum with the transverse component k±: fif act <k\<Q 2 and the 
longitudinal momentum fraction y: x < y < 1. The soft contribution is characterized by 
the overlap of the multi-parton wave functions in different Fock states: 



(PS, n 2 fact \ni){n 2 \PS, n 2 fact ) 



(4.18) 



which contains all the quantum correlations and interference effects of multi-parton (quarks 
and gluons) dynamics in the low energy domain with k\ < ^ ct - Since we assume all the 
internal partons in the relevant DIS processes to carry high momenta /if act < k\ < Q 2 , 
the mixed Hamiltonian Hq CD in the time-ordered expansion of Uh in eq. (|4.17|) has the 
contribution only when it comes at the extreme left or extreme right in the expansion, 
where it can act on the state containing soft parton producing nonvanishing result. In 
effect, it picks up a soft parton from |m) and puts it into a high energy state. It is this 
effect that connects the hard contribution of eq. (|4.17| ) to the soft contribution in eq. (|4.18|) . 
The simple parton picture in deep inelastic processes corresponds to the case when |m) = 
\ri2) in eq.( [4.16| ) with only one parton in |m) actively participating in the high energy process, 
all others being spectators. This immediately leads to 

f 1 Q 2 

Fi(x, Q 2 ) ~ Yl e l / dyPpp'Ay, x , -^—)Qod{y, /4ct) » ( 4 - 19 ) 



act 



where the hard scattering coefficient V vv ^i is determined by 



I pp',i{y> x, 



A^fact 



d V e-^(y, k ±} s\U^ Va(Dr^ a (0) =F h.c. U h \y } k ± , s). (4.20) 



Here we have denoted \y, k±, s) (y = k + /P + ) as the active parton state. Eq. (|4.20|) means that 
we have suppressed all references to the spectators in the states |rii). The hard scattering 
coefficient is directly related to the so-called splitting function whose physical interpretation 
is the probability to find a daughter parton p' in the active parent parton p. The quantity 
q a i{y, /-4ict)> usually called the parton distribution function, is given by 



Qaiiy^Lt) = Y,\( ps ^Lt\ n )\' 



(4.21) 



where n runs over all the Fock states containing the active parton with momentum fraction y. 



Theoretically, the parton distributions are determined by solving eq. (|4.14 ). Physically, they 
contain only the quantum correlations of multi-parton dynamics but no quantum interference 
effects. Example of such distribution functions is given by eqs.( f4.9| - f4.10| ) for F 2 (x)/x (except 
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the transverse momenta are now restricted to be soft), which manifestly exhibits the simple 
parton picture. 

The above discussions indeed constitute a presentation of factorization scheme in the 
light-front Hamiltonian formulation. The leading hard contributions to the structure func- 
tions are given by the the hard scattering coefficient V pp ',i(y, x, -4r— ) and a detail calculation 

"fact 

of V pp ' t i based on the light-front time-ordered perturbative expansion of the multi-parton 
wave functions will be presented later in this Chapter. The evaluation of soft contribution 
to the structure functions given by q a i(x, /if act ), however, remains for future investigations of 
nonperturbative light-front QCD approaches to the hadronic bound states. Thus, a unified 
treatment of both perturbative and nonperturbative aspects of deep inelastic structure func- 
tions in the same framework may emerge which permits one to overcome the obstacles in 
dealing with the nonperturbative QCD dynamics in contrast to the OPE methods and field 
theoretical parton model approaches. To put this factorization scheme on a stronger basis, 
in the next section, we further investigate how this factorization scheme can be realized 
perturbatively. 

C. Factorization: A Perturbative Analysis 

In this section we show in detail, how the factorization picture emerges in a perturbative 
analysis carried over to all orders in the case where the bilocal operator involved does not 
change the particle number. The analysis leads to the concept of the structure function of 
a dressed parton in DIS. 

To explicitly demonstrate the factorization picture on the light-front, we consider the F 2 
structure function as a specific example here. For simplicity we drop reference to the flavor 
and take e n — 1, then 



F 2 (x,Q 2 



x 



1 

4vr 



-v-hP+Z- 



d£ e 



(ps I f(^ + ) t (D^ + (o) - (V> + ) f W + (r)l I ps). 



(4.22) 



From the discussion of the last section, we have 

(ni I C/,- 1 [(V + ) t (DV' + (0) - W> + )*(0)i/> + (C)]Uh I n 2 ) , 
where 



(4.23) 
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U h = T + exp { - - J dx + Pr nt {x" 

= 1 - l - f dx + Pr nt (x + ) + {- l -f f dxiPr nt (xt) H dxtP^M) + •••• (4-24) 

Z J —CO <£ J— OO J — CO 

and P~ nt = Pfrd+Pi^t 1 is denoted as the hard and mixed light-front interaction Hamiltonian. 
Let us use eq.( [4.24| ) in the expression eq.( [4.23| ) and consider the first few terms in the 



order-by-order expansion in QCD coupling constant. The lowest order term (~ g®) yields 
the low energy non-perturbative distribution function 

q ^( Xl Q") = q(x,^) = J2\(PS^Un}\ 2 

n 

= £ /Vfc^PV | bl(yP + ,k ± )b s (yP + ,k ± ) | PSn 2 ) • (4.25) 

s J 

The terms linear in the coupling constant (~ gl) is of the form 

q^(x,Q 2 )= - E(^V I n) f dxf(n \ Pr nt (4) I m)(m \ O \ p)(p \ PS^ 2 ). 

(4.26) 

Here we have put complete set of states in appropriate places and denoted 

O = ^JdCe-i p+ t- x [(iP + y(CW + (0) ~ (V> + ) f (0)V> + (r)]- (4-27) 

Since the plus component of the bilocal operator conserves particle number on the light-front, 
eq. (|4.26|) will have non- vanishing contribution only when | m) and | p) contain equal number 
of particle in the same momentum range (soft). Even if there can be possibility of having 
equal number of particles in these two states, the momentum range for individual particles 
will be different, since | m) is produced after the action of Pint{ x t) an d should contain at 
least one particle with high momentum. Therefore, by orthogonality of the states, q^'(x, Q 2 ) 
vanishes. 

Next we consider the second order contribution given by, 



q {2) (x,Q 2 )= l - ]T (PSfi 2 | n) f dxt(n | P^xf) \ m){m \0\p) 

nmpk x 

(p\ f dxiPr nt (xt)\k)(k\PS^ 2 ). (4.28) 

J — OO 

Here we have not considered the contributions where intermediate states involve vanishing 
energy denominators. These contributions are most conveniently included by introducing 
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the wave function renormalization constant associated with the parton active in the high 
energy process. Using, 

PU* + ) = e* P ^ (x+) P^(0)e-* i W* + >, (4.29) 

we perform the xf and x\ integrations producing energy denominators and the resulting 
expression is 

q^ 2 \x,Q 2 ) = £ l _ p *_ (n I PfnM | m) 

(m\0\p)(p\ p- t (0) | k)(PSfi 2 \n)(k\ PV>, (4-30) 

where the P ~ denotes the light-front energy of the state \n) and so on. Note that the states 
| n), and | k) are forced to be low energy states with (A; -1 ) 2 < /i 2 for otherwise the overlap 
(k | PSfi 2 ) will be zero. We can restrict the states | m), \ p) to be high energy states with 
{k L ) 2 > fi 2 . The bilocal operator O picks an active parton in a high energy state whose 
longitudinal momentum is forced to be xP + . Further we need to keep only terms in P~ nt 
which cause transitions involving the active parton. (Transitions involving spectators lead to 
wave function renormalization of spectator states which are cancelled by the renormalization 
process as shown explicitly later in this Chapter.) 

As we shall see shortly, to order a s (i.e., considering only q^ and q^), a straightforward 
evaluation leads to, 

g(*,Q 2 )=A/-{g(*,/i 2 ) + ^CVln^ f ^-P( x /y) q (y, M (4.31) 

where M is the wave function renormalization constant of the active parton and P is the 
splitting function. Including the contribution from the wave function renormalization con- 
stant to the same order (a s ), we get the factorized form, 

q{x, Q 2 ) = Jdy V(x, Q 2 ; y, /i 2 ) q(y, /i 2 ), (4.32) 

where the hard scattering coefficient 

a O 2 r 1 

V(x, Q 2 ; y, fi 2 ) = S(x - y) + -?-C f In ^ / dzS(zy - x)P(z) (4.33) 

271 /i z JO 

with P(x) = P(x) - 5(1 - x) Jo 1 dyP(y). 

We note that the above analysis can be carried over to all orders in perturbation theory 
with the result 
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V(x,Q 2 ;y,^ 2 )= (yP + ,k ± ,s\ U h x OU h \ yP + ,k ± ,s) 

= (yP + , k 1 , s; (dressed) | O \ yP + , k x , s; (dressed)), (4.34) 

In evaluating the above expression, only in the interaction Hamiltonians in the extreme 
left and extreme right of the time ordered product we need to keep mixture of soft and hard 
partons. This is governed by P^f ■ They are needed to cause the transition of a soft parton 
to a hard parton. In the rest of the interaction Hamiltonians occurring in the chain, the 
partons are restricted to be hard, i.e., they are determined by P^jr only. For the leading 
logarithmic evolution we are discussing here, they appear ordered in transverse momentum. 

4.2. UNPOLARIZED DRESSED PARTON STRUCTURE FUNCTIONS 

Now we turn our attention to the calculations of hard scattering coefficients, 
V(x,Q 2 ;y, fi 2 ), given by Eq. (4.34). If we set k = P, then y = 1 and the hard scatter- 



ing coefficients just become the structure functions of dressed quark and gluon targets in 
DIS, 

ff(x, Q 2 ) = ^J dve~ ivx p (ks\ [^(r)rW(O) =F h.c] \ks) p . (4.1) 

As a matter of fact, we can compute the perturbative QCD correction to the hadronic 
structure functions by calculating the structure functions of the dressed quarks and gluons. 
Here we only outline the necessary tools for calculating the perturbative contribution 
to the structure functions first (for details of two-component formalism of light-front field 
theory, see 0). In old-fashioned light-front perturbation theory, the dressed quark or gluon 
states can be expanded as follows: 

I Ps) P = U h \Ps) = VJ7{ I Pa) + £ | n) ^LgglZf) 

^ . . (m \Pr* \n)(n \ Pr™ \ Ps) 1 fA o . 

van \ r m)\ n ) 

where \Ps), the bare single particle state, and | n), the two-particle state, | m), the three- 
particle state, etc., are eigenstates of the free Hamiltonian. Introducing the multi-parton 
amplitudes (wave functions), 

, _ (n | jgg | Ps) 

s ^ (m\Pr n ?\n)(n\Prf\Ps) 

^ V {P- ~ P m ){P- ~ Pn) ' { ] 
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the expansion in Eq. ( |4.2| ) takes the form 

| Ps) p = VM{ \Ps} + Y,^n\n) + J2 ®m I m) + ...} . (4.4) 

n m 

In the above expressions, P~ nt and P^ nt are the interaction parts of the canonical light- 
front QCD Hamiltonian as given earlier, but the former contains the mixed soft and hard 
partons and the latter only has hard partons. Notice that the amplitudes in eq.( |4.^| ) are 
given in terms of the vertex function and the energy denominators which are the main 
ingredients of the old-fashioned perturbation theory. Also notice the fact that amplitudes 
for the Fock-states containing more particles involve more interaction Hamiltonian. Thus, 
in the perturbative calculations in the high energy domain, we can truncate the expansion 
in eq.( |4.4| ) reliably depending on the desired order of the the calculation (more on this 
later). With this background we now proceed to calculate the F 2 structure function for 
dressed quark and gluon targets in the perturbative region. Let us mention beforehand that 
the calculations are straightforward and we have omitted the details whenever it comes to 
putting some expression into another as will be mentioned and thereby calculating the matrix 
elements by using standard commutation relations for the creation-annihilation operators; 
other details are provided as much as is necessary. 

A. Dressed quark structure function 

The F2 structure function of a dressed quark is given by, 



4t(x,Q 2 ) 1 
1 



x 4tc 

' dr l e-^\{ks\^{C)l + ^) ~^)l + ^C)\ks) q . (4.5) 



4vrP- 

Here we are using the relation involving the plus component of the bilocal current only. 
Note, we have added an extra subscript q above (g in the equation to follow) in F 2 to 
remind ourselves that it measures the quark (gluon) distribution in the target state \q). For 
later purpose, we also introduce the gluon structure function (which measures the gluon 
distribution in the target) as defined in Ref. ||, 

F ^ >(X ' Q2) = 4^r/^ e ~^ (a ' (-) F+ " a (r)F +a u (0) + (£<- 0) I k\) g . (4.6) 

Here, the superscript \q) in F 2 (in both the equations above) denotes the the target to be a 
dressed quark. 
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In view of the discussion above, explicitly we can expand the dressed quark state in terms 
of bare states of quark, quark plus gluon, quark plus two gluons, etc, 



Pa), 



fM q W(P,a)\0) 

r — ^ / CfcrC-i CI rv-i 



01 



d /Co CL /Co , . . 



x yj2(2ir) 3 P+6 3 (P - h - A; 2 )6 t (A; 1 , aOa^fe, A 2 ) | 0) 



(4.7) 



where we have truncated the expansion at the two-particle level. The factor J\f q is the wave 
function renormalization constant for the quark and the function ip2{P, & \ ki<J\, /C2A2) is the 
probability amplitude to find a bare quark with momentum k± and helicity a\ and a bare 
gluon with momentum hi and helicity A2 in the dressed quark. 



Let us introduce the Jacobi momenta (xt, K,f- 



h + — t ■ P+ b 1 - — k 1 - -4- T - P J 
ft- — Xji , ft/- — i%- -|- dj % r 



(4.8) 



so that 



$> = 1,]T^=0. (4.9) 

i i 

The amplitude ip2 is related to the corresponding boost invariant amplitude $2 as 



/ P+lk{kt,kt) = *2(x if Kt)- 



(4.10) 



Using the notation x = X\, K\ — K and using the facts x\ + x<i = 1, K\ + k 2 = 0, we have 
(see Appendix C for details) 

$2 1 ' P2 (x,k ± ;1-x,-k ± ) ' 



M 2 - 


m 2 + (re x ) 2 
a; 


(« X ) 2 1 
1-x 



X 



£ _ T a. 1 



2(2tt) 3 V^ 



x 



Xsi 



.k — im 



1 — X 



-a — a im 



x 



Xa-[e 



f>2> ' 



(4.11) 



where M and m are the masses of dressed quark and bare quark respectively. 

Evaluating the expression in eq.( |4.5| ) explicitly, noting that in the present case the con- 
tribution from the second term (which has non-vanishing contribution if anti-quarks are 
present, see later) in this expression is zero, we get the quark structure function of the 
dressed quark 
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f tf ( '- q, w«{f(l-, 



X — X2 



+ J2 fdx 2 fd 2 ^fd 2 45(l 

x 6 2 (ni + 4) I $^' A2 (a;,4;x 2 ,4) | 2 }. (4.12) 



""l,A2 



The above equation is just a special case of eq.( |4.10| ) with the target being a dressed quark 
which is truncated at the two particle level. It makes manifest the parton interpretation of 
the quark distribution function, namely, the quark distribution function of a dressed quark 
is the incoherent sum of probability densities to find a bare parton (quark) with longitudinal 
momentum fraction x in various multi-particle Fock states of the dressed quark. Since we 
have computed the distribution function in field theory, there are also significant differences 
from the traditional parton model |7| . Most important difference is the fact that the partons 
in field theory have transverse momenta ranging from zero to infinity. Whether the structure 
function scales or not now depends on the ultraviolet behaviour of the multi-parton wave 
functions. By analyzing various interactions, one easily finds that in super renormalizable 
interactions, the transverse momentum integrals converge in the ultraviolet and the structure 
function scales, whereas in renormalizable interactions, the transverse momentum integrals 
diverge in the ultraviolet which in turn leads to scaling violations in the structure function. 
Using eq.( [4.11| ) and taking the bare and dressed quarks to be massless M = m = 0, we 
arrive at 



J2 Jd 2 ^\^(x,K' 



C1,A2 

= #!/7l±f! f d 2 K ± _J_ (4.13) 

(2vr) 3 ' 1-x J (z^) 2 V ; 

where C/ = N 2 Z l - Recalling that | $2(x, fir 1 ) \ 2 is the probability density to find a quark 
with momentum fraction x and relative transverse momentum k x in a parent quark, we 
define the probability density to find a quark with momentum fraction x inside a parent 
quark as the splitting function 

P qq (x) = C/-^. (4.14) 

The transverse momentum integral in eq. (|4.13|) is divergent at both limits of integration. 
We regulate the lower limit by \x and the upper limit by Q. Thus we arrive at, 



FJf(x,Q 2 ) r a 1+x 2 Q 2 
=Af q [5(l-x) + ^C fT -^\n- 



^(1-*) + ^— In^j. (4.15) 
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Note that \x has to be large enough so that perturbative calculation is not invalidated. The 



normalization condition (PS\P'S') 



2(27r) 3 P+(5 3 (P - P')S SS > in this case reads 



^ 



, l+x\ Q 2 
ax— In — - 



1. 



1 — x \x* 
Within the present approximation (valid only up to a s ), we can write 



(4.16) 



M„ 



2vr • 



, l + x\ Q 2 
ax In — -. 



(4.17) 



1 — x ^ 

In the second term we recognize the familiar expression of wave function correction of the 
state n in old fashioned perturbation theory, namely, J2' m Ie-e \ 2 • 



Thus, putting eq. (|4.17|) in eq.( |4.15| ), we get to order a 

■1 + x 2 



Fll\x,Q 2 ) 



2q 



X 



'P-») + S h ? C ' 



X 



5(1 



x 



dy 



1 + y 2 



(4.18) 



Note that eq. fl4.18| ) can also be written as 



-l«> 



F%(x,Q> 



X 



8(l-x) + ^C f \n® 



1+x 2 3 r/1 . 



(4.19) 



which is a more familiar expression. Note that, by construction, | $2(2^, ^ _L ) | 2 is a probability 
density. However, this function is singular as x — > 1 (gluon longitudinal momentum fraction 
approaching zero). To get a finite probability density we have to introduce a cutoff e [x g i uon > 
e), for example. In a physical cross section, this e cannot appear and here we have an explicit 
example of this cancellation. Note that the function P qq = Cf ( ^\ — h|<5(l — x) does not have 
the probabilistic interpretation since it includes contribution from virtual gluon emission. 
This is immediately transparent from the relation 



dxP qq {x) 



0. 



(4.20) 



We also note that the divergence arising from small transverse momentum (the familiar mass 
singularity) cannot be handled properly in the present calculation. This is to be contrasted 
with the calculation of the physical hadron structure function where the mass singularities 
can be properly absorbed into the non-perturbative part of the structure function. On 
the other hand, In Q 2 dependence indicates towards the logarithmic scaling violation of F2 
structure function. Of course, this is not the end of the story since in our leading order 
calculation a s is small but constant. In reality, a s ~ 1/lnQ 2 and one has to consider all 
(a s Yn.Q 2 ) n terms as in LLA. Our interest here does not permit such a discussion. 
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In view of eq.( |4.14| ), the probability density to find a gluon with momentum fraction x 
inside a parent quark is defined as the splitting function 

P Gq (x) = C, + \ X) , (4.21) 

which is obtained by replacing x with (1 — x) in P qq , since in the two particle Fock-state, rest 
of the longitudinal momentum fraction is carried by the associated gluon. This can also be 
obtained directly from the gluon distribution FJjjJ as given in eq. ( (4.6|) in the dressed quark 
state. A similar calculation as above gives us 

t 2g~ T" ill -j G/Z- 



2tt fi 2 J x 

a s . Q 2 
— In — 
2n jj, 



^\n^C f xP Gq {x). (4.22) 



and serves as a clarification. It is easy to check that 



1 

dx x 
o 



P qq {x) + P Gq {x) =0. (4.23) 



Our discussion and results for F^} follows here from eq.( |4.5| ) (or, to be more precise 
eq.( |4.6| )), which involves only the plus component of the bilocal current matrix element. 
Eq.( f4.7p suggests that one should obtain the same information regarding F<^ q from the 
matrix element of the transverse components of bilocal current as well. This is not obvious 
from the operator structure of the transverse component 

^(DtV(o) =^urhl<Mo)+? + (rhl<Mo), (4.24) 

which explicitly depends on the interaction in QCD. However, it is shown explicitly |§ with 
a similar calculation as above that although the operator structures are different for different 
components of the bilocal current, F^J comes out to be the same and have simple partonic 



picture, contrary to the popular notion [[|. Thus, the equivalence of eq.( |4.6| ) and eq.( |4.7| ), 
which we found for the first time in our work, has been established. 

B. Dressed gluon structure function 

Now we consider the calculation of structure functions where the target state is a dressed 
gluon. The dressed gluon state can be expanded as 
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P\) g = VAf g {a^P,\)\0) 



+ £ 



&f£-\ & fC-i (XK<2 CL Kn 



2(2vr) 3 P+5 3 (P-fc 1 -A; 2 



n J ^2{2-nfkt ^/2(2tt) 3 ^ 
i>2(qq)(P, A | k 1 a 1 ,k 2 (T2)b\k 1 cr 1 )d l (k2,cr2) | 0) 



+ iE 



Gt/Ci (X /Ci CLKo CI rCo 



2(2vr) 3 P+5 3 (P - fci - fc 2 ) 



2 A 1 A 2 - / ^2{2n) 3 kf ^2{2n) 3 k 
ip2(gg){P,\ | fciAi,/c2A 2 )a t (A;iAi)a t (A;2,A2) | 0)} 



(4.25) 



where two-particle sector now contains gg-states as well as two-gluon states. Note the 
symmetry factor | (in the last line) for the state having identical bosons. 
As before we introduce the boost invariant amplitudes 



A P T ^2(gg)(^ + ,^ ± ) = $2(tf)(&i,«i"), 



(4.26) 



In terms of energy denominators and relevant vertex functions, the qq wave function of the 
dressed gluon is given by 

1 



$ s 2 1 ' S2 (x,k x ;1-x,-k ± ) 



M 2 _ m2+(^)2 
x(l—x) 



X 



2(2tt) ; 



= T °xL 



a .k 



x 



a .k 



m 



-a — a 



1 — x x(l — x) 



-0~ 



X- 



S2\^p 2 l ' 



(4.27) 



where M and m are the masses of dressed gluon and the bare quark respectively. Henceforth, 
we shall work with M — m — 0. The gg wave function of the dressed gluon state is given by 

,x(l — x) 1 1 



$ 



2{gg)[X,K 



II 



-2if 



afec 



2(2vr) ; 



44(4)* 



(K^) 2 v^^T^ 



- «*5y H <J« + Sij 

x 1 — X 



(4.28) 



Now explicit calculation of gluon structure function Fi g given in eq. (|4.6| ) (but now for the 
dressed gluon target) gives the following. The contribution from the first term in eq.( 4.25 ) 
comes out to be 



7 ls)(i) 
29 



A/^l-x). 



(4.29) 



Since the the operator in F 2g contains the creation and annihilation operator (a, a)) for 
gluons only and the qq component of the state contains those of quarks, the corresponding 
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contribution to structure function is a disconnected one which we omit. The contribution 
from the gg component of the dressed gluon state comes out to be 



? |9>(3) 
23 



2ir /i 2 



X J. X , . 

H h x(l — X) 



x. 



(4.30) 



1 — x x 

We define the probability density to a find a gluon with momentum fraction x in the dressed 
gluon, P GG (x) by 



P GG {x) = 2N 



x 



+ 



1 — X 



1 — X X 

Collecting the above contributions together, we have 

,| fl ) , ..•>. . , r -. .. '■<■. , V" , , .., X 



+ x(l — x) 



<(x,Q 2 )=A^[<5(l-x) + ^ln^2iV[- 



1 — X 



+ x(l — x)]x 



(4.31) 



(4.32) 



/i^ 1 — x x 

The coefficient j\f g is determined from the longitudinal momentum sum rule for the 
dressed gluon target, namely, we require, 



dxF\ 9) (x 



dx 



/o Jo 

Thus we need to evaluate 



'FH{x)+FJ${x) 



2 (p+)2 9 



(P 1 e++(o) I P>, 



2(P+) 



2 9 



<^ 1 e + (o) 1 p> £ 



Explicit evaluation leads to 
1 



(P I e ++ (0) I P\ = —In—- / dr 
2( P +)29^l t ' g W l^/s 2 7T % 2 2 ' 



x 2 + (l-a;) 2 



X, 



(4.33) 



(4.34) 



(4.35) 



We define the probability density to find a quark with momentum fraction x in a dressed 
gluon as the splitting function P qG (x): 

1 



P qG (x) = ^[x 2 + (1 - xf 



(4.36) 



jV g \l + ^ In ^J / rfx{[x 2 + (1 - x) 2 ] + 2iV[-^- + ^^ + x(l - x)]x)l = 1. (4.37) 
2tt /// ./ L 1 — x x > J 



From eq.( |4.33|) we arrive at 

a s , Q 
— in — 

2n fi 

Thus to order a s , we have 

a s ,_ q 2 



A/; = 1 - |i In ^ / dx{[x 2 + (1 - x) 2 ] + 2N[-^— + — + x(l - x)]x\ 

— /I iX i/ J. X X 



(4.38) 



Performing the ^-integration in j\f g and putting it back into eq.( 4.32 ), we get dressed gluon 
structure function to order a,, 
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a s , Q 2 

ATI /i z 



F^x,Q 2 )=6(l-x) + ^\n- 

ATX 

{2N[[—^ + 1Z£ + a;(i - x)]x + H*(i _ x )] _ i 5 (l - x)}. (4.39) 



Including the end point (x — > 1) contributions, we define, 



x 1 — x 



11,,. ,i 1 



P GG (x) = 2N{ [j^^ + — r + x(l- x)] + -6(1 -x)}- -5(1 - x). (4.40) 

To the best of our knowledge, this is the first time gluon splitting function has been calculated 
using multi-parton wave-functions. There exist some discussions in the literature regarding 
the calculation of DIS splitting functions using the language of multi-parton wave-functions 
mainly due to Lepage and Brodsky pi. But for the gluon splitting function, they have 



simply quoted the result from Altarelli-Parisi paper [LTJj. It is easy to verify that 



\JbiAj *aj 





2P qG (x) + P GG (x) =0. (4.41) 



Note that in the last two subsections we have presented the calculations of F 2 structure 
function for a dressed parton (quark or gluon) target. These calculations can also be ex- 
tended for longitudinal structure function Fl as well as the polarized structure functions. 
In fact, some of these extensions have already been performed giving rise to interesting new 
results. For details see some of the works done by us and our collaborators |]12|| , [O]. Our 
discussion here together with these references mentioned, constitute a complete presenta- 
tion of how to perform the perturbative calculations (to the leading order) in a simple and 
straightforward way within the Hamiltonian framework of light front QCD. 

4.3. STRUCTURE FUNCTION OF HADRON: PARTON PICTURE, SCALE 
EVOLUTION AND FACTORIZATION 

In the previous sections in this Chapter, we have seen how the nonperturbative con- 
tribution to the structure functions and the scaling violations from the perturbative QCD 
corrections can be unified and treated in the same framework in our formalism. This is 
done by studying how the factorization of soft and hard part of the structure function can 
be realized, leading to the introduction of dressed parton structure functions and through 
the perturbative calculation of them. In this section, we shall address the issues associated 
with scaling violations in the structure function of the "meson-like" bound state without 
invoking the proposed factorization. This gives us the opportunity to explicitly demonstrate 
the validity of the factorization outlined before. 

64 



A. Parton picture 

Let us first discuss the emergence of parton picture for the structure function of a com- 
posite state. We expand the state | P) for qq bound state in terms of the Fock components 
qq, qqg, ... as follows. 



. __ y^ f x x f 

£rLi ..h(0^\3h+ J 



ou«2 J J2(27r) 3 k? J J2(2ir) 3 kl 



MP I fo, <ri', h, a 2 )^2((27r) 3 P+<5 3 (P - h - k 2 )tf{h, <n)S{k 2 , «r 2 ) | 0) 



+ E 



lJ/Ci & rC-t I Qjhjn d rCto I CUvq CI fvo 



M J ^2(2tt) 3 A;+ J ^2(27r) 3 kt J ^2(27r) 3 k 



MP I fci, ^i; fo, ^2; fe 3 , A 3 ) V2(2tt) 3 P+^(P - fej - A; 2 - fc 3 ) 
6^*1^1)^(^2,^)0^(^3^3)10) 

+ ... . (4.1) 

Here ip 2 is the probability amplitude to find a quark and an antiquark in the meson, ip 3 is 
the probability amplitude to find a quark, antiquark and a gluon in the meson etc. 
We introduce the boost invariant amplitudes as before, 

v / P T ^(^ + ,^)=$ 2 (x 4 ,^), 

and so on. Notice that the P + -dependences of ip 2 arid ^3 axe different. Now we evaluate 
the expression in eq. Q4.6P explicitly with the "meson-like target" given in eq. Q4.il ). The 
contribution from the first term (from the quark) in eq. Q4.6|) , in terms of $s, comes out to 
be the following. 



x 



y^ / dx 2 I d 2 K^ \ d 2 K,25(l — x — x 2 )5 2 {ki + k 2 ) I & 2 1,C72 (x, Kf; X 2 K 2 ) \ 2 

0"1,0"2 

+ ^2 / dx 2 / dx% / d 2 K^ / d 2 K,2 / d 2 K^5(l — x — x 2 — x^)5 2 {ki + k 2 + k 3 ) 

°"1i ct 2,A3 

|^ ,,,a,A8 (^«i-;^,4;a:3,4)r+- • ( 4 -3) 



Again, the partonic interpretation of the F 2 structure function is manifest in this expression. 
Using different techniques and approximations, the same result has also been obtained by 
Brodsky and Lepage |10| . 
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FS(x) 



Contributions to the structure function from the second term (from the anti-quark) in 
Eq. (|4.6|) is given by, 

y^ / dx\ \ d 2 K,\ / (i 2 K^"(5(l — x 1 — x)5 2 (ki + k 2 ) | $2 1 ' Cr2 ( x i5 K i'i x > K %) | 2 

0-1,0-2 ^ J •> 

+ ^2 dxi / dx 3 / d 2 /?^ / d 2 /^ / d 2 K 3 L 5(l — X\ — x — x 3 )5 2 {ki + k 2 + k 3 ) 

01,0-2,^3 

|$r' CT2 ' A3 (x 1 ,^;x,4;x 3 ,4)r+- • (4-4) 



Notice that eq.( |4.3|) and eq. (|4.4|) are the special case of eq.( |4.10|) with Fock expansion trun- 
cated at the three particle level. The normalization condition guarantees that 

f M ¥M + m. ]=a (4.5) 

J x x 

which reflects the fact that there are two valence particles in the meson. Since the bilo- 
cal current component J + involves only fermions explicitly, we appear to have missed the 
contributions from the gluon constituents altogether. Gluonic contribution to the struc- 
ture function F 2 is most easily calculated by studying the hadron expectation value of the 
conserved longitudinal momentum operator P + . 

From the normalization condition, it is clear that the valence distribution receives contri- 
bution from the amplitudes $2; ^3; ...at any scale fi. This has interesting phenomenological 
implications. In the model for the meson with only a quark- ant iquark pair of equal mass, 
the valence distribution function will peak at x — |. If there are more than just the two 
particles in the system, the resulting valence distribution will no longer be symmetric about 
x = | as a simple consequence of longitudinal momentum conservation and the peak shifts 
to the lower value of x. 

The eq.( [4.3|) as it stands is useful only when the bound state solution in QCD is known 
in terms of the multi-parton wave-functions. The wave-functions, as they stand, span both 
the perturbative and non-perturbative sectors of the theory. Great progress in the under- 
standing of QCD in the high energy sector is made in the past by separating the soft (non- 
perturbative) and hard (perturbative) regions of QCD via the machinery of factorization. It 
is of interest to see under what circumstances a factorization occurs in the formal result of 
eq. (fOI) and a perturbative picture of scaling violations emerges finally. We shall explicitly 



address this issue in the following section where we consider only the — as in eq (|4.3| ), which 
is sufficient to explain all the relevant issues. All the arguments can be copied for anti-quark 
contribution — . 

X 
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B. Perturbative picture of scaling violations in a bound state 

To address the issue of scaling violations in the structure function of the "meson-like" 
bound state, it is convenient to separate the momentum space into low-energy and high- 
energy sectors. Such a separation has been introduced in the past in the study of renormal- 



ization of bound state equations [H| in light-front field theory. The two sectors are formally 
defined by introducing cutoff factors in the momentum space integrals. How to cutoff the 
momentum integrals in a sensible and convenient way in light-front theory is a subject under 
active research at the present time. Complications arise because of the possibility of large 
energy divergences from both small k + and large k x regions. In the following we investigate 
only the effects of logarithmic divergences arising from large transverse momenta, ignoring 
the subtleties arising from both small x (x — > 0) and large x {x — > 1) regions and sub- 
sequently use simple transverse momentum cutoff. For complications arising from x — > 1 
region see Ref . Ql0| . 



1. Scale separation 

We define the soft region to be k 1 < /i and the hard region to be /i < k 1 < A, 
where /i serves as a factorization scale which separates soft and hard regions. Since it is an 
intermediate scale introduced artificially purely for convenience, physical structure function 
should be independent of \i. The multi-parton amplitude <3> 2 is a function of a single relative 
transverse momentum k\ and we define 

** = ($' 0<K \ <fI : (4.6) 

The amplitude $3 is a function of two relative momenta, k£ and k^ and we define 

r $| s , o<4,4</x 

I $f , < 4 < fi, ix < 4 < A . , 

3 " J $t s , fi < k± < A, < 4 < fi K } 

Let us consider the quark distribution function q(x) = defined in eq.( [4.3| ). In presence 

of the ultraviolet cutoff A, q(x) depends on A and schematically we have, 

gC*,A 2 ) = W A l* 2 | 2 +E ff I $3 I'" (4.8) 

Jo Jo Jo 

For convenience, we write, 
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q(x,A 2 )=q 2 (x,A 2 )+q 3 (x,A 2 ). (4.9) 

where the subscripts 2 and 3 denotes the two-particle and three-particle contributions re- 
spectively. Thus, schematically we have from eq.( |4.6| - pL8D , 



g(x,A 2 ) = g(:r,4) + W I $ 



h |2 
2 
A* 



+£f/>3*i 2 + i:/;/>ri 2 

r A /-A 



A* J i>- 



+£/ / 1*5* I 2 , (4-10) 



where we have defined q(x, /i 2 ) = /^ | $ 2 | 2 + Jo Jo I ^2 S | 2 - Now we investigate the 

3^ and $3" 



contributions from the amplitudes § 3 h and $3* to order a s in the following. 



2. Dressing with one gluon 

Here we consider the dressing of bare gg-state with a single gluon, i. e., wish to trun- 
cate the Fock-state expansion after three-particle state. We substitute the Fock expan- 
sion eq.( |4.1| ) in the bound state eq. Q4.5p and make projection on a three particle state 
tf(k x , ai)cft(k 2 , <7 2 )at(/c 3 , a 3 ) | 0) from the left. In terms of the amplitudes $ 2 , $3, we get, 

$^ 1CT2A3 (x, «i; x 2 , k 2 ;1-x- x 2 , « 3 ) = Mi + M 2 , (4.11) 

where the amplitudes 

Mt = h-)^—T a l V x $$ n (l - x 2 , -4; x 2 , 4) (4.12) 

E ^2(2^)3 y/l-x-x 2 

and 

M 2 = I-0=T»-=J== y 2 ^ K ±. 1 _ _^.j (4 13) 

^^2(27r) 3 V 1 -^-^ 



with the energy denominator 



_ m 2 + (4) 2 ^ 2 + (4) 2 (4^ 2 



£ = [M 2 ^^ ^^ ^ ], (4.14) 

X x 2 1 — x — x 2 

and the vertex functions are given by 

Vl = x », E ,-^i — (g± - K --'"V + ^ (gi ;'- im) ]v,.Hr (4.i5) 

1 , 1 — x — x 2 x 1 — x 2 :1 
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and 

, \ — X — %2 X 2 1 — X 21 



V 2 = X l V ' 



-CT2 



Physically eq( 4.11|) shows the relation between $2 and $3, i.e., how a three particle state 



goes to two particle state, where J\A.\ is obtained when the hard gluon attaches to the quark 
and J\A 2 is obtained when it attaches to the antiquark. 

3. Perturbative analysis 

Now we consider the various possible contributions to the wave-function $3 for different 
ranges of momenta as discussed earlier. For k± hard and k 2 soft (i. e., we concentrate on 
$3 S ), we can safely assume k^ + k^ ~ k±. Now, with this consideration, the contribution 



from M.i to $3 in eq.( f4.11| ) can be simplified and reduces to 



$3 1 1 ,<T2 ' 3 (x, Ki ; x 2 , k 2 ; 1 - x - x 2 , -k 2 ) = , T 



g x^l-x-x 2 1 



2(2vr) 3 l-x 2 (« 



^2 



— 2«i o- .k x ± . 

XaiZ r [ l-x-x 2 x J-W^ 

^(l-x 2 ,-4;x 2 ,4). (4.17) 

We see that the multiple transverse momentum integral over m and k 2 of $3 factorizes into 
two independent integrals (k 2 only occurs in the argument of $2), while the longitudinal 
momentum fraction integral over x 2 (relabelled as y in the following equation) does not and 
becomes a convolution. Thus, the contribution from .Mi to the structure function is (see 
eq.(E9) 

V Jo Ltx /i z Jx y y 

where 



E T f I *& T= ^/^ f^PU-My,^), (4.18) 

J a Jo ' 27T n z Jx y y 



On the other hand, for the same configuration (/Cf hard, k^ soft) contribution from .M2 
does not factorize and the asymptotic behaviour of the integrand critically depends on the 
asymptotic behaviour of the two-particle wave function $2- To determine this behaviour, we 
have to analyze the bound state equation which shows that for large transverse momentum 
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$ 2 (k- l ) « ( k I)2 - Thus contribution from M. 2 for scale evolution is suppressed by the bound 
state wave function. Analysis of the interference terms between .Mi and A4 2 (while calcu- 
lating | $ 3 | 2 , see eq.( |4.11| )) shows that their contribution also is suppressed by the bound 
state wave function. Same is true for ^ h . 

For the other configuration k^ soft and k 2 (i. e., now concentrating on Qf 1 , which is 
the only thing left in eq.( |4.1U| ) to be considered), contributions from .Mi and the inter- 
ference terms are suppressed by the wave function. Contribution from A4 2 factorizes both 
in transverse and longitudinal space and generate a pure wave function renormalization 
contribution: 

Notice that the soft distribution q 2 (x,/i 2 ) is a function of x and hence the ^/-integral is not 
a convolution. 

Thus, we see that even though the multi-parton contributions to the structure function 
involve both coherent and incoherent phenomena, in the hard region coherent effects are 
suppressed by the wave function and we are left with calculable incoherent contributions. 

4- Corrections from normalization condition 

In the dressed quark calculation, we have seen that the singularity that arises as x — > 1 
from real gluon emission is cancelled by the correction from the normalization of the state 
(virtual gluon emission contribution from wave function renormalization). In the meson 
bound state calculation, so far we have studied the effects of a hard real gluon emission. In 
this section we study the corrections arising from the normalization condition of the quark 
distribution in the composite bound state. 

Collecting all the terms arising from the hard gluon emission contributing to the quark 
distribution function, we have (see eq. fl4.lUp , ( [4. 18| ) and ( |4.2U| )), 



q(x,A )= q 2 (x,n ) + g 3 (x,// 



a s A 2 f l dy ,.., 

+ 7r c f ln — / — p qq(-m{y, fJ- 

27r jjl z Jx y y 



A_ 2 
7 2 

We also have a similar expression for the antiquark distribution function in F 2 . 



+ ^-C f \n^q 2 (x^ 2 )JdyP(y). (4.21) 
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The normalization condition on the quark distribution function should be such that there 
is one valence quark in the bound state at any scale Q. We choose the factorization scale 
\i = Qo- Let us first set the scale A = Q . Then we have (in the truncated Fock space) 



'o Jo 

Next set the scale A = Q. We still require 



/ dx q 2 {x, Ql) + f dx q 3 {x, Q 2 ) = 1. (4.22) 

Jo Jo 



[ dx q 2 (x,Q 2 ) + [ dx q 3 (x,Q 2 ) = 1. (4.23) 

jo Jo 

Carrying out the integration explicitly, we arrive at 

r 1 2a O 2 r r 1 

/ dx q 2 (x, Q 2 )[l + -±C S In ^ / dyP(y)] + / dx q 3 (x, Q 2 ) = 1. (4.24) 

Jo 2n Vq J Jo 

Thus we face the necessity to "renormalize" our quark distribution function q 2 (x, Ql) in 
such a way that the renormalized one now represents the true quark distribution at the scale 
Qo and contains all the dynamics up to the scale Q only and independent of any higher 
scale Q. Let us define a renormalized quark distribution function 

q?(x,Q 2 ) = q 2 (x,Q 2 )[l + 2-^C f \n^ / dyP(y)]. (4.25) 

Zn Qq Jo 

We note that the evolution of q 3 requires an extra hard gluon which is not available in the 
truncated Fock space. Thus in the present approximation q 3 (x, Q 2 ) = q 3 (x, Ql) = q 3 (x, Ql), 



so that, to order a s we get from eq. (|4.24j) 



/ dx q*(x, Ql) + f dx q 3 {x, Q 2 ) = 1, (4.26) 

Jo Jo 

which is the renormalized version of eq. (|4.22|) . We have, to order a s from eq. (|4.25 ), 



a O 2 r 1 

q2(x,Ql) = q 2 H (x,Q 2 )[l-2-^C f \n^ dyP(y)}. (4.27) 

2n Qq Jo 

Collecting all the terms, to order a s , we have the normalized quark distribution function, 

q(x,Q 2 )=q*(x,Ql) 

a O 2 r 1 r 1 

+ -iC, In ^ / dy q*(y, Q 2 ) / dz 5{zy - x) P(z) 

Z7T C^Q JO JO 

+ q 3 (x,Q 2 ), (4.28) 

noindent with P(z) = P(z) —8(z — l) / dyP(y). It is also easy to see that the normalization 
condition eq.( |4.23| ) is satisfied with the q(x,Q 2 ) given in eq.( 4.28| ), since as shown earlier 
J 1 dzP(z) = 0. 
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We see that just as in the dressed quark case, the singularity arising as x — »■ 1 from 
real gluon emission is cancelled in the quark distribution function once the normalization 
condition is properly taken into account. From this derivation we begin to recognize the 
emergence of the Altarelli-Parisi evolution equation. In fact, this is the solution of Altarelli- 
Parisi equation to the leading order in a s , with a s being constant and unrenormalized. 

Let us summarize what we have done here in a more physical terms. In this section 
we have carried out an analysis of the scale evolution of structure functions of a meson-like 
composite system. We have separated the parton transverse momenta into soft and hard 
parts. The three body wave function which is a function of two relative momenta has soft, 
hard and mixed components. The mixed components of the three body wave function which 
are functions of soft and hard momenta are responsible for the scale evolution of the soft 
part of the structure function in the perturbation theory. 

In the analysis with wave functions, there are two contributions to the three body wave 
function: One where the gluon is absorbed by the quark and second where the gluon is 
absorbed by the anti-quark (spectator). There appears a non- vanishing contribution when 
the hard gluon is absorbed by the antiquark. This corresponds to the transition caused by 
the interaction Hamiltonian when the active parton remains soft, while a hard spectator 
makes transition to a soft spectator state. This leads to wave function renormalization of 
the spectator anti-quark but this is eventually cancelled by the normalization condition as 
discussed here in detail. This justifies a posteriori the prescription given in Sec. 4.1 that we 
need to keep only those terms in P~( H ^ which cause transitions involving the active parton. 

In the wave function analysis here, we see that there are also contributions that are 
omitted a priori in the proposed scheme which lead to factorization in Sec. 4.1. All of these 
contributions are suppressed by the asymptotic behaviour of the bound state wave function 
as we have explicitly shown. In summary, the detailed analysis carried out with the help of 
multi-parton wave-functions in Sec. 4. 3. B justifies the approximations made in Sec. 4.1 which 
lead to the emergence of factorization to all orders in perturbation theory and to the simple 
scale evolution picture. 

4.4. PHYSICAL INTERPRETATION OF THE STRUCTURE FUNCTIONS 

FROM SUM RULES 

In this last section, we shall explore the physical meaning of the deep inelastic structure 
functions in our framework of light-front QCD. The physical meaning of the structure func- 
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tions can be easily understood from the sum rules they obey. These sum rules generally 
arise from the existence of conservation laws. Here we consider the sum rules for F<i in 
unpolarized deep inelastic scattering in detail which is an attempt to show that the struc- 
ture functions (in general) can be connected to the matrix elements of some physical light 
front operators and thereby, making their physical meaning obvious. We first consider the 
energy-momentum tensor in QCD, spatial integration of what gives rise to various physical 
operators which are connected to the unpolarized structure functions. 

A. Energy-momentum tensor in QCD 

The symmetric, gauge-invariant energy-momentum tensor in QCD is given by 

0"" = Ih^D u + -fD^ - F» Xa F\ a + -^(F Xca ) 2 

-g^(i-i X D x -m)i). (4.1) 

The last term vanishes using the equation of motion. Formally, we split the energy momen- 
tum tensor into a "fermionic" part 9%" and a "gauge bosonic" part 0g": 

0£" = hfi ^fD v + YD"] i>, (4.2) 

and 

0f = ~F» Xa F\ a + \g^{F X(Ta )\ (4.3) 

with F v Xa = d u A Xa — d\A v a + gfab c A v h A\ c . To be consistent with the study of deep inelastic 
structure function which is formulated in A + = gauge, we shall work in the same gauge. 

The fermionic and gauge bosonic part of the longitudinal momentum densities are given 
respectively by, 

0++ = i^ + d + ijj. (4.4) 

6+ + = -F +X F+ X = d+A l d + A l . (4.5) 

Thus the total longitudinal momentum density, 

9 ++ = i^+d+ip + d+A^+A*, (4.6) 

and is free of interactions at the operator level itself. 
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Next consider the transverse momentum density 

P 
q 2 



Of = l^h + D* + YD + U = 6+U + 9+l 2 , (4.7) 



with 



6*! = -tlrij+Dty and 6+i 2 = ^ifd+ip. (4.8) 

which depend explicitly on the interaction (see later). Now, the longitudinal and transverse 
momentum operators are obtained by spatial integration as 

P + ' i = - fdx-d 2 x ± 6 ++ '\ (4.9) 

Similarly we can take out the (H — ) component of the energy-momentum tensor whose 
spatial integration gives the Hamiltonian and can be used to find the physical meaning of Fl 
|T2|| . Here we concentrate on F 2 and work out the connection of it with 6 ++ in the following. 

B. Longitudinal momentum sum rule 

The content of the momentum sum rule is known for a long time. The sum rule simply 
says that if we add up the longitudinal momentum fractions carried by all the quarks, 
antiquarks, and the gluons in the nucleon, we should get one. Here we shall rederive it in 
our framework, paving the way for finding connections between the structure functions, in 
general, and various physical light-front operators. From the expression of F 2 in terms of 
the plus component of the bilocal current matrix element given in Eq.( [4.6|) , we have, 

Fi = ^Jd V e-^Y, e l( P I fe(D7 + i(0)-t(0)7 + i(D] I P) (4.10) 

Now for the first term in the integral, we use translational invariance, 

t(f)T + ^(o) = $ a (Ph + M-C), (4-ii) 

and change the variable £~ — > — £~, to obtain, 

F ^^j d ^ nx - ^E4(^ i vaoh + van i p) 

= 4^W/^ [e ^ +e "**]IX< P I ft.(°)7 + lMr) I P). (4-12) 

Notice that we have replaced x by -M-, which changes the sign of the second exponential. 
Now, we integrate over all possible values of x, so that, 
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_ dxF 2 (x)= ^T J ^WrtU-oo ^^ + ^^ ? el{P ' ^a(0)7 + ^«(D I P> 

= 4^/ rf ^( 2(27r)5(r?) )? e ' (P ' ^°^ + ^~) I P ) 

= 4jjE^ I ^(0) T + 9+Va(0) I P). (4.13) 

Last step is obtained first by partial integration and ignoring the surface term assuming the 
fields to vanish at £,~ — > ±00; and then doing the ^-integration by using 5(rj). Note that we 
have used ■§- = -p+d + . Now, using the crossing symmetry 

F 2 {x) = F 2 {-x) 

in LHS and the fact that physical structure function is nonvanishing only for x lying between 
to 1, we finally obtain the sum rule 

jf 1 dxF 2 (x) = ^y 2 £ e 2 a (P I 0g I P) . (4.14) 

The above sum rule shows that F 2 is connected to the fermionic part of the light-front 
longitudinal momentum density Og . 

Similarly, we can define the "gluon structure function" || as mentioned earlier, 

F 2g (x) = -J^Jdrie-**{P I (-)F +w (DF + :(0) + (C - 0) | P) , (4.15) 

where the second term obtained by interchanging the longitudinal coordinates (£~ <-> 0), 
ensures the crossing symmetry F 2g (x) = F 2g (—x). Now a similar exercise as in the case of 
F 2 gives, 

o ' dxF 2g (x) = i^y 2 {P I *++ I P) ■ (4-16) 

This shows the connection between F 2g and the gluonic part of the 6 ++ , similar to the quark 
case. Now, if we assume e a = 1, the physical picture of the F 2 structure function functions 
become even more clear, since the sum of eq.( |4.14j) and eq. (|4.16|) then is nothing but the 
statement of the longitudinal momentum sum rule. 

Jl dx[F 2 + F 2g }= ^y 2 [(p 1 eg 1 p) + (p 1 eg 1 p) 

= ^ky {Ple++lP) = 1 > (417) 

since from eq.( |4.9| ) it follows (P \ 9 ++ \ P) = 2P + P + , where we have used the normalization 
of state (P I P') = 2(27r) 3 P+(5 3 (P - P'). 
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Similarly, starting from Eq. (|4.7|) in terms of the transverse component of bilocal current 
matrix element and following exactly the same procedure used earlier, we obtain, 

1 



J^pf^eKPSmPS) . 



(4.18) 



Notice that 9p a does not appear to be the same as 9+£ as defined in eg. (P~7|) . However, 
explicit demonstration (see Chapter 6) shows that 



tl)^ l iD + ip = i/j^ l id + i(j = i/j~/ + iD l ijj 



(4.19) 



up to certain surface terms which do not contribute to the transverse momenta operator 
P % as given in eq. (|4.9|) . Thus, if we are interested only in P l , 6p l a is indeed equivalent to 
9q£ an d the sum rule in eq.( |4.18| ) suggests that the integral of F 2 over all x measures the 
expectation value of the fermionic part of transverse momentum density in the target state. 
From eq.( [4.14|) and eq. (|4.18|) we see that 



> 



1 



3+i 

q 



P), 



(4.20) 



P) = —(P\ 8n 

which shows the underlying covariance of the theory at the level of matrix elements. The 
demonstration here shows that (as also is suggested from the explicit calculation in Ref. 
||) interaction dependence of the transverse bilocal current matrix element is spurious and 
both plus and transverse components contain the same information and have same partonic 
interpretation. It clearly shows that drawing conclusions by looking at the operator structure 
is quite misleading, as is the case with transverse component of bilocal current operators 
that are twist three in the working definition of twist introduced by Jaffe || . 

Since F 2 involves quark charges in specific combinations, it does not give the direct test 
of the above momentum sum rule. To test the sum rule experimentally, one can combine the 
data for both the electron-proton and electron-neutron deep-inelastic scattering and assume 
that the sea is flavor symmetric, then 

5 1 



dx 



*?(*) + Fl 



x 



9(P+) 2 
5 1 
9P+P] 



Y,(ps\&£\ps) 

a 

Y,(ps\<\ps)- 



(4.21) 



This shows that I / dx 



F?\ 



x) + F" 



IX 



is the total longitudinal momentum fraction carried 
by the all the quarks in proton and neutron. If the quarks carry all the momentum, then 
we expect that 
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dx 



F?(x) + FS 



x 



5 
9 



(4.22) 



Experimental data shows that the above integral is 0.28. In other words, as is well-known, 
there are gluons and sea-quarks, and half of the momentum in hadrons are carried by them. 
Similar physical interpretation for the structure functions F^ and gi are already worked 
out by showing their connections with the matrix elements of appropriate physical light- 
front operators (for details we refer to the original works 1 12] , [[R|). Thus, in the light-front 
Hamiltonian formulation, physical interpretation of the DIS structure functions are most 
conveniently obtained as they can be shown to be connected with some physical light-front 
operators via the sum rule. 
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Chapter 5 



After obtaining satisfactory results in the leading order calculation in bare coupling 
as presented in the previous Chapter, it is natural to extend the dressed quark structure 
function to the next higher order, namely, to the fourth order 0(g A ). However, a complete 
fourth-order calculation is always very involved one ||. In light-front, it is particularly 
so for the renormalization becoming more complicated due to lack of covariance. Instead, 
we recall that in the dressed parton structure function calculations to the lowest non-trivial 
order, we truncated the Fock-expansion of the dressed parton state at the two-particle level 
and obtained 






Also recall that the coupling constant a s is a bare parameter in our theory. For the purpose 
of fourth-order extension, we need to include three-particle sector. For example, in the 
dressed quark state, one includes states where a bare quark is associated with two gluons or 
a pair of quark and anti-quark. To this order, we expect to obtain (from our experience in 
LLA) 

FJM An x ■ <Q 2 ) „ , Q'p ( , a 1 f a 'Y^ Q 2 Y f'^r p ( \n p ( x \ 
—^ = ^--) + ^rC f \n-P qq ( X ) + -^-j (in-j I jC f P qq (y)C f P qq (-). 

Notice that, in the leading order term, we have the renormalized running coupling in QCD 
(which, strictly speaking, contains fourth order term in bare coupling), while the last term 
corresponds to the well known ladder contribution involving two rungs where the coupling is 
still a bare parameter and needs further higher order calculation for its running. Thus, the 
coupling constant renormalization being a part this calculation, one should have a complete 
understanding of it in light-front theory as a first step towards working out the fourth order 
extension of the structure function and we shall do that in this Chapter. 

Now, the coupling constant renormalization as required for the fourth order extension of 
structure function calculation is also very important from a different viewpoint as we discuss 
now and we shall put major emphasis on that in our discussion. In light-front field theory, 
at present, many higher order calculations need to be performed using time (x + ) ordered 
perturbative techniques in order to overcome several conceptual and practical problems [|3],|J . 
Investigations have revealed a very complex structure for the coefficient functions (accompa- 
nying the divergences) which emerge at the end of notoriously long and tedious calculations 
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of individual time (x + ) ordered diagrams which are higher order in the coupling. It appears 
that almost no guidance is available to look for possible mistakes in these structures. In 
contrast, in covariant field theory, the structures accompanying the divergences are quite 
simple. The complexity of the former is due to the fact that power counting is different on 
the light-front ||. In the latter case, simplicity of the structure is due to the underlying 
Lorentz symmetry (rotational and boost invariance) which can be maintained at every stage 
of the calculation. Since the light-front formalism do possess some kinematical symmetries, 
it is worthwhile to investigate whether they can provide some constraint on the possible 
structure of coefficient functions for individual x + ordered diagrams. 

Two of the most important kinematic symmetries in light-front field theory which are 
relevant especially for phenomenological concerns are the longitudinal and transverse boost 
symmetries. As we discussed in Chapter 2, the longitudinal boost symmetry is a scale sym- 
metry on the light-front whereas transverse boost symmetry is simply Gallilean symmetry in 
two dimensions in non-relativistic dynamics || . The implications of the Gallilean symmetry 
for the structure of the interaction vertices resulting from the light-front Hamiltonian at tree 
level are known. For example, the symmetry can be utilized || to reduce the number of free 
parameters in an interaction Hamiltonian constructed at tree level purely from light-front 
power counting. The implications of this symmetry beyond tree level is not well-understood. 

Incidentally, we mention that previous calculation J7| of vertex corrections have employed 
the four-component representation of Brodsky and Lepage || and only the final answers after 
summing different time orderings have been presented which is of no use to us. A calculation 
@ (more suitable to our purpose) of the vertex correction using the two-component repre- 
sentation |T0|JTT| have studied only the two specific cases of helicity-flip part of the vertex 
(proportional to quark mass) and zero momentum (q +,± = 0) limit for the gluon. But for 
the structure function calculation, it is necessary that the particles connected to the vertex 
have generic momenta, which makes the calculation complicated exactly in a way as men- 
tioned earlier. Therefore, the calculation of vertex correction with general kinematics gives 
us the opportunity to investigate whether and how the Gallilean symmetry manifests itself 
beyond tree level and whether the symmetry can provide some guidance in understanding 
the complex structure of coefficient functions accompanying the divergences. As we shall see 
here, explicit calculation reveals that the accompanying structures are either proportional 
to the canonical vertex or independent of the total transverse momentum and thereby ex- 
hibiting transverse boost invariance. Since, these processes are relevant for the calculation 
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of asymptotic freedom in light-front QCD we also present the /^-function calculation for the 
most general kinematics in the two-component formalism. 

We reemphasize the fact that our motivation here is to study the role played by Gallilean 
boost symmetry in ensuring the correctness of the structure of coefficient functions appearing 
in the calculation beyond tree level for each x + ordered diagram separately. Such nontrivial 
checks are going to be extremely helpful in any higher order calculation using old-fashioned 
time-ordered perturbation theory used in LFQCD, as is the case in extending the calculations 
of structure functions presented earlier to the fourth order in coupling. 

5.1. TRANSVERSE BOOST SYMMETRY: CANONICAL CONSIDERATIONS 

In light-front theory, the generators of transverse boost are given by 



E l = M +i = l - I (Lv-(l 2 .v J 



x* e ++ - x + e +i 



(5-1) 



where 9^ u is the symmetric energy- momentum tensor. As discussed in Chapter 2, the gen- 
erators E l leave x + = invariant and hence are kinematic operators. We recall that the 
motion in the transverse plane generated by E l and J 3 are Gallilean in nature. This is 
ensured by the fact that the commutation relations satisfied by them, namely, 



E\E 3 



, J 3 , E i = ie ij E j , (5.2) 



where e iJ the two dimensional antisymmetric tensor, ressemble the corresponding Gallilean 
generators in the non-relativistic dynamics. Also the commutation relation 

" E\ P j ] = - i5 ij P + , (5.3) 

shows that the light-front longitudinal momentum plays the role of non-relativistic mass, 
which occurs instead of P + in the corresponding commutation relation in the non-relativistic 
case. Thus the generators E l s act just like Gallilean boosts in the transverse plane, familiar 
from non-relativistic dynamics. 

In light-front theory involving fermions and gauge bosons, the interaction vertices have a 
nontrivial structure (see Vi below). Gallilean symmetry implies that the interaction vertices 
in the theory (in momentum space) are independent of the total transverse momentum in 
the problem. Since the issues associated with Gallilean invariance are most transparent in 
the two-component representation (which we have always been using in this dissertation), 



it is most convenient to use this representation of light-front QCD in contrast to the more 

familiar four- component representation ||. 

The canonical quark-gluon vertex in our notation is (see Fig. 1) 
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Fig.l The canonical quark-gluon vertex in light-front QCD. 
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(5.4) 



Since the mass term (helicity-flip interaction) is irrelevant for the Gallilean invariance, we 
drop it in the following. Note that the canonical vertices and energy denominators in 
x + ordered diagrams in the two-component representation are manifestly invariant under 
Gallilean boost. (See Appendix D for an explicit example at the one loop level). 

5.2. ONE LOOP CALCULATIONS 



In the massless limit, the helicity-flip contribution vanishes and the canonical vertex has 
the structure 
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(5.1) 



q^- p x pi 

In this section we consider corrections to this vertex at one loop level in LF Hamiltonian 
perturbation theory. Specifically we consider the corrections arising from quark-gluon vertex 
and the three-gluon vertex. Note that the corrections arising from instantaneous vertices in 
the theory do not contribute to the divergent structure of the vertex for zero quark mass at 
one loop level and hence will not be considered here at all. 

In order to perform the calculations beyond tree level, we need to regulate the loop mo- 
menta. How to introduce regulators in light-front theory is, at present, an active subject of 



research fL2 |. One may (1) introduce cutoffs on the sum of light-front energies (the so-called 
boost invariant cut-off), or (2) choose to cutoff energy differences at vertices (which emerge 
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naturally in similarity renormalization perturbation theory) or (3) simply cutoff single par- 
ticle momenta. We employ the third choice for the regulators, namely, we put cutoffs, such 
that kf > e, \x < kf- < A, which is simple to implement but obviously violate both lon- 
gitudinal and transverse boost invariance. Since the vertices and energy denominators are 
explicitly invariant under the Gallilean boost (see Appendix D), the violation of the symme- 
try can occur only through the explicit appearance of total transverse momentum P x in the 
limits of integration. From power counting, the vertex corrections at one loop level are only 
logarithmically divergent in the transverse plane. Hence we expect the coefficient functions 
accompanying the logarithms to still exhibit the symmetry. 

Consider the one loop corrections to the vertex involving two quark-gluon vertices. There 
are two time-ordering contributions shown in Figs. 2(a) and 2(b). 
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Fig. 2 Contribution to the quark-gluon vertex from contributions involving two quark-gluon vertex. 



The contribution from Fig. 2(a) is 
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(5.2) 



Pi fci 

Here k^ = q(+^ + k^'^ and ^ +,±) = ^ +1±) - k^. 

The calculation is cumbersome but straightforward. Here we mention the important 
points that one should keep in mind while actually performing the calculations. First of 
all we notice that, there are two sources of divergences: one coming from the lower limit 
of integration over longitudinal momenta k + and the other coming from the high value 
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of transverse momenta k^-. The /c + -integration for various terms yields either finite or 
logarithmically divergent contribution. The leading divergence coming from the transverse 
momenta integration is only logarithmic. But to obtain the correct structure of the co- 
efficient functions accompanying the log-divergence, one has to be careful. We briefly discuss 
the general procedure that can be pursued in evaluating the expression in eq. ( p.2| ) as well as 
for the other diagrams in Fig.2 and Fig.3. Notice that all the vertex functions involved in 
eq. (|5.2|) depend on the integration variable k and are linear functions of k x . Similarly, both 
the energy denominators are quadratic functions of k x . Thus, the general structure of the 
transverse momenta integration becomes 

(A± ■ k 1 )^ 1 ) 2 + B 1 (k 1 ) 2 + (Ai ■ k 1 ) + B 2 



dk + d z k A 



(5.3) 



[(^) 2 + (ai-k ± ) + 61] [(A^) 2 + (4 • k 1 ) + b 2 ] ' 

where the coefficients Aj~s and BiS as well as a^s and 6jS are independent of k x and com- 
pletely determined by the external particles' quantum numbers (momenta, helicity etc.) 
and the longitudinal momenta k + . To determine large k x behavior we expand the energy 
denominators as 
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[(A;^) 2 + (of- • fc-L) + &i] (k^y 
Putting this back into the general expression and dropping the terms which vanish either 
due to large values of k 1 - or being an odd function of A;- 1 , we get 
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where we have used 



{Ai-k L ){al-k L ) = C l {k L )' 



(5.5) 



(5.6) 



Now k + dependence of (Bi —C\ —C 2 ) is such that the integration over k + yields either finite 
or logarithmically divergent contribution, i.e., 



fdk + (B 1 -C 1 -C 2 ) = A + B\n— + C\n— + D\n — 
Je e e e 



(5.7) 



where pf, p 2 and q + are the longitudinal momenta of the external particles (see the Fig.2(a)) 
and e is a small cut-off used for the lower limit of fc + -integration. Here A, B etc. are known 
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once B\, C\ and Ci are determined through explicit calculation. For example, we can 
easily read out C and A from eg. (|5~8|) and eq.( |5.9|) respectively where we presented the 
contributions from the diagram in Fig. 2 (a); B and D are zero for this diagram. Thus, in 
general, we expect two types of divergence structures: one involving single logarithm such 
as ~ In - and the other involving a product of logarithms such as ~ In — In %—. 

Following the above procedure, after a long and tedious calculation, we arrive at two 
types of divergent contributions, as mentioned, viz, one containing product of logarithms 
and other containing a single logarithm. We have taken jjl and A to be much larger than the 
external momentum scales in the problem. Divergent contributions that contain products 
of logarithms: 
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Since the coefficient of the divergent factor is proportional to the canonical vertex, the 
transverse boost invariance of the above result is manifest. Divergent contributions that 
contain single logarithm are 
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In this case the coefficient of the divergent factor involving transverse momenta is not pro- 
portional to the canonical vertex. However, in terms of the internal momenta (see Appendix 
D), the quantity inside the square bracket can be rewritten as 
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which satisfies the constraint from Gallilean invariance, namely, independent of the total 
transverse momenta P . 

Contribution from Fig. 2(b) is 
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Here fci" 1- '" 1 "'' = pi^'^ ; — /4 an d ^3 = g ( - + '~ L ' ) — /4 ' ■ The overall negative sign arises 
from the anti symmetry property of fermionic states. Note that this negative sign is missing 
fromEq. (A8) of Ref. @. 

As in the previous case, explicit evaluation leads to terms containing two types of diver- 
gences. Divergent contributions that contain products of logarithms are 
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Again, the transverse boost invariance of this result is manifest since the contribution is 
proportional to the canonical vertex. Divergent contributions that contain single logarithm 
are 
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The transverse boost symmetry of the terms inside the square bracket is not manifest but 
becomes explicit once we express the result in terms of the internal momenta. Alternatively, 
by subtracting and adding the term — 6^+ to these terms we can rewrite the terms inside 
the square bracket as the canonical term plus the terms contained in the square bracket in 
eq. ( |5.9| ) which again shows the boost invariance of the result in eq. ( |5.13| ) . 
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Fig. 3 Contribution to the quark-gluon vertex from contributions involving one quark-gluon and 

one thrcc-gluon vertex. 
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Consider, next, the one loop contributions to the quark-gluon vertex involving one quark- 
gluon vertex and one three gluon vertex. There are two time ordering contributions shown 
in Figs. 3(a) and 3(b). The contribution from Fig. 3(a) is 
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Here k[ = p[ — k^^ and k 2 = p 2 — k^^. Divergent contributions that 
contain products of logarithms are 
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The boost invariance of this result is again clear. Divergent contributions that contain single 
logarithm are 
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Expressing the terms inside the square bracket in terms of the internal momenta we get 
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which makes boost invariance explicit. 



The contribution from Fig. 3(b) is 
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Here k^ = p^ - k[ + ' X) and k^ = g(+.-D - k{ + ^ . 

Divergent contributions that contain products of logarithms are 
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which is manifestly boost invariant. Divergent contributions that contain single logarithm 
are 
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As in the case of eq.( [5~T3"l) , an addition and subtraction of -4p- in the square bracket, renders 
itself as a combination of canonical vertex and the quantity already encountered in eq. fl5.16|) 
and thereby manifestly boost invariant. 



5.3. COUPLING CONSTANT RENORMALIZATION 



For the sake of completeness, we present here the results for the other diagrams which 
are relevant for the coupling constant renormalization. We also calculate the /3-function 
which exactly matches with the well known results and therefore extends the results arrived 
at in the Ref. |5|, to the most general kinematics in the two-component formalism. 

The sum of divergent contributions from Figs. 2(a) and 2(b) is 
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where we observe the emergence of the canonical vertex structure. 

The sum of divergent contributions from Figs. 3(a) and 3(b) is 
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where we again observe the emergence of the canonical vertex. 

The diagrams in Figs. 4(a), 4(b), 5(a) and 5(b) correspond to the renormalization of 
the external quark and gluon legs that are connected to the vertex. Their contributions are 
given below. 
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Fig. 4 Contribution to the quark-gluon vertex from quark wave-function renormalization. 
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Fig. 5 Contribution to the quark-gluon vertex from gluon wave- function renormalization. 
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Now, to evaluate the contributions to the coupling constant, we have to multiply V4 and 
V5 with I in order to take into account the proper correction due to the renormalization of 



initial and final states |I3| . Thus adding the contributions we get, 
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Note that all the mixed divergences cancel. The correction to the coupling constant is given 
by 
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We compute the /3-function as 
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which is well known result to the one-loop order. 

Running of the QCD coupling can be depicted as follows. Notice that the regularized 
coupling constant depends on the ultraviolate cut-off A as well as fi : 



g R = g R (A,fj) 
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To obtain the renormalized coupling, we have to remove the cutoff dependence by adding 
suitable counter term and then take A —> 00. But we must notice that one energy scale 
/1 has crept in the coupling constant and should be kept there with large enough value to 
ensure the validity of the perturbative calculation. The counter term is chosen to be the 
following : 
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Notice that the counter term introduces another arbitrary energy scale fi r which is large for 
obvious reason, but need not be the same as /1. We call /i r the renormalization scale. So, 
we get the renormalized coupling constant as, 



90 



0ren(/J,A*r) = A lim [fl^A, fi) + 0ct(A, fj, r )] 



4 



8tt 2 V 6 3 S ! > fi 



(5.12) 



Thus the renormalized coupling constant at some scale \x always depends on the renormal- 
ization scale as well as its value at that scale obtained by setting \i = fj, r in eq. (|5.12| ) as the 
renormalization condition: 
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Therefore, from eq. (|5.12|) and eq.( |5.13| ) we get the renormalized coupling constant at any 
perturbatively large scale Q as 



g(Q) = <?ren(Q,/V) = 9 fa [l + ^^(^ - \NjTj) In ^ 



Q 



g{nr 



or, in terms of a s (Q 2 ) = 9 ^ ^ we have 



(5.14) 
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(5.15) 



Expressions in eq.(5.14) and eq.( |5.15| ) are the well known results for the running QCD 
coupling constant. 



5.4. DISCUSSION 

In this Chapter, we have presented the calculations necessary for coupling constant 
renormalization in light-front QCD with the most general kinematics, which is the first 
step towards extending the structure function calculation to fourth order. This calculation 
and also other similar ones employing time (x + ) ordered perturbative techniques in light- 
front theory are known to be straightforward but long and tedious. A lot of effort has 
to be invested in the calculation of coefficient functions accompanying the divergences for 
individual diagrams. No clue seems to be there in the intermediate steps regarding the 
correctness of the calculations due to lack of covariance. Other than actually performing 
the calculations, we tried to investigate in parallel the role played by the Gallilean boost 
symmetry present in light-front formulation here. To the best of our knowledge, this is 
the first work to investigate the utility of Gallilean boost symmetry in determining the 



91 



correctness of the structure of the coefficient functions accompanying the divergences in 
light-front perturbation theory beyond the tree level. 

In this initial investigation we have employed the simplest choice of regulators that 
cutoff single particle momenta. One should note that in addition to possible violations of 
boost invariance, such simple minded cutoff procedure could in principle even introduce 
non-analyticities in the structure of counter terms (see Sec. VI of Ref. H for an explicit 
example). However, in the case of vertex diagrams, we encounter only logarithmic transverse 
divergences. Even with finite cutoffs, violations of transverse boost invariance can appear 
only inside the logarithms and we expect the symmetry to be present in the non-trivial 
structure of the coefficient functions that accompany the divergences. We are primarily 
interested in understanding the complex structure of these coefficient functions on the basis 
of Gallilean symmetry. Incidentally we note that, in contrast, longitudinal boost invariance 
is a scale invariance in light-front theory. The implication of longitudinal boost symmetry 
for the coefficient functions is trivial, namely, simple scaling behavior. 

Let us summarize our findings. Out of all the x + -ordered diagrams relevant for our 
calculation, four involve wave-function renormalization correction and have the structure of 
the canonical vertex. For the remaining diagrams which correspond to vertex corrections, 
the divergent contributions from each of the them contain terms that involve (I) product 
of logarithms and (II) single logarithm. For contributions that belong to (I), we find that 
for each diagram separately, the coefficient of the divergent factor is proportional to the 
canonical vertex and hence Gallilean boost invariance is manifestly maintained. For contri- 
butions that belong to (II), for each diagram, the coefficient of the divergent factor is not 
proportional to the canonical vertex. Nevertheless, in each case, rewriting the coefficient 
in terms of the internal momenta explicitly shows that the coefficient is independent of the 
total transverse momentum P x . Hence for the contributions that belong to (II) the con- 
straint from transverse boost invariance is maintained, even though the canonical form is 
not reproduced. 

Our results show that two-dimensional Gallilean invariance which is manifest at tree level 
is also exhibited in the coefficient functions accompanying the divergences in the regulated 
theory at the one loop level in the case of quark-gluon vertex in light-front QCD even with 
a regulator that violates the symmetry. Since the symmetry is only a part of the complete 
Lorentz symmetry, we expect the constraints which follow from the invariance to be less 
restrictive. Indeed, our results show that the structure of the vertex that satisfies transverse 
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boost invariance is not unique. 

Even-though the canonical vertex structure is not reproduced in the coefficient of the 
single logarithms, it still has some usefulness in practical calculations since it obeys constraint 
from Gallilean boost invariance. The coefficient functions accompanying single logarithms 
are obtained after isolating the leading double logarithms and they exhibit a complicated 
structure. It is quite easy to make a mistake in the sign in one of the terms for individual 
x + ordered diagrams. Our calculations show that using the underlying transverse boost 
symmetry one can easily recognize the mistake in the calculation and hence correct it. 

Finally, of course, we summarized the results for the complete set of diagrams con- 
tributing to the coupling constant renormalization for the massless quark case. We have 
extracted the /3-function and obtained the running QCD coupling renormalized at a partic- 
ular energy scale, both of which match with the well-known results and therefore extends 
the results arrived at previously in the literature to the most general kinematics. Using 



the two-component representation |10] we have presented for the first time the results sep- 
arately for each x + ordered diagram with arbitrary external momenta which is essential to 
study the renormalization of the helicity-non flip parts of the vertex. Present calculations 
together with the calculations presented in Ref. II explicitly show that linear divergences 
of the type - where e is the cutoff on longitudinal loop momenta occur in individual time- 
ordered diagrams only in radiative corrections to the chiral symmetry breaking part of the 
quark gluon vertex. This divergence is a special feature of non-abelian gauge theory. At one 
loop level, this divergence cancels with our choice of regulators when different time-ordered 
diagrams are summed up. Since intermediate states involved are, in general, different in 
different time ordered diagrams, the cancellation may no longer be operative once more so- 
phisticated regulators that explicitly depend on the intermediate states are employed. This 
needs to be investigated in detail in the future because of its nontrivial consequences for the 
renormalization of chiral symmetry breaking terms in the QCD Hamiltonian. 

As far as the structure function calculation to the fourth order goes, one has to essentially 
embed the coupling constant renormalization calculation in another loop, making the final 
calculation a two-loop one. Which seems to be a straightforward extension of the calculations 
presented here (which are all one-loop); but one has to pay proper attention to the intricacies 
involved in the complete two-loop perturbative calculations. And hopefully it would yield 
the solution of Altarelli-Parisi equation to the fourth order in coupling providing greater 
validity of calculating structure functions using our approach. 
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Chapter 6 



In this Chapter, we try to investigate the nucleonic spin structure in the light-front QCD 
Hamiltonian formulation that we have been using throughout to describe DIS structure 
functions. Like the longitudinal momentum, it is of importance to know how the spin of the 
nucleon is shared among its various constituents. Individual constituents are allowed to have 
orbital motions and hence orbital angular momenta which then combine with their intrinsic 
angular momenta (or, spin) to give a net total angular momentum (or, simply the spin, if the 
orbital angular momentum quantum number is zero) for the nucleon. As of now, we can only 
measure the intrinsic helicity part carried by the fermionic constituents (namely, the quarks 
or antiquarks), through the measurement of polarized structure function g\. The orbital 
helicity parts or even the intrinsic helicity for the gluons are not known how to be measured 
in the experiments. Since helicity is the object measured in the experiment, we shall only 
be concerned with the third component of the angular momentum operator throughout this 
Chapter. Thus, as it turns out, the contribution from constituents' orbital helicity to the 
total helicity of the nucleon is more of theoretical interest at present. In fact, it is still 
theoretically uncertain how the total helicity operator (relevant for nucleon), as defined in 
the field theory, separates into a sum of corresponding pieces describing intrinsic and orbital 
helicity of the constituents. The problem seems to lie in the gauge invariance. The quark or 
gluon fields in a particular gauge becomes a combination of both by a gauge transformation. 
So, an operator defined to be, for example, the quark orbital helicity operator solely in terms 
of quark fields in a particular gauge loses its characteristics in another gauge. Thus, the 
separation we are interested in seems to be unrealistic due to lack of gauge invariance. We 
shall discuss this problem in detail in Sec. 6.1, which hints towards the resolution that the 
gauge fixing may be necessary in order to define such operators. There we shall derive such a 
gauge fixed (in A + = gauge) light-front operators, which necessarily does the separation we 
are looking for. Next, in Sec. 6. 2, with the help of these operators we shall define the relevant 
structure functions, which can measure the orbital (and intrinsic) helicity contributions of 
quarks or gluons to the total helicity of the nucleon. Then we calculate these structure 
functions for a dressed parton target, which shows the utility of our definitions by obtaining 
various anomalous dimensions necessary to study their Q 2 -evolution in Sec. 6. 3. Validity of 
our calculations are checked by showing that the helicity sum rules are satisfied in Sec. 6. 4. 
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Thus, the plausibility of defining such structure functions and the calculations leading to 
various anomalous dimensions in a transparent way (as will be discussed here) exhibit one 
of the major triumph of our formulation presented in this dissertation. 

6.1. LIGHT-FRONT HELICITY OPERATOR J 3 FROM THE MANIFESTLY 
GAUGE INVARIANT ENERGY MOMENTUM TENSOR 

A. A Brief Review 

Before we show the detailed derivation of the light-front helicity operator J 3 , let us 
first review the situation which necessitates such an exercise. The well-known proton spin 
crisis that emerged after the the publication of EMC data on the measurement of polarized 
structure function gi(x, Q 2 ), drew a lot of theoretical interest on how the total helicity of the 
nucleon is distributed among its constituents. In the parton model, partons are assumed to 
be moving collinearly (i. e., with zero relative transverse momenta kx = 0) and do not have 
any orbital motion L = f x k^ = 0. Thus, the total helicity of the nucleon comes solely from 
the intrinsic helicity of the constituents (quarks and gluons). Now an integral of gi(x,Q 2 ) 
over all possible x can be shown (see later) to measure, out of the total nucleonic helicity, 
only the intrinsic helicity part (AS) coming from the fermionic constituents. EMC results 
showed that this AE contribution was very small, which appeared to be quite puzzling and 
became known as the proton spin crisis. Incidentally, the role of orbital angular momentum 
in deep inelastic scattering was first emphasized by Sehgal |]] and then by Ratcliff [[| in 
the context of Altarelli-Parisi equation for real partons in QCD. But it should be noted 
that, in the interpretation of EMC data, orbital helicity parts of the constituents were 
ignored as some higher twist contribution, even if the partons in reality might have some 
non-zero transverse momenta kx ^ 0. On the other hand, it is understood that an anomaly 
contribution coming from the gluonic sector is responsible for a cancellation to occur, giving 
rise to the small measured value of AS. It also made ambiguous whether to call AE as the 
intrinsic helicity contribution from quarks and anti-quarks or the one before the cancellation 
occurred. Whichever be the situation, it is now realized that to understand the nucleonic 
spin structure properly, one should put more emphasis in studying the helicity sum rule, 
which anyway has to be satisfied and free from anomaly due to the conservation of angular 
momentum. 

The total helicity operator J 3 acting on a nucleon state | PS) gives, 
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J 3 \PS) = ±~\PS), (6.1) 



depending on the value of 5 = ±~. Thus, we have, 



1 



'q(i) -r "q(o) -r J g(i) -r ^g(o) 



1 



<^| ^m + ^ + ^m + ^o) l^> = ±-, (6.2) 



2 



with A/" = (PS | PS 1 ). Here we have assumed for the time being that the total J 3 can be 
separated into corresponding intrinsic and orbital helicity operators for quarks and gluons. It 
is this separation which becomes ambiguous due to lack of gauge invariance of the individual 
parts. 

Before we proceed further, let us give a brief account of recent developments in this 
regard. Jaffe and Manohar || first noted that angular momentum operator constructed from 
the gauge invariant, symmetric, energy-momentum tensor fails to display the distinction 
between intrinsic and orbital angular momentum. They suggested the use of free field theory 
form which is interaction independent and can be separated unambiguously into quark and 
gluon orbital and spin parts. Ji, Tang and Hoodbhoy [f|], starting from the Jaffe-Manohar 
choice, studied the asymptotic fraction of the nucleon spin carried by quarks and gluons at 
the one loop level. They mentioned that this separation is gauge variant and supported the 
choice of light-front gauge and light-front coordinates for their calculation. Later on Ji || 
introduced a gauge invariant definition of J q and J g starting from gauge invariant symmetric 
energy-momentum tensor. But these are now interaction dependent contrary to the well- 
known kinematical nature of the angular momentum operators (only J 3 in light-front theory 
is so). No justification has been given why they are called angular momentum operators. 

In fact, Singleton and Dzhunushaliev || claim to show by explicit calculations that the 
gauge invariant orbital angular momentum operator proposed by Ji || do not obey the 
angular momentum algebra and hence do not qualify as the angular momentum operator. 
It seems to be in agreement with the observation made by Chen and Wang |7[ regarding 
the gauge invariant definition (that they do not obey the angular momentum algebra) and 
they further put more emphasis on the matrix elements, which are important as far as 
the experiment is concerned. From their argument it appears that the gauge dependent 
operators may have gauge invariant expectation values in the hadronic eigenstates having 
definite angular momentum. There are other works (for example, see Refs. 0,0, [IC|, [IT 



|IH?|) which show that the issue of gauge invariant separation is still under hot debate. In our 



case, we always work with the gauge fixed theory which is free from any of these ambiguity 
as we discuss next, in detail. 
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B. J 3 in light-front gauge 

It is well-known that the energy-momentum density (which gives rise to Hamiltonian 
and three-momentum) and the generalized angular momentum density (which gives rise to 
angular momentum and boosts) can be expressed in a manifestly covariant, gauge invariant 
form. But, in order to define the Poincare generators in quantum field theory, one has 
to choose a particular hypersurface over which these densities are integrated. Thus, if 
one chooses x° = as the surface over which integration is to be performed, one gets 
using standard notations, P^ = Jd 3 x6 0fl . On the other hand, the choice x + = gives 
P^ L = J dx~d 2 x ± 8 +tl . Therefore, the the Poincare generators explicitly depends on the 
frame of reference. It may not be surprizing that it also depends on the gauge choice in the 
case of gauge theory. This of course does not imply that the theory has lost Lorenz and 
gauge symmetry. The symmetries are no longer manifest, but the physical observables in 
the theory still obey the consequences of the symmetries. 

Poincare generators can be further classified as kinematical (which do not contain interac- 
tions and do not change the quantization surface) and dynamical (which contain interactions 
and change the quantization surface). Which operator is dynamical and which is kinematical 
of course depends on the choice of quantization surface. It is well-known that in light-front 
field theory, on which our formalism of deep inelastic scattering is based on, the generators 
of boosts and the rotation in the transverse plane (light-front helicity) are kinematical like 
three momenta whereas the generators of rotations about the two transverse axes are dy- 
namical like the Hamiltonian. Thus, the operator in light-front field theory relevant to the 
proton spin crisis is the light-front helicity operator which belongs to the kinematical sub- 
group. In light-front literature, it is customary to construct this operator from the canonical 
symmetric energy momentum tensor and one explicitly finds that this operator is indeed 
free of interaction and has the same form as in free field theory |13| . 



In non-Abelian gauge theories like QCD, one should be extra cautious since such theories 
are known to exhibit non-trivial topological effects. In this work, we restrict our attention to 
the topologically trivial sector of QCD which is relevant for DIS. In this sector, interactions 
do not affect kinematical generators [13 j. In view of the prevailing confusion in the literature 



as we mentioned, we provide an explicit demonstration of this fact in this section in the case 
of the light-front helicity operator. 

We start from the manifestly gauge invariant, symmetric energy momentum tensor in 
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QCD: 



&»= -i,\YD u + YD"]ip - F^F^ 



- g^{ - \{F Xaa ) 2 + ^ X D X - m)ij}, (6.3) 

where iD* = id 11 + gA», F» Xa = <9M Aa - d x A^ a + gf^A^A^ , F u £ = d ua A x - d x A ua + 
gf abc A ub A c x . 

We define the light-front helicity operator 

J 3 = - fdx-d 2 x ± [x 1 Q +2 -x 2 Q +1 ]. (6.4) 

J 3 is a manifestly gauge invariant operator by construction. However, it depends explicitly 
on the interaction through @ +J and does not appear to be a kinematical operator at all. 
Furthermore, it is not apparent that J 3 generates the correct transformations as an angular 
momentum operator. Thus at this stage, we are not justified to call it a helicity operator. 
Explicitly, we have, 



J 3 = - f dx~d 2 x ± {x 1 [-^{ 1 + D 2 + 7 2 D + )^ - F +Xa F 



2al 

aJ 



- x 2 [-ij^D 1 + 7 1 J D+)^ - F +Xa F^] } . (6.5) 

Notice that the last term in + * does not contribute for g +l = 0. The fermion field can be 
decomposed as usual, ■?/ ,± = A^ifj, with A 1 * 1 = |7 =F 7 ± and we shall work in the gauge A + = 0. 
In this gauge, we still have residual gauge freedom associated with ^"-independent gauge 
transformations. Note that only ip + and A 1 are dynamical variables whereas ip~ and A" are 
constrained. 

Let us proceed to calculate +2 explicitly in terms of dynamical fields. We have, 

-^( 7 + D 2 + 7 2 J D + )V> = ^ +t id 2 ^ + + gtp +t T a tl) + A 2 a + -^ 7 2 i<9 + ?/>, (6.6) 

where the last term depends on the interaction through the constrained field ip~ and the 
interaction dependence in the second term is explicit. We use the constraint equation 

id + ip- = [a 1 - ■ (id 1 - + gA L ) + 7 °m]^ + , (6.7) 

to eliminate the constraint variable ip~ as follows. 
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-^ 2 d + ^= -^ +t 7 V(9 + ^-) + -V- t 7° 7 2 (9 + ^ + ) 



-^a 2 (d + ^) - -{d + ^)a 2 ^j + + -d+(>~W + ) 



2^ 



a 



(id + gA L ) + 7°m ip + + -ip +j 



a 



{-id + gA L ) + 7°m| a 2 ^ 



+ V(fW + ) 



-?/> +t aV(<9^ + ) - -(d^K'crVi + -^ +t (a 2 a J + a j a 2 )Ty + A ja 



2 



+ V(^"W + ) 



^ +t 9 2 V + + -9 1 (V' +t E 3 V' + ) + # +t T>+A 2a 

+ V(fW + )-V(^V) 



(6.8) 



Notice that d +)2 are like total derivatives (but not d 1 ), since G +2 term is multiplied by x 1 
and we have used this fact to switch around <9 + ' 2 producing surface terms (see the last line 
in eq. (|6.8|) ) . First and second terms in the last line follow from the first term in the previous 
line for j '• = 2 and j '• = 1 respectively with a corresponding surface term. Also note that the 
terms involving m has canceled and we have used S 3 = ia 2 a l . 

Now we restrict ourselves to the topologically trivial sector by requiring that the dynam- 
ical fields (ip + and A 1 ) vanish at x~' 1 — > oo. The residual gauge freedom and the surface 
terms are no longer present and so we drop total derivatives of d + and d 2 . Thus we obtain, 

-^(7+D 2 + 7 2 J D+)V = 2iV> +t 9 2 V + + -9 1 (^ +t S 3 ^ + ) + 2g^ +] T a ^ + A 2a . (6.9) 

Notice that the last term explicitly depends on the interaction and rest of the expression is 
the same as one would have got in the free theory. 

Similarly, we can calculate the gluonic contribution to + \ In the gauge A + = 0, we 
have, 

_ F +\a F 2a = _( d +A Xa )(d 2 A a x - d x A 2a + gf abc A 2b A c x ) 

= -{d + A~ a ){d + A 2a ) + d+A ja {d 2 A ja - d j A 2a ) + gf abc (d + A ja )A 2b A jc 

= --(d + ) 2 A~ a A 2a + d + A ja (d 2 A ja - d j A 2a ) + gf abc (d + A ja )A 2b A jc 



V-d + [d + A- a A 2a 



(6.10) 



where we have done partial integration using d + . This still depends on the constrained field 
A~ a . We have the constraint equation for the elimination of the variable A", 
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-{d + ) 2 A- a = d + d i A m + gf abc A ib d + A ic + 2gi) +] T a i) + . (6.11) 



Thus, we obtain, 



_ p+\ap2a = d + A J^ d 2 A Ja _ gj jtfaj + (& A ja )(d + A 2a ) - 2g^T^~ 

+ -d + (d + A~ a A 2a ) - d + (d l A ia A 2a 
= d l A ia d + A 2a + d + A ja (d 2 A ja - d j A 2a ) - 2gi) +] T a i) + A 2a 



+ -d + (d + A- a A 2a ) - d + (d j A ja A 2 



\2a^ 

I i ( / . -i ,-i i ■ ■ ( / i ( / . -i • . - 

= 8 + A la d 2 A la + 8 + A 2a d 2 A 2a + d\A la d + A 2a ) - 2gi) +] T a i) + A 2a 

+ ]-d + (d + A- a A 2a ) - d + (d J A' a A 2a ) - d + (A la d L A 2a ) . (6. 12) 

Again we have performed appropriate partial integrations. Note that the terms correspond- 
ing to the interaction among the gluon fields got cancelled already. As in the earlier case, 
we drop the surface terms involving d + . Thus, collecting together the results in eq. (|6.9| ) and 
eq.(|6.12|), we get, 

9 +2 = 2# +t <3> + + V(^ +t £ 3 ^ + ) 

+ d + A la d 2 A la + d + A 2a d 2 A 2a + d l (A la d + A 2a ) , (6. 13) 

where we see that the interaction dependent terms get cancelled completely among them- 
selves. By a similar calculation, 

+1 = 2«V +t 0V + - -d 2 (ip +t X 3 iP+) 

+ d + A la d l A la + d + A 2a d 1 A 2a + d 2 {A 2a d + A la ). (6.14) 

From the above two equations it is clear that G +1 and +2 agree with the free field the- 
ory form at the operator level. This shows that in light-front quantization, with A + = 
gauge, J 3 = J 3 (the naive canonical form independent of interactions) at the operator level, 
provided the fields vanish at the boundary. Using the form of +4 in eq.( |6.5| ), we obtain 

J = J f(o) + J f(i) + J g(o) + J g(t), (6-15) 

with 

73 _ f J„.-j2„.±„/,H-t„Y~lp2 _ „2ol 



jj {o) = / dx-d 2 x ± ^ +l t(x 1 d 2 -x 2 d 1 )^ 
J z m =\jdx-d 2 x^ +] ^ + , 
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JJ (o) = -Jdx~d 2 x ± {x 1 [d + A 1 d 2 A 1 + «9 + A 2 <9 2 A 2 ] - x 2 [d + A l d l A l + d+A 2 ^ 2 ]}, 
J 3 g{t) =^fdx-d 2 x ± [A 1 d + A 2 - A 2 8 + A 1 ]. (6.16) 



The colour indices are implicit in these equations. 

Using canonical commutation relations, we can explicitly find that, 



(x 1 d 2 -x 2 d 1 )ij + (x), 
1 



T 3 


ip + (x) 


T 3 


ip + (x) 


T 3 

J g{o) 


,A\xj 


T 3 


,A\xj 



1 1 i> (x), 
(x 1 d 2 -x 2 d 1 )A l {x), 



(6.17) 
(6.18) 
(6.19) 
(6.20) 



Thus, these operators do qualify as angular momentum operators (generators of rotations 



in the transverse plane) in the theory [13 



To summarize, the helicity operator constructed from manifestly gauge invariant, sym- 
metric, energy momentum tensor in QCD, in the gauge A + = 0, and after the elimination 
of constraint variables, is equal to the naive canonical form of the light-front helicity oper- 
ator plus surface terms. In the topologically trivial sector, we can legitimately require the 
dynamical fields to vanish at the boundary. This eliminates the residual gauge degrees of 
freedom and removes the surface terms. Thus we have a gauge fixed Poincare generator 
which we consider in the following sections. 

6.2. ORBITAL HELICITY DISTRIBUTION FUNCTIONS 

In this section, we use the gauge fixed helicity operator J 3 and its physical separation 
as shown in the previous section and define various structure functions containing the in- 
formation regarding nucleonic helicity structure. The polarized structure function gi(x, Q 2 ) 
is connected to the intrinsic light-front helicity content of the quarks and anti-quarks as 
is evident from the following expression. We recall the expression for g\ given earlier in 
Chapter 3, 

1 



gi(x,Q 2 



8ttS- 



d V e-^^e 2 a (PS\ ^ a (r)7 + 7 5 V<a(0) + (£ ~ 0) | PS), (6.1 



with 7] = |-P + £ • By following exactly the same procedure as that used in obtaining the 
sum rule for F 2 , we get, 
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l dx 9l (x,Q 2 ) = -L^e 2 a {PS | ^t (0)S 3^ (0) | ps) 



(6.2) 



where S 3 = ij l j 2 . 

One can define the intrinsic helicity distribution function Aq(x, Q 2 ) to be the same as 
gi with the weight factors e a — 1, 



Aq(x,Q 2 



1 



8nS+ 
It follows from eq.( |6.2|) that 



dr]e- irix (PS | ^(D7 + 7 5 ^(0) + (f <-> 0) | PS). 



(6.3) 



1 da;Ag(x,Q 2 )= i-(PS | ^ +t (0)£ 3 ^+(0) | PS) 
o 2o + 

(PS | J 3 (t) | PS) 
2S+(2tt) 3 <5 3 (0) 



(6.4) 



where A/" = 2(27r) 3 P + 5 3 (0). In the last step, we have taken into account the fact that 
S + = ±P + for longitudinally polarized target (only in which case the above formula makes 
sense) . 

In analogy with eq. (|6.3| ), we define the orbital helicity distribution for the fermion and 
gluons as follows: 

1 



Aq L (x,Q 2 
Ag L (x,Q 2 



4ttP h 
-1 

4ttP^ 



d7]e- lvx (PS | ^(^-)7 + z(ar 1 ^ - x 2 d l )ip(Q) + h.c.\ \ PS) 
drie- ir > x (PS\ \x 1 F +a (C)d 2 A a (0)-x 2 F +a (C)d 1 A a (0)} \ PS) 



(6.5) 



Here | PS) denotes the hadron state with momentum P and helicity S. Similarly one can 
define the gluon intrinsic light-front helicity distribution [ |15| 



as 



Ag{x,Q 2 



4tt(P+) 2 x 
The dual tensor is given by 



dT 1 e~ irix {PS | F +a (C)F + a (0) | PS). 



(6.6) 



F"" = -e» vpa F pai with e +1 ~ 2 



2. 



(6.7) 



Note that the above distribution functions are defined (as we always do) in the light-front 
gauge A + = 0. 
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The above distribution functions are denned in such a way that integration of the above 

distribution functions over x is directly related to the expectation values of the corresponding 

helicity operators analogous to eq.(|6.4j) as follows: 

1 



dxAq L (x,Q 2 )- — r 
o Jv 



(PS | JJ o) | PS) 



£dxAg{x,Cf)= ^(PS\J 3 g{l) \PS) 



dxAg L (x, Q 2 



1 



(PS | J 3 g{o) | PS) 



(6.8) 



/o M 

Eq.( |6.5|) constitutes the new definition we proposed for the structure functions containing 
the information regarding orbital helicity of quarks and gluons. Now, in the next section we 
aim to demonstrate the utility of these definitions by explicit calculations. 



6.3. PERTURBATIVE CALCULATION OF ANOMALOUS DIMENSIONS 



Before showing the explicit calculations, we first concentrate a little more on the vari- 
ous gauge fixed angular momentum operators as defined in eq.( |6.16| ). We use the Fourier 
decomposition of the dynamical quark and gluon fields as given in Chapter 2, 

dk + d 2 k J 



and 



ip (x) 



A ± (x) 



2(27r)Vfc+ 

dk + d 2 k x 
2(2vr) 3 fc+ 



Y,Xs[b(k, 



s)e 



-ik.x 



+ (ft(k, 



-s)e 



Ik.x 



ERm) 



e^e 



-ik.x 



+ a\k,X)(ei)*e 



* ik.x 



(6.1) 



(6.2) 



in terms of which the J 3 s become as follows: 



J f{o)~ l Y. 



T 3 



*E* 



dk + d 2 k A 
2(2nfk- 
dk + d 2 k 



b\k 



k 2 







k 



i d 



b(k,s) + d t (fc, 



A 



J 



A 



2(2ir) 3 k+ 
dk + d 2 k 



dk 1 dk 2 

tf(k,\)b(k,\) + S(k,\)d(k,\) 



k 2 L - k' d 



dk 1 



dk 2 



d(k, s) 



7 3 " 



EA 



2(2tt) 3 ^ 
dk + d 2 k~ 



a\k,\) 



k 2 JL-k^ 9 



dk 1 



dk 2 



a(k, A), 



■a\k, X)a(k, A). 



(6.3) 



A . 2(2 V 

In actual calculation we generally replace the momenta of individual quarks and gluons (in 
the above operators) in terms of relative internal momenta and centre of mass momenta. 
We first show how special features of light front transverse boost symmetry simplifies the 
calculation and helps understanding some of the features of the calculation very easily, which 
had been claimed to be very surprizing in the literature H. 
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A. Internal orbital helicity: Non-relativistic versus light-front (relativistic) case 



Here we address the connection between non-relativistic situation and the light-front 
relativistic case. We need to decompose the total orbital angular momentum of a composite 
system as a sum of the orbital angular momentum associated with internal motion and the 
orbital angular momentum associated with the center of mass motion. We are interested only 
in the former and not in the latter. For illustrative purposes, consider a two body system 
consisting of two particles with masses mi and 7712 and momenta ki and k2- Let P denote 
the total momentum. In the non-relativistic case, let q denote the relative momentum, i.e., 
_ m 2 ki-r?iik2 ^ j.j- j g we h_] inown I|i7|] that the contribution of particle one (two) to the third 

^ mi +7712 " u v ' 

component of internal orbital angular momentum is given by 



r 3 

^1(2) 



; m 2(l) 

mi + m 2 



Q 



d , d 

q 



(6.4) 



dq 1 "* dq 2 - 

Next consider the light-front case. Let k\ = (kf, k x ) and k 2 = (k 2 , k 2 ) denote the single 
particle momenta and P = (P + , P^) denote the total momentum of the two particle system, 
i.e., fei +&2 = P + '\ Light-front kinematics allows us to introduce boost-invariant internal 
transverse momentum q L and longitudinal momentum fraction Xi by 



k l L = q ± + x 1 P ± , k+ = Xl P + , k 2 L = -q ± + x 2 P ± , k+=x 2 P + . 



(6.5) 



Note that x\ + x 2 = 1 and q L = x 2 k x — X\k 2 . For the first particle, we have 



i_d_ ,1 9 
dk{ l dk 2 
d 



ix 2 



Q 



, ,, i d 



dq 1 

+ %X\X 2 



Q 



dq 2 
d 



+ ix\ 
,1 9 



P' 







P 



1 d 



P 2 — - P 

dq 1 dq 2 



dP 1 dP 2 - 

d x d 



+ t 



Q 



dP 1 



q 



dP 2 



For the second particle, we have 



rtn ~ : -, rv- 



2 dk\ 2 dk 2 



IX 1 



2_^__ 1_^_ 

dq 1 dq 

d 



— IX\X 2 



+ ix 2 
,1 d 



P< 



d 



op 1 



p 2 — — p 

dq 1 dq 2 



— 1 



-P 

, 9 

op 1 



1 d 



dP 2 



JL 9 



OP 2 



Total orbital helicity 



(6.6) 



(6.7) 
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t . = i ^ = j[ ^_^] +j[ P^_^]. (6.8) 

Thus we have decomposed the total orbital helicity of a two particle system into internal 
orbital helicity and the orbital helicity associated with the "center of mass motion" . 

Note that the internal orbital helicity carried by particle one is the total internal helicity 
multiplied by the longitudinal momentum fraction carried by particle two and vice versa. 
This factor can be understood by comparison with the situation in non-relativistic dynamics 
and recalling the close analogy between Gallilean relativity and light-front dynamics in the 
transverse plane. In non-relativistic two-body problem, the center of mass coordinate is 
defined by R= m ^ 1 l™ 2r2 • The generator of Gallilean boost is B— — J2i m i r %- Thus in non- 
relativistic dynamics, R— — j| with M = mi + m 2 . In light-front dynamics, the variable 
analogous to B L is E^, the generator of transverse boost and the variable analogous to M 
is P + . Thus in light-front theory, the transverse center of mass coordinate R = y, ' + ' = 
x\Ti + X2T2 ■ Thus we recognize that instead of m "T m2 ( m"+m ) * n n on-relativistic theory, 
xi (xi) appears in light-front theory. (See the the discussion at the end of this Chapter how 
this helps understanding our result very easily.) 

By comparing light-front (relativistic) and non-relativistic cases, we readily see that the 
role played by particle masses in individual contributions to the third component of internal 
orbital angular momentum in non-relativistic dynamics is replaced by longitudinal momen- 
tum fractions in relativistic (light-front) theory. This also shows that the physical picture 
of the third component of internal orbital angular momentum is drastically different in non- 
relativistic and relativistic cases. We stress that it is only the latter, in which parton masses 
do not appear at all, that is of relevance to the nucleon helicity problem. Lastly, we empha- 
size that it is the transverse boost invariance in light front dynamics that makes possible 
the separation of dynamics associated with the center of mass and the internal dynamics. 
In equal-time relativistic theory, this separation cannot be achieved at the kinematical level 
since boosts are dynamical. 

B. Dressed Parton Calculations 

In this section, we evaluate the internal helicity distribution functions for a dressed quark 
in perturbative QCD by replacing the hadron target by a dressed quark target. We have pro- 
vided the necessary details of the calculation which may serve as the stepping stone for more 
realistic calculation with meson target. From this simple calculation, we have illustrated how 
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easily one can extract the relevant splitting functions and evaluate the corresponding anoma- 
lous dimensions. Note that, since we are not interested in exhaustive calculation of various 
anomalous dimensions and the purpose of this section being illustrative, we can safely drop 
the derivative of delta function which will come naturally in the following calculations and 
work explicitly with forward matrix element. 

The dressed quark state with fixed helicity can be expressed as 

\k\ k ± , a) = **(*)&! (^lo) + £ / dkt<fk ^ d -tj£hg ^^ 6 -\k -h- k 2 ) 

x $ x XiX2 (k; k u h^U^alihm + ■ • • , (6.9) 

where the normalization of the state is determined by 

(k' + , k' ± , X'\k + , k±, A) = 2(27i) 3 k + 5 x ,yS{k + - k' + )5 2 (k ± - k' ± ), (6.10) 

We introduce the boost invariant amplitudes ip* and ip* Xi (x, k x ) respectively by $ A (/c) = ip* 
and Q XiX2 (k; ki, k 2 ) = -7=?/^ lA2 (x, fir 1 ). From the light-front QCD Hamiltonian, to lowest 
order in perturbation theory, we have (see Appendix C for details), 



./57^i VT^^+mm-xY^\ i 



w / 2(2vr) 3 VI -xk? ± + m 2 q (l-x) 2 ai { 1-x 



+ -(o-i. • k^)^ - im 3 a*— W*(A 2 ) Vi- (6.11) 

X X J 

Here x is the longitudinal momentum fraction carried by the quark. We shall ignore the 
m q dependence in the above wave function which can lead to higher twist effects in orbital 
helicity. In the following we take the helicity of the dressed quark to be + |. Note, due 
to transverse boost invariance, without loss of generality, we take the transverse momen- 
tum of the initial quark to be zero. Kinematical nature of the transverse boost invariance 
also suggests that the wave-functions ip^ lX2 (x, k±) are independent of the target momenta 
k l (transverse component of centre of mass momenta), as is also evident from eq.( 3.11| ). 



Therefore, in the calculation of structure functions involving orbital motions, only the first 
term in eq.(|6.6|) or, eq.( |6.7| ) contributes. 

First, we evaluate the intrinsic helicity distribution functions Aq(x, Q 2 ). As in the calcu- 
lation of unpolarized structure functions, we use the Fourier decomposition of the dynamical 
fields and use the dressed quark state eq(|6.9j) to work out Aq(x, Q 2 ) given in eq.(|6.3[). Ex- 



plicit calculation using the standard commutation relation for creation and annihilation 
operators of dynamical quarks and gluons fields, gives 
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Aq(x,Q 2 ) = {- | 4 | 2 8(1 -x) + Y, *ijd 2 ^ | ^ lM {x,K L ) | 2 } 

o"i,A2 



(6.12) 



where we have used S + = k + for A = | in our case. Here and in the following, we have 
replaced the superscript A by f in the wave- functions if)'s, since in our calculation it is always 



l 

"2 - 



In the massless case, using eq. fl6.ll ) with m 

1 9 2 



5>i 

o"i,A2 



d K I V'ii.Aa ( X >« ) 



2(2tt) 



C 



0, we arrive at 

1 + x 2 /■ ,„ I 



/■ 



1 — x 



dV 



W< 



^ll 2 



(6.13) 



where C 



f 



N 2 -l 



, x . It is instructive to compare this equation with the corresponding one 
in unpolarized case. The transverse momentum integral in Eq. ( |6.13| ) is divergent at both 



limits of integration. We regulate the lower limit by \i and the upper limit by Q. Thus we 
have 



x jJ? 



Aq{x,Q 2 ) = -\i,\\ 2 [5{l-x) + 1 ^C f - 



(6.14) 



Note that jjl has to be large enough so that perturbative calculation is not invalidated. Now, 



ipi \ 2 is evaluated using the normalization condition eq. (|6.10| ) and proceeding exactly the 



same way as in the case of unpolarized structure function F 2 9 , we get to order a s 



Aq(x,Q 2 ) = 1 -[5(l-x) + ^C f ln^ 



1 + x 2 3 r/1 . 



(6.15) 



(6.16) 



Notice that 

f 1 Aq(x,Q 2 ) =\. 
Jo 2 

Next, we evaluate the gluon intrinsic helicity distribution function Ag(x, Q 2 ) as given in 

eq. (|6.6|) in the dressed quark state. Using the the Fourier decomposition of the dynamical 

fields and the dressed quark target given in eq. (|6.9| ) it easy to see that the non vanishing 

contribution in Ag(x, Q 2 ) comes from the quark-gluon state. A straightforward evaluation 

gives, 



As(l-x,Q 2 )= £ A 2 j d 2 K^i iX2 \x,K^i iX2 {x,^) 



o-i,A2 

^ln-.r'll 



a, „ , Q 2 2 



x 



x 2 1 



X 



1-x) 



^1 



(6.17) 



2tt j '" fi 2 " v_ "' 1-x 
The first (second) term inside the square bracket arises from the state with gluon helicity 
+1 (-1). To order a s calculation, we take | ipi | 2 to be unity. So, we have the gluon intrinsic 
helicity contribution in the dressed quark state, to order a s , 
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a O 2 

Ag(l-x,Q 2 ) = ^C f \n^(l+x). (6.18) 

Note that the gluon distribution function has the argument (1 — x) since we have assigned 
x to the quark in the dressed quark state. 

Next, we proceed in a same way to evaluate the quark orbital helicity distribution func- 
tion given in eq.( |6.5| ) in the dressed quark state. Again the non vanishing contribution comes 
from the quark- gluon state. We get, 

d 
ni 1.1 i . 



^ C/ln ^ (1 _^ (1 _ l) ^[_l___J_j |^,v (6 . 19) 



The first (second) term inside the square bracket arises from the state with gluon helicity +1 
(-1). Note that for convenience in calculation, we have changed the variable k -1 : (k 1 , k 2 ) — > 
(| /r 1 |,0). Thus we have the quark orbital helicity contribution in the dressed quark state 
to order a s (and therefore, taking | ip\ | 2 = 1), 



Aq L (x, Q 2 ) = -£!C, In *- (1 - x )(l + x ). (6.20) 



2tt t /i 2 

Similarly, we get the gluon orbital helicity distribution defined in eq.( |6.5| ) in the dressed 
quark state as 



J d*K^ X1p^ lX2 {X,K^)(-l 



Ag L (l-x,Q 2 )= J2 j d 2 ^ xi/jl^ix^^i-i—^l^ix^K- 



a s n , Q 2 

— Of In — 

27T /i 



C/ln-^- x{l + x), (6.21) 



where we have already taken \%pi | 2 = 1 here for the order a s calculation. 

We note that the helicity is conserved at the quark gluon vertex. For the initial quark 
of zero transverse momentum, total helicity of the initial state is the intrinsic helicity of the 
initial quark, namely, +| in our case. Since we have neglected quark mass effects, the final 
quark also has the same intrinsic helicity +~ as is also evident from eq.( |6.16| ). Thus total 
helicity conservation implies that the contributions from gluon intrinsic helicity and quark 
and gluon internal orbital helicities have to cancel. This is readily verified using eqs. ( |6. 18[ ) , 
( |Oq ), and flOT] ). 

From eqs. ( |6.18 ), ( |6.2(J ) and (6.21) we extract the relevant splitting functions. The 



splitting functions are 
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Pss( 9q ){l-x) =C f (1 + x) 

PLS {qq) (x) = - C > (1 - 



2^ 



^5(«)(l-x)=-C>x(l + x). (6.22) 

We define the anomalous dimension A n = J dxx n ~ 1 P(x). Thus, we can easily work out the 
set of corresponding anomalous dimensions using eq. ( |6.22j ) and these are given by 



n + 2 A n r 2 ,n __ r ^ + 4 

n ( n + 1) ' Als M 0/ n (n + 2) ' L5(9<?) " J n(n + 1) (n + 2) ' 



It is to be noted that these anomalous dimensions agree with those given in the recent work 



of Hagler and Schafer fll 



6.4. VERIFICATION OF HELICITY SUM RULE 



Helicity sum rule for the fermion target is given by 



1 



>q(i) -r ^ q {o) T ogft -r o g{o) 



1 



(PS\ J (t) + J(o) + J a) + J (o) \PS) = ±-. (6.1) 



2 



For boson target RHS of the above equation should be replaced by the corresponding helicity. 

Here we verify the correctness of our definitions of distribution functions in the context 
of helicity sum rule for a dressed quark as well as a dressed gluon target perturbatively. This 
is necessary as a consistency check of our results as well. We use the operator J 3 s given 
in eq.( |6.3|) and calculate their matrix elements in the dressed quark and gluon states. For 
simplicity, we take the external transverse momenta of the target to be zero so that there is 
no net angular momentum associated with the center of mass of the target. 

Here we present the final results omitting all the details of calculations (see Appendix 
C). For a dressed quark target having helicity +~ we get, 



jfftWn 


PA) q = Jdx 

1 


-5(1 -x)-\ C f In -2- — + -6(1 - 

-2 v ; 2tt / n 2 \l-x)+ 2 v 


-x)\ 


]\f( P ^\ J f(o) 


2 
P,]) q = -^C f \n%Jdx(l-x)(l + x) 




Jf( P ^\ J 3« 


p,T>«= ^ c f l ^J dx ( 1 + x ) 




Jf( P >1\ J 9(o) 


p t) = -—( 


Of In — - / dx x (1 + x). 
jjr J 





(6.2) 
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Adding all the contributions, we get, 

Jf < P > 1 1 J 3 m + 4o) + Jh) + ■&) \PA) q =\- (6-3) 

For a dressed gluon having helicity +1, the corresponding expressions are worked out to be 
the following. 



^\ J m 


\P>l)g = 







^( p , Tl ■//(„) 


\P,l)g = 






^( P >Tl Jg(i) 


\P,1)g = 


</>;</>! 






= 


1-£W*$ /*& + (!- 


- x) 


1 1 3 
Jf\ P ^\ J 9(o) 


I^T>* = 








(6.4) 



Adding all the contributions, we get 
1 



j?(p, Tl J 3 m + J 3 m + J 3 g{i) + 4o) I p, T> 9 = i • (6.5) 



Thus, eq.( |6.3| ) and eq.( |6.5| ) clearly show that the sum rules are satisfied. Note that in 
evaluating the above expression, we have used the Fourier decomposition of the dynamical 
fields and the Fock-expansion for the target states. For the dressed quark we have used 
eq. (|6.9|) , while for gluon we have used similar expansion but ignored two-gluon Fock sector 
for simplicity. 

Before concluding, let us summarize what we have presented in this Chapter. We have 
presented a detailed analysis of the light-front helicity operator (generator of rotations in 
the transverse plane) in QCD. We have explicitly shown that, the operator constructed 
from manifestly gauge invariant, symmetric energy momentum tensor in QCD, in the gauge 
A + = 0, and after the elimination of constraint variables, is equal to the naive canonical form 
of the light-front helicity operator plus surface terms. In the topologically trivial sector, we 
can legitimately require the dynamical fields to vanish at the boundary. This eliminates the 
residual gauge degrees of freedom and removes the surface terms. 

Next, we have defined non-perturbative quark and gluon orbital helicity distribution 
functions as Fourier transform of forward hadron matrix elements of appropriate bilocal 
operators with bilocality only in the light-front longitudinal direction. We have calculated 
these distribution functions by replacing the hadron target by a dressed parton providing all 
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the necessary details. From these simple calculations we have illustrated the utility of the 
newly defined distribution functions in the calculation of splitting functions and hence the 
anomalous dimensions in perturbation theory. We have also verified the helicity sum rule 
explicitly to the first non-trivial order in perturbation theory. 

We also have compared and contrasted the expressions for internal orbital helicity in non- 
relativistic and light-front (relativistic) cases. Our calculation shows that the role played 
by particle masses in the internal orbital angular momentum in the non-relativistic case 
is replaced by the longitudinal momentum fraction in the relativistic case. Although four 
terms appear in the expression of L3 for individual particles in two body system, only the 
term proportional to the total internal L3 contributes due to transverse boost invariance of 
the multi-parton wave-function in light-front dynamics. We also note the occurrence of the 
longitudinal momentum fraction 22 (^i) multiplied by the total internal L 3 in the expressions 
of L 3 for particle one(two). This explains why one needs to take first moment with respect 
to x as well as (1 — a;) for the respective distributions in obtaining the helicity sum rule ||]. 

It should be emphasized that our explicit demonstration here, that the gauge-fixed light- 
front helicity operator is exactly equal to the naive canonical form, is facilitated by the fact 
that in light-front theory only transverse gauge fields are dynamical degrees of freedom. The 
conjugate momenta (color electric fields) are constrained variables in the theory. Thus we 
were able to show explicitly that the resulting gauge fixed operator is free of interactions. 
The question naturally arises as to whether this result is valid in other gauges also. Several 
years ago, in the context of magnetic monopole solutions, it has been shown |nj that in 
Yang-Mills-Higgs system, quantized in the axial gauge A s = using the Dirac procedure, 
the angular momentum operator constructed from manifestly gauge invariant symmetric 
energy momentum tensor differs from the canonical one only by surface terms. In the study 
of QCD in A 3 = gauge, it has been shown p0| that in the presence of surface terms, 
Poincare algebra holds only in the physical subspace. The situation in A = gauge or 
in covariant gauges where unphysical degrees of freedom are present is to be investigated. 
Another interesting problem to be studied is the helicity conservation in the topologically 
non-trivial sector of QCD and its implications, if any, for deep inelastic scattering. 
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SUMMARY AND CONCLUSION 

In this thesis, we have tried to put forth some of the aspects of light-front field theory 
through their successful application in a physical problem. The physical problem that we 
have chosen is DIS, which is strongly believed to have a convenient description in light- 
front language but seldom pursued consistently We have developed a light-front QCD 
Hamiltonian description of the DIS structure functions consistently working in the light- 
front co-ordinates and light-front gauge and using light-front QCD Hamiltonian as dictating 
the underlying dynamics. In retrospect, the importance of our work can be summarized 
from two different viewpoints as follows. 

As far as DIS is concerned, our approach is convenient in the sense that it closely fol- 
lows our physical intuitions as is the case in parton model where DIS has a simple partonic 
interpretation. Our approach differs from the usual covariant one in an essential way that 
it addresses directly the structure functions, which are experimental objects, instead of its 
moments which come naturally in the usual way. More importantly, our approach gives 
a complete description of the structure functions in the sense that it has the potential of 
incorporating the non-perturbative contents of the structure functions. We have shown in- 
troducing a new factorization scheme that the non-perturbative contents of the DIS structure 
functions can be obtained by solving the light-front bound state equations, which seems vi- 
able due to the ongoing research activity for solving QCD bound state problem in light-front 
using the similarity renormalization scheme. On the other hand, the perturbative contents 
can be extracted by calculating the dressed parton structure functions as we showed by 
explicitly working them out to the leading order in QCD coupling. In contrast, the usual 
approaches are so designed that one only parametrizes the nonperturbative information and 
the whole concentration is put on the Q 2 -evolution of the moments of the structure functions 
using perturbation theory. 

Simplicity of our approach becomes evident when we try to describe structure functions 
in the context of the nucleonic helicity by defining new structure functions in our approach. 
The ambiguity of gauge invariance and interaction dependence in defining various parts of 
the helicity operator for quarks and gluons in the usual way are absent in our gauge fixed 
theory and the well known helicity operator in light-front seems to provide the consistent 
physical information as we have shown explicitly. 

On the other hand, although the study in light-front field theory has been around for 
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quite some decades, there is still some doubt whether it produces the same result as the 
more familiar approach, since the formulation and the methods used in two cases are so 
different. Perturbative calculations of the dressed parton structure functions presented here 
show that, to the leading order in coupling, our approach yields equivalent results as in 
the usual way. This constitute a proof of equivalence of light-front field theory and the 
more familiar equal time field theory, even though the differences in formulation and the 
methods used in the two cases make the equivalence not so apparent. Our investigation in 
the context of coupling constant renormalization in light-front QCD hamiltonian methods 
shows the importance of Gallilean boost symmetry in understanding the correctness of any 
higher order calculation using time(x + )-ordered old-fashioned perturbation theory. 

Our study, of course, leaves open several issues that needs further investigation. Firstly, 
to add more substance to our proof of equivalence between light-front field theory and the 
equal time one, the dressed parton structure function calculation can be extended to the 
next higher order. This may turn out not only to serve as a check between various methods 
but also to give us the opportunity to investigate the intricacies involved in the two-loop 
calculation using light-front perturbation theory.^ 



* We prefer not to put any epilogue simply because it is a dissertation of an on going research 
work. And, in general, the quest for knowledge goes on and on. 
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APPENDIX A: NOTATION, CONVENTIONS, AND USEFUL RELATIONS 



Light-Front variables are defined in terms of x M = (x°, x 3 , x 1 , x 2 ) as 



x ± =x°±x 3 , x ± = (x l ,x 2 ). 



Let us denote the four- vector x^ by 



Scalar product x^x^ is defined as 



1 + - l - + ± ± 
x-y = -x^y + -x y^ - x .y . 



The metric tensors are as follows 

fO 2 0\ 

2 

0-1 

\0 -1/ 

so that 



jif 



9p.v 



f Q \ ° °\ 

| 

0-10 

\0 -1/ 



(6-1) 



(6.2) 



(6.3) 



(6.4) 



rp ry> <^ rv> rp 



Similarly for the partial derivatives: 



d + = 2<9_ = 2 



d 



«T = 2d, = 2 



dx~ 



d 

dx + 



Four-dimensional volume element: 



d 4 x = dx°d 2 x ± dx 3 = -dx + dx d 2 x A 

2 

Lorentz invariant volume element in momentum space: 

dk + d 2 k^ 



[d 3 k] 



2{27i) 3 k+' 



The step function 



(6.5) 

(6.6) 
(6.7) 

(6.8) 

(6.9) 



(x)= 0, x < 



1, x>0. 



(6.10) 
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The antisymmetric step function 

e(x) = 0(x)-0(-x). (6.11) 

where 5(x) is the Dirac delta function. 

| x | = x e(x). (6.13) 

We define the integral operators 

^p/OO = \j dy~e(x- - y-) f(y~), (6.14) 

(^r) 2 ^) = U dy ~ i x ~ ~ y~ i f(v~)- ( 6 - 15 ) 

Unless otherwise specified, we choose the Bjorken and Drell convention for the gamma 
matrices: 



where as are Pauli matrices 



1 / -i\ _ I 1 Q 

io' a2 -\i o ' as - -1 



(6.17) 



7 5 = , 7 o 7V7 3 = ^0 ^ (61g) 

a = 7°7. (6.19) 

7 ± =7°± 7 3 . (6.20) 

Projection operators are defined as 

A ± = ^7V = ^V = \(I±^). (6.21) 

such that 

A + + A" = I, (A±) 2 = A ± , A+A" = 0, (A ± ) t = A ± . (6.22) 

and they satisfy the following relations 

7 ± A ± = AV", 7° a± = A V, 7 5a± = A V- (6-23) 
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APPENDIX B: AN EXPANSION OF T^ IN INVERSE POWER OF 
LIGHT-FRONT ENERGY OF THE VIRTUAL PHOTON 

In this appendix we show how the virtual photon hadron forward Compton scattering 
can be obtained as an expansion in the inverse power of light front energy of the virtual 
photon q~ = q° — q 3 . This is the light front version of what is known as Bj or ken- Johnson- 
Low limit. Explicitly, the forward virtual photon hadron Compton scattering amplitude is 
given by, 

T"" = i J d^e iq< {PS\T{J^)J v ^))\PS) 

= i [d 4 t(^zd-e i <'<)(ps\e(t)[J ti (Z),r(o))}\PS) (6.1) 

Here, in the second line we have changed the time ordered product into a commutator and 
introduced a derivative operation d~ = 2-S+ on the exponential which does not alter the 
expression at all. Now, doing a partial integration we get, 



T^ = — [d*£e iq *d-(PS\0(£ + )[J>'(£),J' / {0)]\PS) 
a J 

>£ e ^(PS|2<$(^)[J^(£),J y (0)]|PS'} 



or 

1 

qr 
-— [ d i & iq *(ps\e(z + )[d-ji i (£), r(o)]\PS) (6.2) 

q J 
First term is integrated using the delta function and the fact that (i 4 £ = ^d^ + d^~d 2 ^ ± and 
with the second term we follow exactly the same procedure applied earlier, to obtain the 
following result. 

T» v = -— J dCd 2 ee iq< (PS\[J^O,J u ms+=o\ p S) 

-— fd^e iq ^{PS\26^ + )[idJ"{^,r{0))\PS} 
q J 

_J_ f d^e iq -HPS\e(c + )[d-(id-j^)),r(o)]\ps) (6.3) 

q J 
Notice that the first line of the RHS contains the equal-x + commutator and constitutes 
the first term in the BJL expansion, whereas second and third line ressemble with that of 
eq. (|6.2|) . Thus, iterating the above procedure we finally get an expansion of T^ u in the 
inverse power of q~ as follows. 

1 \n+l 
n=0 "9 

The above expansion shows that the time-ordered matrix element can be expanded in terms 
of an infinite series of equal light-front time commutators. 
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00 1 . . 
T " V = - E (-F / dC<Pt±e i9< (PS\W[) n J' t {t), r(0)y=o\PS) , (6.4) 



APPENDIX C: DETAILS OF ONE LOOP CALCULATIONS 

In this appendix, we discuss how the vertex functions (VF) and energy denominators 
(ED) in the light-front Hamiltonian perturbation theory are used in one loop calculations 
providing necessary details. In particular, we shall work out the VF and ED relevant in 
the context of F 2 calculation and then show various summation procedure which may be 
necessary in other cases as well. 

Recall that the dressed quark state | Pcr) q truncated at the two particle level, is given 
by 

\Pv) q = \[M q ib\P,a)\Q) 

/ i = = / i 2 = MP, v I K <n; k 2 , A 2 



x p(2-KfP+5\P -h- k 2 )tf(h, a x )a\k 2 , A 2 ) | 0)1, (6.1) 



which satisfies the light-front version of the Schroedinger equation 

P\ + M 2 
~P 



Pqcd\P°)= ^z l^>. (6-2) 



where Pqcd = Pa + V with P and V being the free and interaction parts of the LFQCD 
Hamiltonian respectively. Introduce the Jacobi momenta (x,-,kj~) 

so that 

5> = 1,£X = 0. (6.4) 

i i 

Also, introduce the boost invariant amplitude $2 as 

v^Mkt^i)=$2(x l ,^). (6.5) 

Now substituting the dressed quark state eq.(|6.1|) into eg. (|6T2|) and taking a projection on a 
bare one quark- one gluon state | 2) = tf(pi, ai)a){p 2) X 2 ) | 0) we get, 

(2 I [ ^p^ ~ Pol I P*) = (2 I V qqg I Pa) (6.6) 

Notice that the operator on the LHS being free of interaction | 2) projects out ip 2 an d to the 
order g RHS gets contribution only from V qqg . Straight-forward calculation starting from 
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eq.(|6.6|) using the Fourier decomposition of the dynamical fields and using the standard 
commutation relation for dynamical quarks and gluon fields, we get, 

■M 2 + (P 1 ) 2 m 2 + {pi) 2 (p^) 2 - 



9 



2(2tt) 



rriQ 



Pi 
1 



V~2 



ip 2 (P,o- I Pi,o"i;p 2j A 2 ) 



X.O\ 



pt_ °_ -Pi 

+ 

?2 



im 



a^.P J 



im 



-a 



o" 



Pi 



P^ 



X*-{£\ 2/ 



(6.7) 



where M and m are the masses of the dressed quark and bare quark respectively. The factor 
multiplying the ip 2 is the known as the ED and the RHS constitute the relevant VF here. 
We rewrite the above equation in terms of Jacobi momenta (pf = X{P + , nj- = pj- + XiP L ) 
and the wave-functions $j which are functions of Jacobi momenta. Using the notation 
x = Xi, K\ = k and using the facts x\ + x 2 = 1, n\ + k 2 = 0, we have 

1 



®2<n,\ 2 ( x > K ;i -»,-«" 



M 2 - 


m a +(«- L ) a 

X 


1-x 



X 



2(2tt) 3 y/l-x 



K.<T\ 



.K — im 



-a 



a im 



x 



x 



xAeiS- (6-8) 



Now, we have mostly used $2 in the massless cases M = m = 0, where it simplifies to 

g m „ 1 x(l — x* 



K,^,^) 



rjia 



2(2tt) 3 y/l-x (k 



±\2 A °"i 



1 — X 



+ Or"' 

X 



X*-{ e \ 2 , 



(6.9) 



The above example shows how easily one can obtain various two-particle wave-functions in 
terms of VF and ED. 

Now, to calculate various structure functions with a dressed quark target, we encountered 
this wave-function $ 2 where we needed to perform different summation over helicities. For 
example, 

1 



for 



F Vl)(„ ^>2> 



2q 



{x,Q< 



for Aq(x, Q 2 ) 
for Ag(l-x,Q 2 ) 



-V V I $ ff x 

o C t / ' I Ti,A2 
O cri,A2 

E^ I rT.cr |2 

0I5A2 

E a 2 I ^ ljA2 1 2 



(6.10) 
(6.11) 
(6.12) 



Here we provide the necessary details how these are actually worked out after calculating 
$2 explicitly for various helicity configurations. Here, we denote a = +| as (f) and so 
on which are self-indicating. To perform the summation over <j\ in the above expressions, 
we notice that in f 2 , 2 x 2 matrix in helicity space sandwiched between x and \ a is a 
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diagonal matrix. Thus, the matrix element will be zero unless o\ = or irrespective of the 
gluon polarization e\ 2 . That implies, taking o = +| 



$; t (x,q 2 )=o 



(6.13) 



For convenience, we change the variable as {k,\ — » |k|cos0, K2 — > |fi:|sin0}, so that k\ = \k\ 2 
and k- et — K\ — IK2 — \k\ e - *^, and so on. For a% = o = +| and A2 = 1, $2 can be written 



as 



TT 



•)xj 



V2 



:---)v / 2 



Similarly, for A2 = — 1, we have, 



tr 2(«-ef) («-<r)(ff-e|) 



1 — x 



X\ 



2\k\ e 



1 — x 



+ 



XT 
|«| e~^ + cr 3 |ft| e" 



a; 



XT 
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k| e" 


-«</> 


V2 


x(l- 

|k| e" 


a?) 



2x + (l-x)(l + l) 



x(l — x) 



(6.14) 



*n- 



■)x\ 



2(« • el) , (K-a)(a-e|) 



v/2 



1 — x 

tr 2 N e+ ^ 



X 



XT 



3 +* 



+ 



/d e" 1 " v + (7°k e 



+40 



x 



XT 
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|«| 


e +i<P 


v/2 


x(l 


— x) 



2x + (l-x)(l- 1) 



(6.15) 



x(l — x) 

Here we have used x' s an d polarization vector e's as are given earlier and the fact that 
Xt cr3 XT = 1 an d Xt°" 3 Xt = ~~ 1) while XaX/3 = 4*/3- The omitted common factor in the above 
expressions is (••■) = ; -* -75 / X n 2 ■ 



Now, from eq.( |6.13|) , eq. (|6.14j) and eq. (|6.15|) , it is straight forward to perform the required 

2\ 



summations as follows. For F 2q (x,Q 



1 1 



5 -5?< 



2k| 2 2|k| 2 x 2 

+ 



x 2 (l — x) 2 x 2 (l — x)' 



3 c^{\)- 

1 — X v K, J 



2(2tt) 



(6.16) 

(6.17) 
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Note that here in the unpolarized case averaging over initial spin a = ±| does not matter, 
since both are same. Next, in the polarized case of Aq(x, Q 2 ), 



V o-i I $ T , I 2 = +-( I $1, I 2 + I $!, I 2 \ 



2 

0"1,A2 



1 9 2 r 



2 2(2^-1^(4)- (6 - 18) 

(6.19) 

Since the contribution from <j\ =| is zero, the result in the last equation is simply helicity +| 
times the expression in unpolarized case( eq. fl6.16j )). It trivially follows that if we considered 



the dressed quark to be in helicity —\ state, i.e., o —\ the result would have been negative 
of what we have in eq.( |6.18| ). On the other hand, for Ag(l — x, Q 2 ), 



ea 2 i$ui 2 = {(+i)i$j t i 2 + (-i)i$;j 2 } 

2kl 2 , ^ 2|«IV 



o"i,A 



/ \l{ ( i\ II / i\ II 1 

*- x 2 (l — x) 2 x 2 (l — x) 2 > 



U ~ rCV(l + x)(A). (6 .20) 



2(2tt) 3 rv ^k 



APPENDIX D: MANIFEST BOOST SYMMETRY OF ENERGY 
DENOMINATORS AND VERTICES 

In this appendix we verify the Gallilean boost invariance of vertices and energy differences 
that occur in light-front time-ordered loop diagrams. First consider the canonical vertex 
given in eq.(5.4). Let P + and P -1 denote total longitudinal and transverse momentum in 
the problem. We introduce the momentum fractions X{ and the relative transverse momenta 

*i b y 

p + = x p+ ; p ± = K ± + xP X , q + = (l-x)P + , q 1 - = -K^ + (l-x)P ± . (6.1) 

The longitudinal momentum fractions x,i and the relative transverse momenta nf obey the 
constraints Yl x i = 1 an d J2 K j~ — 0. The canonical vertex takes the form 

Vi = 9 T a yfr xl [2-^- + a -^± a ± + im (l - -)} Xsi . (ef )*. (6.2) 

In terms of the internal momenta, the boost invariance of the quark-gluon vertex is clearly 
manifest . 
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Next consider loop diagrams. As an example we consider the diagram shown in Fig. 
2(a). Parameterize the single particle momenta in terms of the internal momenta as follows. 



yP" 






/"» ■-} 



yP ± , kf = (l- y)P + , fcf = -4 + (1 - y)P J 



Then 



(6.3) 



k+ = kf - q + = (x - y)P + , k 1 - = k± - q 1 - = k± - 4 + ( x ~ y) pi 
The energy difference appearing in the two energy denominators are, then, 
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(6.4) 



(6.5) 
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(6.6) 



y i-y 

Thus the vertices and energy denominators appearing in Fig. 2(a) are manifestly invariant 
under the Gallilean boosts in the transverse plane and this is a general property of any x + 
ordered diagram in light-front perturbation theory. 
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